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The induced pseudoscalar coupling g p is the least well known of the weak coupling constants of the 
proton's charged— current interaction. Its size is dictated by chiral symmetry arguments, and its 
measurement represents an important test of quantum chromo dynamics at low energies. During 
the past decade a large body of new data relevant to the coupling g p has been accumulated. This 
data includes measurements of radiative and non radiative muon capture on targets ranging from 
hydrogen and few— nucleon systems to complex nuclei. Herein the authors review the theoretical 
underpinnings of g p , the experimental studies of g p , and the procedures and uncertainties in 
extracting the coupling from data. Current puzzles are highlighted and future opportunities are 
discussed. 
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I. INTRODUCTION 

It is well known that the weak interactions of leptons 
are governed by a current-current interaction, where the 
currents are given by a simple V — A form, 7 M (1 — 75). 
When hadrons are involved, the interaction is still 
current-current and still V — A but the individual vector 
and axial vector currents become more complicated, pick- 
ing up both form factors and new structures involving 
the momentum transfer q. The couplings associated with 
these new structures are the so called 'induced' couplings 
which have been a topic of investigation, both theoretical 
and experimental, for a long time. We are concerned in 
this review with one of these, the induced pseudoscalar 
coupling, g p . 

The most general weak current, actually the matrix 
element of that current, for say a neutron and proton 
can be written as 

— ft u 1 l 9m au 9s u 
Un( + 9vT + n a Q» H 1 

- 5a7 M 75 - ^V75 - ^^<Z„75K (1) 



where t he notation for gamma m atrices, spinors, etc. is 
that of Bjorken and Drcll (1964) and m M and are 
respectively the muon and nucleon masses. The momen- 
tum transfer q = p n — p p , with p p , p n respectively the 
proton and neutron momenta. This most general form 
for the current contains six coupling 'constants', which 
are actually functions of q 2 , namely g v and g a , the vector 
and axial vector couplings, g mi the weak magnetism cou- 
pling, <? p , the induced pseudoscalar coupling, and g s and 
g t , the second class induced scalar and induced tensor 
couplings. 1 

Of these six, the second class terms g s and g t trans- 
form differently than the others under G-parity, which 
is a combination of charge conjugation and a rotation 
in isospin space. In the standard model they are gener- 
ated only via quark mass differences or electromagnetic 
effects and are thus predicted to be small ( jGrenacs"^ 



1985 ; Shiomi , 1996] ). Thus although the experimental 



evidence is not un equivocal ( Minamisono et ai, 2001 
Wilkinson , 2000a| jb|), they are normally assumed to be 
zero and will be ignored in this review. The vector terms 



1 There is an unfortunate confusion about the signs of the axial 
couplings which has arisen as conventions have changed over the 
years. In very early calculations in a different metric, the signs 
were such that a n was positive. Wit h the widespread use of the 



metric of Bjorken and Drell (19^4 ) it became conventional to 
write the weak current like Ea. f yTonlv with all signs positive . 
This was the co nvention used in Reder and Fearing (1987, 1989); 



Fearing (1980) and most other modern papers, and implies that 
g a and g p are negative. With the advent of chiral perturbation 
theory, for which g a is normally taken to be positive, the con- 
vention changed again. We have adopted this latter convention. 
This means that the axial current of Eq. (hi) has an explicit over- 
all minus sign and that g a and g p are positive numbers. 
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g v and g m are related by the Conserved Vector Current 
(CVC) theory to the isovector electromagnetic current. 
Thus their static values are related to the charge and 
anomalous magnetic moments of the nucleons and their 
q 2 dependence is given by measurements of the isovec- 
tor electromagnetic form factors of the proton and neu- 
tro n. The axial vector coupling g n can be determined 



II. THEORETICAL PREDICTIONS FOR g p 

A. PCAC, g p , and the Goldberger-Treiman relation 

Historically the measurement of g p has been considered 
interesting and important because there is a very definite 



prediction for its value deri ved many years ago (Gold- 



pre cisely from neutron beta decay and its q 2 dependence berger and Trciman , 1958|) based on the very funda- 



from neutrino scatterin g (Ahrens et ai, 1988) or from 
pion electroproduction (Bernard et ai, 1992; Choi et 



1991 Del Guerra et al, 1976; Esaulov et al, 1978) 



This leaves the induced pseudoscalar coupling, g p , 
which is by far the least well known. In principle it can 
be predicted using the Partially Conserved Axial Current 
(PCAC) theory vi a a relation known as the Goldb erger- 
Treiman relation ( Goldberger and Trciman , 1958| ). The 
details of this derivation will be discussed in the next sec- 
tion. Physically the dominant diagram is one in which 
the nucleon emits a pion which propagates and then cou- 
ples to the fj, — v or e — v vertex. The coupling g p thus 
contains a pion pole. In practice however it has been 
difficult to verify this relation and there are some situa- 
tions, notably radiative muon capture in liquid hydrogen, 
where experiment and theory do not agree. 

Thus the main purpose of this review is to discuss in 
detail the information available, both theoretical and ex- 
perimental, regarding the value of g p . In Sec. [fi] we dis- 
cuss the theoretical prediction for g p arising from PCAC 
and from the more modern approaches using chiral per- 
turbation theory. Sections |n| and |l^ describe in general 
some of the sources of experimental information and some 
of t he atomic and molecular physics required to undcr- 



mental notion of the Partially Conserved Axial Current 
(PC AC). Detailed derivation s are given in textbooks, 
e.g. Bjorken and Drcll (1964), so we simply outline the 
basic ideas here. The underlying assumption is that the 
divergence of the axial current is proportional to the pion 
field. Applying this idea to the divergence of the axial 
current of Eq. (Q), taking matrix elements, and evaluat- 
ing at four-momentum transfer q 2 — 0, gives the relation 
known as the Goldberger-Treiman relation: 



(2) 



where g^NN is the pion nucleon coupling constant, = 
92.4 ± 0.3 MeV is the pion decay constant, and uin is 
the nucleon mass. This equation is rather well satisfied 
when one uses modern values of g n NN- A measure of the 
difference, which is known as the Goldberger-Treiman 
discrepancy, is given by the equation 



1 



m N g a (0) 

girNNimDFn' 



(3) 



The value of Agt was about 6% using the older (larger) 
value of g^NN = 13.4, and there were a number of pa- 
pers discusshig_j30j^ 

e.g. Coon and Scadron (1981, 1990); Jones and Scadron 



star jd muon processes in hydrogen or deuterium. The 
remainder of the review is divided into three main parts 
which discuss the three main areas which have led to in- 



(1975) and references cited therein. With the newer and 

13.05 ± 0.08 



somewhat smaller value of 



forn |tation on g v . The first of these, Sees. V ■ VIII , deals 
wit h ordinary, non radiative, muon capture fOMC) and 
radi ktive muon capture (RMC) on the proton, deutcron, 



g^NN\m v ) 



(Arndt et ai, 1995; de Swart et ai, 1997; stoks et al. 



and 3 He. Here the experiments are difficult, but the re- 
sults presumedly are not obscured by unknown details of 
nuclear structure. In the next major part, Sec. DC, exclu- 
sive OMC on complex nuclei is discussed, with particular 
emphasis on spin observables in partial transitions. Here 
the sensitivity to g p is high, but the experiments are of- 
ten difficult and their interpretation is model dependent. 
Then in the final part, Sec. [x], inclusive RMC on com- 
plex nuclei is examined. RMC is quite sensitive to g p 
and nowadays measurements of inclusive RMC on Z > 2 
nuclei are relatively straightforward. However model de- 
pendence in the theory makes the extraction of a unique 
value of g p very difficult. Finally in Sec. XI we will sum- 
marize the current situation with regard to g p and make 

somp siiggpstioms for fiirt/hpr work 



1993), and updated values of g a and F n the di screpancy 
is now 2% or less. See e.g. poity et al. (199E ); Nasrallah 
(2000 ) and references cited therein. Thus at q 2 — the 



Goldberger-Treiman relation is quite well satisfied. 

Now consider the matrix element of the divergence 
of the axial current for non-zero q 2 . This leads, using 
Eq. (^), to an expression for g p given by 



9p(q 2 ) 



m 2 — q 2 



9a(0) 



(4) 



with m-n- the charged pion mass. This is known as 
the Goldberger-Treiman expression for g p . Observe the 



explicit presence of a pion pole. At q 



-0. 



Rela ted review articl e s include those of Mukhopadhyay 
(19l j7|),|Grenacs (1985|), pmitro and Truol (1987| ), |Mcas- 
day (2001|), andjBernard et al. (2001a|). 



which is the relevant momentum transfer for muon cap- 
ture on the proton, this formula gives g p {— 0.88m^) = 
6.77g a — 8.58, where we have use d the latest value of 
g a (0) = 1.2670 ± 0.0035 from the |Particlc Data Group 
(2000| ) and taken for rajy the average of neutron and 
proton masses. This is what we shall refer to in sub- 
sequent sections as the lowest order PCAC value g p GAC '. 
At q 2 = the result is g P {-m 2 ) = 6A7g a = 8.20. 
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Somewhat later the first order correction to this, pro- 
portional to the derivative of g a , was derived using cur- 
rent algebra techniques (Adler and Dothan , 1966). Nu- 
merically however the correction is rather small, as de- 
scribed in the next section. 



B. ChPT derivations of g p 

In the years since the original derivations of the 
Goldberger-Treiman relation there have been major ad- 
vances in our understanding of the way to include chiral 
symmetry in such analyses, particularly in the framework 
of what is known as chiral perturbation theory (ChPT), 
or, when nucleons are involved, heavy baryon chiral per- 
turbation theory (HBChPT). This approach provides a 
way of incorporating the symmetries of QCD into a sys- 
tematic low energy expansion, where the expansion pa- 
rameter is something of the order m^lm.M. Thus one can 



If one wants to express g p in terms of g a , as is conven- 
tional, rather than f^g-nNNi some simple manipulation 
of Eq. (||) using Eq. (||) gives 

, a (0)(l + g)- m ^f (0) ^ (9) 



. 2m fli m N 
[mi - q z ) 



where (1 + e) = [sWiViv(g 2 )/ fl , 7rJVJv(0)]. Note that (1 + e) = 
1 if one neglects the q 2 dependence of g^NN, which we will 
do in this review. This gives 5 P (—0.88m^) = 8.58—0.45 = 
8.13 where the slight difference from the result of Eq. (||) 
originates in the difference between left and right hand 
sides of Eq. (|^) when experimental values are used. It is 
this latter value, 8.13, which we will take as g p CAC when 
the constant term is included. 

It is important to note that the specific results in terms 
of the LEC's of ChPT depend on the specific choice of 
starting Lagrangian and the details of the calculation 
The expressions of Eqs. (||), (0), and (||) come from Fear 



repr bduce all the old current algebra results, but more ing et al. (1997|), but similar re sults were subsequentl y 



importantly calculate in a systematic way the corrections 
to these results. 



This approach was applied ( Bernard et al, 1994 ) to 
obtain for g p the result 



9p(Q 2 ) = 



2m ll g vNN (q 2 )F^ 



^ga(0)m^m N r 2 A , (5) 



where r\ is the axial radius of the nucleon, defined in the 
usual way via g a (q 2 ) = g a (0)(l + q 2 r 2 A /6 + 0{q 4 )). This 



is e ssentially the re sult obt ained much earlier by Adler 
and |Dothan (1966j ) and by |Wolfcnstcin (1970| ), but the 
systematic approach allows us to be confident that the 
corrections are of higher order. 

Using the value of the axial radius r\ = 0.42 ± 0.04 fm 



obt ained fr om anti neutrino- nucleon scattering (Ahrcns 
et a 



0.45 



1988J) , this formula leads to g p (- 



-0.88m 2 ) 



8.70- 



8.25, so the correction term is indeed rather small. 
An alternative approach, still within HBChPT, is sim- 
ply to calculate the amplitude for muon capture, as was 
done by Fearing et al. (1997), and identify g p by compar- 



ing with Eq. (111). This gives, up to corrections of 0(p 



9 P (q 2 ) 



2m M mjv 
(ml - q 2 ) 



rf) 



(4ttF t ) 



r(26i 9 -6 23 ) 



(6) 



where big and 623 are low energy constants (LEC's) of 
the basic Lagrangian. At first glance this appears to be 
different than the result above, but can be put in the form 
of Eq. (|5|) by noting that in the same approach g n N n and 
the axial radius squared r 2 A are given by 



2 m N 
TNN(q ) = 



9a(0) 



8tt 2 F 2 



and 



r A = -6 



g a {0){^F^f 



(7) 



(8) 



obtained by |Bernard et al. (199q ) and |Ando et al. (2000[ ) 
However the expressions in terms of the physically mea- 
surable quantities, as given in Eqs. Q or (|^) are inde- 
pendent of the detailed conventions of the approach. 

Finally we can summarize this section by observing 
that theoretical prediction for g p as given in Eqs. (||) 
or (^|) and based on chiral symmetry seems to be quite 
robust. The original Goldberger-Treiman relation is un- 
derstood as the first term in an expansion and correc- 
tion terms have been evaluated and are understood via 
a ChPT calculation carried out through the first three 
orders. Thus a test of this prediction should be an im- 
portant test of our understanding of chiral symmetry and 
of low energy QCD. 



III. SOURCES OF INFORMATION ON g p 

At the simplest level g p appears as a phenomenologi- 
cal parameter in the fundamental definition of the weak 
nucleon current, Eq. (|l|), so it should be obtainable from 
any process which directly involves this current. This 
would include beta decay, muon capture, radiative muon 
capture, and in principle any of the crossed versions of 
these reactions, as for example processes initiated by neu- 
trinos. The g p term is proportional to the momentum 
transfer however, so in practice beta decay is not a use- 
ful source of information since the momentum transfer 
is so small. Neutrino processes are of course extremely 
difficult to measure. This leaves ordinary muon capture 
and radiative muon capture as the two main sources of 
information on g p . 

Thus for ordinary, i. e., non radiative, muon capture, 
g p is purely a phenomcnological parameter which appears 
in the most general weak current. PC AC is not needed 
and it is quite reasonable to simply fit to data, treating 
g p as a free parameter. This is what has typically been 
done, with the result looked upon as a test of the PCAC 
prediction. 
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At a somewhat deeper level the g p coupling is under- 
stood to arise from the diagram in which a pion, emitted 
from a nucleon, couples to the weak leptonic current. 
Since this pion-nucleon coupling is a component of the 
axial current, other processes involving this current such 
as pion electroproduction for example, also in principle 
allow one to obtain information on g p . The interpretation 
of information from such processes must be somewhat 
different, however, than that obtained from processes 
like muon capture which contain the phenomenological 
weak current explicitly. For processes like electroproduc- 
tion the direct information available is really information 
about the ttNN vertex. Thus the connection to g p is only 
via the theoretical infrastructure of chiral symmetry and 
of PCAC and the interpretation that g p originates in a 
pion exchange diagram. If PCAC were in fact wrong, 
the whole connection would break down and there would 
be no convincing physical reason why fitting data using 
for g p a multiple of the PCAC expression containing the 
pion pole should work. 

It is interesting to note that RMC is somewhere in 
between the situations corresponding to OMC and pion 
electroproduction. The dominant diagrams contain the 
phenomenological weak current directly, albeit with one 
leg off shell, but there is also a diagram which explicitly 
contains pion exchange. 

Finally there is a third level where the weak current 
and one pion exchange do not appear directly but where, 
in the context of HBChPT, some of the LEC's needed for 
g p do appear. For such processes it is, at least in princi- 
ple, possible to determine those LEC's, and thus deter- 
mine Op, at least indirectly, via an equation analogous to 
Eq. (|6j). Again such information must be interpreted in 
a context which accepts the validity of PCAC. 

For all of these processes it may be that measurable 
quantities such as correlations relative to some of the 
particle spins, or capture from hyperfme states, or cap- 
ture to or from a specific nuclear state, or some such more 
detailed observable may provide more information than 
just the overall rate, so these should be considered. 

We will begin by looking at these various processes, 
starting with the simplest, muon capture on the proton, 
and working up through reactions on nuclei to see what 
has been learned and what potentially could be learned 
with regard to g p . 



system. Thus a detailed knowledge of this muon chem- 
istry is required in order to determine g p from H2/D2 
experiments and it is appropriate to discuss this chem- 
istry before considering the capture process. Below we 
denote the two hyperfme states of the muonic atom by 
F± = J± 1/2 where 1=1/2 for the proton and 1=1 
for the deuteron, so that F is the total spin of the atom. 

Before we discuss the details we make a few general 
comments. One important aspect of muon chemistry is 
/z-atom scattering from surrounding molecules which re- 
sults in the hyperfme depopulation of the upper F-state 
into the lower F-state. Another important aspect is col- 
lisional formation of muonic molecules which results in 
additional arrangements of /i-nuclcus spin states. Also 
muon recycling from molecular states to atomic states 
via /i-catalyzed fusion occurs for d/id molecules and pfid 
molecules. How these effects unfold for muons in hydro- 
gen and deuterium, as a function of the density, is the 
focus of ou r discu ssio n in th is section. 

In Sees. IV. A and [V.B, respectively, we discuss the 



chemistry of / ip ato ms in pure H2 and /id atoms in pure 
D2. In Sec. IV. C we describe the chemistry of muons 
in H2/D2 mixtures. Related r e view articles are those o f 
Bracci and Fiorentini (1982|); Brcunlich et oL (1989); 



Froelich (1992| ); |Ponomarev (1973| ) 



A. Muons in pure hydrogen 

To assist the reader a simplified diagram of muon 
chemistry in pure H2 is given in Fig. [l| The figure shows 
the F + and F_ states of the /ip atom, the ortho (1=1) 
and para (1=0) states of the p/zp molecule, and relevant 
atomic and molecular transitions. 

The /ip atom is initially formed in a highly excited 
state with a principal quantum number n ~ 14 and a 
kinetic energy ~ 1 eV. The excited atom then rapidly 
de-excites via combinations of Auger emission 

(np) n + e —> (fip) n > + e, (10) 

radiative decay 

(fj,p) n {p.p)n' + 7) ( n ) 
and Coulomb de-excitation 



IV. MUON CHEMISTRY AND THE INITIAL SPIN 
STATE 

When a negative muon is stopped in matter a muonic 
atom is formed. Unfortunately in hydrogen and deu- 
terium such atoms undergo a complicated sequence of 
chemical processes, changing the spin-state populations 
with time, as the muons eventually reach the level from 
which they are captured. As discussed elsewhere, the 
capture rates for muonic hydrogen and muonic deuterium 
are strongly dependent on the spin state of the /z-nucleus 



{p.p)n + p^ (Hp) n > +P- 



(12) 



Note that in Auger emission and radiative decay the p,p 
recoil has a relatively small kinetic energy, since the e or 7 
carry the released energy, while in Coulomb de-excitation 
the fip recoil has a relatively large kinetic energy, since 
the fip and p share the released energy. Consequently 
when formed the ground state atom has kinetic energies 
of typically about 1 eV but occasionally up to 100 eV. 
For recent experiment s on energy distr ibuti ons of ground 
state np/wp atoms see Bigg et al. (1996) and Schottmiillcr 
et al. (1999| ) . Note that at formation of the ground state 
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F=1 



107 eV 



1=1 



, ortho 



_eV_ 



0.18 eV 



para 

PMP 
S=3/2 

S=1/2 



different make-ups in terms of \ip components with par- 
allel spins (i.e. F + ) and anti-parallel spins (i.e. FJ). 
Specifically the para-state is 3:1 triplet-to-singlet, the 
5=1/2 ortho-states are 1:3 triplet-to-singlet, and the 
5 = 3/2 ortho-states are pure triplet. For further details 
see Bakalov et al. (1982| ). 

The ppp molecules are formed by Auger emission 



0.14 eV 



H 



ppp- 



(15) 



Calculation s of the rates for the process have b een per- 



forme d by Faifman and Men'shikov (1999| ) , Faifman 



(1989|), iPonomarev and Faifman (1976Q , and IZcl'dovicb 
and Gcrshtcin (1959| ). Formation of the ortho-state in 



FIG. 1 A simplified diagram of the various atomic and molec- 
ular states and transitions relevant to muon capture in pure 
H 2 . 



ato m the hyperfme states are statistically populated, i.e. 



3:1 or F+-.F-. 



This 'hot' ground state atom is rapidly thermalized by 
elastic scattering, 



and spin- -nip collisions, 



(mp)ti +p ^ (mp)tt +p> 



(14) 



fron t the atomic nuclei of the neighboring moloculoo. 

Here TT denotes the triplet state and TJ, denotes the sin- 
glet state. Once the pp kinetic energy falls below the 
0.18 eV pp hyperfme splitting the singlet-to-triplet tran- 
sitions are energetically forbidden and triplet-to-singlet 
transitions depopulate the upper F + state. The resulting 
F + state lifetime is about 0.1 ns in liquid H2 and about 
10 ns in 10 bar H2 gas. The short lifetime arises from the 
large pp scattering cross sections due to a near-threshold 
fip + p virtual state. For further details of experimental 



stu dies of u p scattering in H2 environments see Abbott 
et a . (1997 ) and references therein. 

At sufficient densities, i.e. at pressures exceeding 10 
bar, the formation of ppp molecules is important. The 
molecule (see Fig. ^) comprises a para-molecular ground 
state, with total orbital angular momentum I — and nu- 
clear spin 1 = 0, and an ortho-molecular excited state, 
with total orbital angular momentum I = 1 and nuclear 
spin 1=1. Note spin-spin and spin-orbit interactions 
produce a five-fold splitting of the ortho-molecular state 
with two 5=1/2 sub-states and three 5 = 3/2 sub- 
states, where 5 is the total spin angular momentum of 
the pfj,p complex. Importantly the different states have 



volves an El transition with a predicted rate A"^ 10 ~ 
n/n x 1.8 x 10 6 s _1 and formation of the para-state in- 
volves an EO transition with a predicted rate A^™ ~ 
n/n x 0.75 x 10 4 ,s _1 , where n/n is the H 2 target number 
density normalized to the liquid H2 number density. Re- 
cent measurements of the summed, i.e. ortho and para, 
rate are typically 30% grea ter than the calculated rate. 
See Mulhauscr et al. (199E) for further details. 



Note that the El transition feeding the ortho- 
molecular state populates only the 5=1/2 s ub-states 
(i.e. y ieldi ng a 1:3 ratio of F + : F_ spin-states). Weinberg 
( 1960 ) and Ando et al. (200C| ) have discussed the possi- 
ble mixing of the 5 = 1/2, 3/2 levels which would lead to 
changes in the 1:3 triplet-to-singlet ratio for th e ortho 
molecu le. However the available calcula tions of Halpcrn 



(13) (1964a|M ; |Wessel and Phillipson (1964|) and |Bakalov et 
al. (1982|) have suggested such effects are negligible.^ 



At first glance the A J = selection rule for El tran- 
sitions forbids decay of the ortho excited s tate to the 
para grou nd state. However as discussed by Bakalov et 
al. (1982 ), via the small components of the relativistic 
wave functions, the ortho state contains 7 = admix- 
tures and the para state contains 1=1 admixtures. 
Therefore ortho-to-para El transitions occur via cross- 
combinations of the small components and the large com- 
ponents of the molecular wave functions. Note the rate 
A op for this Auger process is a function of the electron en- 



vironment of the ppp molecule. Bakalov et al. (1982) ob- 



tained A op = 7.1±1.2xl0 4 s 1 assuming an environment 
consisting of 75% [{ppp) + 2p2e} + and 25% [(pfip) + e]. 
These pr oportions are a consequen ce of the Hirshfcldcr 
reaction ( Hirshfcldcr et al., 193q ) involving ppp com- 
plexes and H2 molecules as disc ussed by Faifman (1989); 
Faifman and Men'shikov (1999 ). Note that the only pub- 
lishe d experimental value of A ov = (4.1 ± 1.4) x 10 4 s _1 
from Bardin (1982 ) ; Bardin et al. (1981b ) is in marginal 
disagreement with the calculation at the 2a level. 



2 In addition the relative rates for fj, capture in H2 gas (mostly 
singlet-atom capture) and H2 liquid (mostly ortho-molecule cap- 
ture) are consistent with a 1:3 tri plet:s inglet make-up of the 
ortho-molecule. For details see Sec. |V.D| 
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In summary, in H2 gas at pressures 0.1 < P < 10 bar, 
where the F — 1 atoms disappear very quickly and the 
Pfip molecules form very slowly, the capture process is 
essentially dominated by singlet atoms. However in liquid 
H2 the molecular formation rate A PW and ortho-to-para 
transition rate A op are important. Here the rate is a 
superposition of singlet atomic capture, ortho-molecule 
capture and para-molecule capture and depends on A pflp , 
A op , and the measurement time window. 



B. Muons in pure deuterium 

The atomic capture and cascade processes for muons in 
pure D2 and pure H2 are very similar. Most importantly 
for muons in D2 the ground state fid atoms are rapidly 
formed in a statistical mixture of the hyperfine states, 
i.e. 2:1 for 

The fid atoms are then thermalized by elastic scatter- 



ing 

Od)u + d ■ 
and spin-flip collisions 



(fid) 



(M)tJ. + d <-> (/irf)tT + d 



(16) 



(17) 



on surrounding nuclei. When the /id kinetic energy falls 
below the 0.043 eV hyperfine splitting the spin-flip col- 
lisions then depopulate the higher lying F + state. How- 
ever the cross sections are considerably smaller for fid + d 
than fip + p and consequently the F+ lifetime is consid- 
erably longer in deuterium than hydrogen. For example 
in liquid D2 at 34 K, the hyperfine depopulation rate is 
A = 4 2.6 x l O 6 ^ 1 . For further details see |Kammel et al. 
(19^|, |l983| ) 



An interesting feature of [id chemistry is resonant for- 
mation of dfid molecules. For example in liquid D2, while 
dfid formation from F = 1/2 fid atoms involves a non- 
resonant Auger process, the dfid formation from F = 3/2 
[id atoms involves a resonant excitation process, i.e. one 
where the dfid binding energy is absorbed by D2 vibro- 
rotational modes. In liquid D2 at 34 K the effective rates 
are A x J^ d ~ 5x 10 4 s _1 for non-resonant, i.e. doublet, for- 



mation and A 



3/2 

dpid 



4 x 10 6 s 1 for resonant, i.e. quartet, 



formation. Note that the resonant form ation is tcmpcr- 
ature dependent. For further details see Brcunlich et al. 
(19%. 

When dfid molecules are formed, the two deuterons 
immediately fuse via /x-catalyzed fusion 

d+d -^ 3 He + n [B.R. ~ 55%] (18) 



d + d -^ 3 H + v [B.B, ~ 45%] (19) 



d/id molecules, independent of density and temperature. 
However the fi sticking probability in fi catalyzed fusion, 
13% in Eq. @ and 1% in Eq. @, is non-negligible. 
Consequently with increasing D2 target density an in- 
creasing /i 3 He capture background is unavoidable. 

In summary, for muons in pure D2 the capture is a 
superposition of the rates from the doublet atom and 
the quartet atom and fi capture from dfid molecules is 
completely negligible. Actually for the particular c ondi- 
tions of the fid experiments by Bardin et al. (1986) and 
Cargnclli et al. (1989) ) the deuterium c aptur e is almost 
entirely from doublet atoms. See Sec. Vd.Bl for details. 
However due to muon sticking in the dfid fusion process, 
a correction for backgrounds from fi 3 He capture is nec- 
essary in liquid D2. 



C. Muons in hydrogen-deuterium mixtures 

Our interest in hydrogen-deuterium mixtures is two- 
fold. First 'pure-H2' experiments and 'pure-D2' experi- 
ments must inevitably be concerned with contamination 
from other isotopes. Second some early experiments on 
deut erium capture use d hydrogen-deuterium mixtures, 



■JL. 



Wang et al. (1965| ) used 0.3% D 2 in H 2 liquid and 



Bcrtin et al. (1973j ) used 5% D 2 in H 2 gas. 
Muon transfer from fip atoms to fid atoms 



HP 



fid + p 



(20) 



occurs with an energy release of 135 eV and a transfer 



rate of n/n xcjx 1.7 x 10 



where a is the D2 



concentration in the H2 target a nd n/n the target num - 
ber density relative to liquid H2 (Adamczak et al., 1992). 
Consequently, in H2 liquid a 10~ 3 D2 concentration and 
in P > 10 bar H2 gas a 10 _2 D2 concentration, is suffi- 
cient to engineer the transfer in roughly 100 ns. 

Following transfer, the fid atom is thermalized via col- 
lisions with H2 molecules. An interesting feature of fid+p 
scattering is the Ramsauer-Townsend minimum at a ki- 
netic energy 1.6 eV. The tiny fid + p cross section means 
slow thermalization and large diffusion of fid atoms in H2 
gas. In addition the slow thermalization and small deu- 
terium concentration makes hyperfine depopulation via 
spin-flip collisions 



(W)u 



(fid) 



TT 



(21) 



very slow by comparison to fip atoms in pure H2 and fid 
atoms in pure D2. 

At sufficient densities the formation of pfid molecules 
occurs by Auger emission 



fid + H2 — > pfid ■ 
with a rate in l iquid H2 of A pfl d = 



H 



(22) 



5.6 x lOS" 1 (Pctitjean 



which quickly recycles the muon from the molecular 
states to the atomic states. Consequently for muons in 
pure D2 the deuterium capture is from fid atoms and not 



et al., 1990/91 ). The formation of molecules is important 



as (i) capture is consequently a superposition of fid re- 
actions and fip reactions, and (ii) the various pfid states 
have different decompositions into fid -F-states. Further 



8 



the parent distribution of fid atom F-states effects the 
resulting distribution of pfid molecule states, making the 
relative population of pfid states a complicated function 
of target density and deuterium concentration. 

Muon catalyzed fusion from pfid molecules occurs via 
both radiative reactions 



Pfid — > fi 3 He + 7 
and non-radiative reactions 

pfid — > fi + 3 -ffe. 



(23) 



(24) 



The radiative rates ( Petitjean et ai, 1990/91 ) are A1/2 = 
0.35 x 10 6 s^ 1 and A 3/2 = 0.11 x lO^" 1 and the non- 
radiative rates are A ± / 2 — 0.056 x 10 6 s _1 and A 3 / 2 = 0, 
where the subscripts 1/2, 3/2 denote the p-d spin states. 



Fortunately a lot is known about the couplings. The 
weak vector current is completely determined by the well 
established Conserved Vector Current theory which tells 
us that the weak vector current is simply an isospin 
rotation of the isovector electromagnetic current. In 
practice this means that g v (q 2 ) = F\(q 2 ) — F™(q 2 ) and 
9m(q 2 ) = K p FP{q 2 )-K n F?(q 2 ) where k p = 1.79285, n n = 
— 1.91304 are the proton and neutron anomalous mag- 
netic moments and where Ff(q 2 ), F{ l (q 2 ), F% (q 2 ), and 
F^iq 2 ) are the usual proton and neutron isovector elec- 
tromagnetic form factors. The coupling g a is well deter- 
mined from neutron beta decay, g a (0) — 1.2670 ± 0.0035 
( Particle Data Group , 2000] ). The momentum depen- 
dence is known, at least at low momentum transfers, 
from neutrino scatterin g (Ahrcns et ai, 1988) or from 
pion electroproduction ( I Choi et ai, 1993 ; Del Guerra ei 



The | slow rates make fi capture from pfid molecules an a l : 1976| ; |Esaulov et ai, 1978[ ). For a long time there was 



important contribution at high target densities. Fur- 
ther the 100% sticking probability for the radiative re- 
action makes fi capture in muonic 3 Hc a troublesome 



bacl :ground. 



some disagreement between these two sources, but that 
has now been resolved by a more careful analysis of the 
corrections necessary in pion electroproduction (Bernard 



In summary, both rZ 10~ J D 2 admixtures in H2 liquid 
and ~ 10~ 2 D 2 admixtures in H 2 gas have been used 
in the study of fid capture. Unfortunately the hyper- 
fine depopulation of fid atoms in H 2 /D 2 mixtures is slow 
and therefore the doublet -quartet make-up in H 2 /D 2 ex- 
periments is dependent on target density and deuterium 
concentration. Further at densities where pfid molecules 
are formed, the observed rate of muon capture is a com- 
plicated superposition of the jip, fid and fi 3 He rates and 
thus disentangling the outcome is difficult. 



et at, 1992). Thus all of the ingredients for a theoretical 



calculation of the OMC rate on the proton are well deter- 
mined by general principles which have been verified in 
many other situations, except for the value of g p , which is 
given primarily by the theoretical predictions discussed 
in Section |l| above. 



2. ChPT calculations 



V. MUON CAPTURE IN HYDROGEN 



The amplitude for OMC on the proton has also been 



A. theory of ord i nary muon capture 



1. Standard diagram calculations 

The simplest of the muon capture reactions is OMC on 
the proton, \i + p — > n + v, which has been studied theo- 
reticall y for many years. So me of the early work include d 
that ofjPrimakoff "(1959] ) and [Fujii and Primakoff (1959| ) . 
Opat (1964) evaluated amplitudes for both OMC and 
RMC in an expansion in powers of 1/mjy. Many other 
authors have performed similar evaluations of the OMC 
rate. 

The basic physics is completely determined by the 
weak nucleon current given in Eq. (|l]) , which is the most 
general possible form consistent with the known current- 
current form of the weak interaction. Given this current, 
the OMC amplitude is determined by its product with 
the leptonic weak current. One then uses standard di- 
agrammatic techniques to square the amplitude, put in 
pha se space and thus obtain an expression for the rate 



in t ;rms of the coupling parameters gv,g m ,ga, and g v - 
This expression is the same whether the couplings are 
obtained from purely phenomenological sources or from 
some detailed fundamental model. 



evalua ted in the context of HBChPT by [Fearing et al. 



(1997). Such calculations start with the general ChPT 
Lagrangian, in this case through 0(p 3 ), and evaluate the 
amplitude consisting of tree and one loop diagrams. The 
outcome of such calculations are expressions for couplings 
appearing in the most general amplitude of Eq. ([!]) in 
terms of the LEC's appearing in the Lagrangian. This 
approach thus provides a systematic way of calculating 
couplings and their corrections, as was described for g p in 
Section pTB| above. However the vector part of the am- 
plitude must still satisfy CVC, g a must still reproduce 
neutron beta decay, etc., so in actual fact the couplings 
to be used are exactly the same as those which have al- 
ways been used in the phenomenological approach and 
there is no new information arising from a ChPT cal- 
culation, except perhaps for the (small) correction term 
appearing in g p , Eq. (^). What is accomplished by such 
a calculation however is to evaluate some of the LEC's 
which are needed for other calculations, such as RMC. 

Similar ChP T evaluations of OM C we re subsequently 
carrie d out by Ando et al. (200"o| ) and Bernard et al. 



(1998 ). In the latter case the so called small scale ex- 



pansion was used, which is a way of including the A as 
an explicit degree of freedom in the ChPT formalism, 
rather than absorbing its effects in the LEC's. 
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FIG. 2 Capture rates for OMC on a proton plotted versus 
g P {— 0.88m^) for the muon initially in the singlet or triplet 
state of the up atom or in the ortho or para state of the pfip 
molecule. The vertical line corresponds to the PCAC value of 
g p = 8.13, as defined in Eq. (S), which includes the constant 
term. Details of the model used are described in Sec. V.D 



B. Theory of radiative muon capture 

We now want to consider the radiative muon capture 
process /i + p — > n + ^ + 7. The basic ingredients are 
the same as for OMC and the additional coupling of the 
photon is known. However the presence of the photon 
changes the range of momentum transfers available and 
leads to contributions coming from regions much closer 
to the pion pole than for OMC. More specifically, for 
OMC on the proton the momentum transfer is fixed at 
q 2 = — 0.88to^. For RMC however the momentum trans- 
fer for some of the diagrams can approach These 
diagrams, all of which involve radiation from hadronic 
legs, are not the dominant ones. The muon radiating di- 
agram dominates, at least in the usual gauge, and it in- 
volves similar momentum transfer as for OMC. However 
the other diagrams contribute enough in the measurable 
photon energy region k>60 MeV that the overall sensi- 
tivity of RMC to <7 p is significantly increased as compared 
to OMC. 



1. Standard diagram calculations 



3. Spin effects 



When a muon is stopped in hydrogen it proceeds via 
a rather complicated series of atomic and molecular pro- 



cess is, as was described above, to a low level state in ei- 



ther a fip atom or a pfip molecule from which the muon is 
captured. The details of this cascade process and relative 
probabilities for population of the various states depend 
on the density of the target and the time at which one 
starts detecting the captures. For now it suffices to note 
that the capture from any initial state will be a linear 
combination of captures from the singlet and triplet [ip 
states. Thus it is important to calculate separately the 
rate 5 from these two spin states. 



The standard approach to RMC on the proton has 
been a Feynman diagram approach, Fig. |^(a)-(e), which 
includes the diagrams involving radiation from the muon, 
the proton, the neutron via its magnetic moment, and 
from the exchanged pion which generates the induced 
pseudoscalar term. The fifth diagram makes t he re- 
sult g auge invariant using a minimal substitution. Opat 



(1964 ) performed one of the earliest such calculations, 
using however a l/mpt expansion of t he amplitude. T he 
completely relativistic calculation of Fearing (1980| ) is 



an example of a more modern calculation using this ap- 
proach. 4 

There are some possible enhancements to this basic ap- 
proach. For example the intermediate nucleon, between 
weak and electromagnetic vertices, can in principl e be a 
A. Such effects, Fig . 5(f)- (g), were considered by |Bcder| 
and Fearing (1987, 1989). They increase the photon 



Figure ^| shows the individual singlet and triplet cap- 
tures ra tes for OMC, u sing a standard diagram calcu- 
198C), plotted versus the value of g p . 



lation 3 (Fearin 



Clearly the capture rate from singlet state is much larger 
than that from the triplet state. However it is also much 
less sensitive to the value of g p . Unfortunately, although 
one can increase the sensitivity to g p to some extent by 
cho osing conditions that enhance the triplet contribu- 



spectrum by amounts rangi ng from 2-3% at 60 MeV to 
7-8% at the upper endpoint. Truhlik and Khanna (2001 ) 
found a similar sized effect. 

Some additional terms, higher order in an expansion of 
the RMC amplitude in powers of the photon momentum 
k or the momentum transfer q, were obtained originally 
by [Adler and Dothan (1966| ). These arise via the re- 



quirements of gauge invariance and PCAC for the full 
amplitude. Similar terms w ere obtained by Christillin 



tion , the singlet rate is so much larger that it dominates and Servadio (1977 ) and by Klieb (1985 ). They how- 



in essentially all circumstances. 



ever seem to be fairly small. 



3 The calculation has been updated to include form factors and 
modern values of the coupling^ particularly g a . The general 
table of numerical values, Table [fl, gives the values of parameters 
and constants used. 



4 Hwang and PrimakofT (1978) also evaluated RMC in hydrogen 
using what they called a linearity hypothesis. This was s hown t.n 



be inc orrect however by Wullschlcgcr and Schcck (1979 ) , Fearing 
(198(f ), and |Gmitro and Ovdmmikova (1981 ). 



10 



TABLE I Numerical values of the parameters and derived quantities used in the text and in our evaluations of rates for 
comparison with experiment. 



Symbol 



Description 



Value 



Reference 



<?irJViv(m^) 
9a (0) 

rl 

PCAC 
9p 



Kortho I A para 
■L^-pfip I *^-pf_ip 

Ao 
27° 

9m(0) 



'op 

ortho 



pion decay constant 

pion nucleon coupling 

Fermi constant for f3 decay 

axial coupling from /3-decay 

rms radius squared for g a 

PCAC value, ffp (-0.88m 2 ) 

PCAC value, including constant term 

Pfip molecular formation rate 

ratio of ortho to para molecular formation 

ortho to para transition rate 

ortho-molecular overlap factor 

para-molecular overlap factor 

weak magnetism coupling, n p — re n 

rms radius squared for g m 

rms radius squared for g v 



92.4 ± 0.3 MeV 
13.05 ±0.08 
1.13548 x 10~ 5 GeV" 
1.2670 ± 0.0035 
0.42 ± 0.04 fm 2 
6.77 g a {0) = 8.58 
6.42 g a (0) = 8.13 
2.5 x 10 6 s- 1 
240 : 1 
4.1 ± 1.4 x 10 4 s _1 
1.009 ±0.001 
1.143 ±0.001 
3.70589 
0.80 fm 2 
0.59 fm 2 



Particle Data Group (200C ) 



de Swart et al. (1997 



Particle Data Group (2000) 



Particle Data Group (200C ) 



Ahrens et al. (1986 
Eq. (§) 
Eq. (!) 



Faifman and Men'shikov (1999|) 



Bardin et al. (1981a) 



Bakalov et al. (1982) 



average, Wright et al. (1996 



Bakalov et al. (1982) 



Particle Data Group (2000) 



Mergell et al. (1996) 



Mergell et al. (1996) 



This basic diagrammatic approach clearly contains 
most of the important physics. However there are some 
things it does not contain. In particular in the simplest 
approach gauge invariance is imposed only via a minimal 
substitution p — ► p — eA on the explicit momentum de- 
pendence of the operators in the weak vertex, Eq. (Q). 
Thus one picks up only two terms, one from the explicit 
q in the g m term and one from the q in the g p term. In 
principle there may be many other terms, gauge invariant 
by themselves, for example coming from situations where 
the photon couples to some internal loop structure of the 
ttNN vertex. 

Also for RMC it is difficult to put form factors in at 
the various vertices in a general way, though lowest or- 
der form factor effects can be included via a low energy 
expansion as done by Adler and Dothan (1966). This 
is because the momentum transfer at the various ver- 
tices in different diagrams is different and so putting in 
form factors evaluated at these different momenta would 
destroy the delicate cancellation needed for gauge invari- 
ance. This is not a problem unique to RMC. It appears in 
any electromagnetic process described by more than one 



introduce only a small correction in RMC as well. This 
turns out to be the case in the simple approach described 
below. 



2. ChPT calculations 

Just as for OMC, it is possible to calculate the am- 
plitude for RMC using HBChPT. Unlike OMC however 
such a calculation introduces some new physics. The 
ChPT Lagrangian is constructed to satisfy gauge invari- 
ance and also CVC and PCAC. Thus the calculated am- 
plitude will satisfy all of these general principles. This 
means for example that such a calculation may contain 
explicitly gauge terms beyond the simple minimal sub- 
stitution used in the standard diagrammatic approach. 
Figure ^ shows some examples of such terms present in a 
ChPT calculation which are not present in the diagram- 
matic calculation. A second advantage of such a calcu- 
lation is that it automatically incorporates form factors, 
at least to the order of the calculation. One can cal- 
culate weak and electromagnetic form fa ctors explicitly 



diag ram if the momentum transfers arc different. Various (1998 ) within ChPT. The same sub-diagrams responsible 



as wa s done by Fearing et al. (1997) or Bernard et al 



prescriptions have been proposed to address this prob- 
lem, but all are pretty much ad hoc. In OMC intro- 
ducing form factors reduces the rate, but only by a few 
percent. This is because the relevant momentum transfer 
is small compared to the scale relevant for the form fac- 
tor. In RMC all momentum transfers are comparable or 
smaller than that for OMC, so hopefully the form factors 



for these form factors will appear in the RMC calculation, 
so the form factors will be present in such a calculation, 
and in a gauge invariant way. 

The ChPT approach has one disadvantage however, 
not shared by the diagrammatic approach. It is an ex- 
pansion in a small parameter. Thus the calculations are 
done at best to 0(p 3 ) which is sufficient to include tree 
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v n 




v n 



n p / t y \ n p / y 




A similar calculation, but just involving tree level di- 
agrams, and working only to 0(p 2 ) was performed by 
Meissner et al. (1998). They initially found a photon 



spectrum harder by 10% or so than that of the diagram- 
matic approach in the region of photon energi es greater 



l?>r 



than 60 MeV, but this was later attributed ( Myhrcr 



(d) 




v u 



n p 



V u. 



Y v 




(c) 




1999) to unwarranted approximations made in the phase 
sp ace evaluation. A sim ilar 0(p 2 ) calculation was done 
by Bernard et al. (2001b| ) who however used the so called 
small scale expansion which allows one to put the A in 
explicitly. Such an approach is not fundamentally dif- 
ferent from the usual HBChPT approach. It just allows 
one to extract explicitly the contribution of the A to the 
various LEC's. In the usual approach the A degrees of 
freedom are integrated out and their effects absorbed in 
the LEC's. These authors found results similar to those 
of earlier diagrammatic calculations and specifically that 
the A contribution was only of order of a few percent, 



consistent with the findings of Beder and Fearing 
1989] ). 



1987 



(f) 



FIG. 3 Diagrams contributing to the standard diagrammatic 
approach to RMC: (a)-(e) are the standard diagrams, (f) and 
(g) are those involving A contributions 



3. Model calculations 

Recently there has been a somewhat different model 
calculation of RMC (Fruhlik and Khanna , 2001) based 



on a -KpojCLi Lagrangian which also includes A's (Sine 



jkal et al., 199E). This is basically a diagrammatic ap- 




n p_ 




^1 

FIG. 4 Several examples of diagrams which appear in a 
HBChPT calculation but which are not included in the usual 
minimal substitution gauge contribution of the diagrammatic 
approach. 



proach, but based on a Lagrangian which has a number 
of additional pieces, and which thus results in a explic- 
itly gauge invariant expression, which also satisfies CVC 
and PCAC and which generates some of the higher order 
terms analogous to those derived by Adler and Dothan 



(1966| ) . For the best values of the parameters the A con- 



tributions range from about 3-7% which can be increased 
to roughly 7-11% at the extreme edge of the allowed pa- 
rameter space. Thus their results are very simi lar to 
those obtained by Beder and Fearing (1987 , 1989| ), and 
so this model would seem to give a result only a few per- 
cent different from the ChPT result or from the standard 
diagrammatic approach including the A. 



level and the lowest contribution to one loop diagrams. In 
principle a diagrammatic calculation keeps all orders of 
whatever diagrams are included. In practice however di- 
agrammatic calculations are limited to tree level graphs, 
so do not include the loop contributions of ChPT and 
probably at most have only a few higher order terms 
coming from relativistic corrections which would not be 
included in the ChPT approach. 

The most complete of the HBChPT calculati ons of 
RMC has been done by |Ando and Min (1998| ), who 
worked to third order, i.e. 0{p A ) 1 or in their terminology 
NN LO (next to next to leading order). They found that 



the loop contributions contributed less than 5% correc- 
tion to the tree level amplitude and that the result was 
in reasonable agreement with the standard diagrammatic 
approach. 



4. Spin effects 

For RMC, just as for OMC, there are fairly dramatic 
differences in the capture rate from singlet and triplet 
initial states, and likewise in sensitivities to g p . Fig- 
ure ^| shows the integrated photon spectrum above 60 
MeV plotted versus g p for RMC on the proton. Now 
it is the triplet state which dominates and the singlet 
which is most sensitive to g p , just the opposite situa- 
tion from OMC. Furthermore, if one looks at the RMC 



photon spectrum before integration (Beder and Fearing 



1989D one sees that variations in g p about the PCAC 



value tend to affect the singlet rate more at the upper end 
of the spectrum than at the lower end, whereas for the 
triplet state the change is more uniform across the spec- 
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FIG. 5 The photon spectrum for RMC on a proton, inte- 
grat ed from 60 MeV to the upper end point, shown as a func- 



tion of g p (— 0.88m M ). Shown separately are the rates for the 
muon initially in the singlet or triplet state of the fip atom or 
in the ortho or para state of the ppp molecule. The vertical 
line corresponds to the PCAC value of g p = 8.13, as defined 
in Eq. M, which includes the constant term. Details of the 



model used are described in Sec. V.D 



trum. This further enhances the sensitivity to g p , since it 
is the upper end of the spectrum which is experimentally 
accessible. 

Finally one should note that the effects which change 
the relative proportion of singlet and triplet states work 
in different directions for OMC and RMC. This is be- 
cause triplet capture is most important for RMC whereas 
singlet capture is most important for OMC. Thus for ex- 
ample, some effect which increases the proportion of the 
triplet state relative to the singlet will decrease the OMC 
rate and increase the RMC rate, for a given value of g p . 



C. Experimental methods for hydrogen and deuterium 

Herein we discuss the experimental methods for muon 
capture on hydrogen and deuterium. They include mea- 
surements of ordinary capture via neut ron det ection, Sec. 
V.C.3 , Michel electron detection, Sec. V.C.2 , and radia- 
tive capture via 7-ray detection, Sec. V.C.4. 



1. Bubble chamber studies 



The first observations of OMC on hydrogen were re - 
port ed by |Hildebrand (1962|), |Bcrtolini et al (1962; ) 



and Hildcbrand and Docde (1962). In these early bub- 



ble chamber experiments the muons were identified by 
range and curvature and neutrons were identified via 
their knock-on protons. An attractive feature of bub- 
ble chamber studies is the determination of the neutron 
energy from the measurement of the proton kinematics. 
This is helpful in distinguishing the neutrons from OMC 



in H2 from combinatorial backgrounds of uncorrelated 
muons and knock-on protons. Recall the neutrons from 
OMC in H2 are mono-energetic E n = 5.2 MeV. However 
the bubble chamber studies were limited by statistics. 



2. Lifetime method 

The 'lifetime method' involves comparing the negative 
muon lifetime and positive muon lifetime when stopped 
in hydrogen or deuterium. For positive muons the life- 
time is the inverse of the \i decay rate, whereas for neg- 
ative muons the lifetime is the inverse of the sum of the 
fi decay rate and the fi capture rate. The lifetime dif- 
ference thus determines the capture rate. This approach 
was used at Saclay to determine the fi capture rate in 
liquid H2 (Bardin et al, 1981a) and liquid D2 (Bardin ei 



al, 1986| ). 



The Saclay experiments used a pulsed beam with typ- 
ically a 3000 Hz repetition rate and typically a 3.0 lis 
pulse width. The beam entered via a lead collimator and 
a copper degrader, stopping in a 24 cm diameter Cu- 
walled target. The target was filled with LH2 or LD2 
with high elemental and isotopic purity. The target was 
viewed by six Michel electron telescopes each compris- 
ing a three-ply plastic scintillator sandwich. A quartz 
oscillator was used to determine the time between the 
beam pulse start signal and the decay electron stop sig- 
nal. Note that the stop signals were recorded in time 
windows of about 1 lis < t < 16 lis following the end of 
the beam pulse. 

The major challenge in the lifetime experiment is the 
required accuracy. Since the fi decay rate is about 1000 
times the \i capture rate, the fi + -fi~ lifetime difference 
is only about 0.1%. Consequently a determination of the 
capture rate A to about 5% requires a measurement of 
the lifetime r to about 5 x 10~ 5 . Therefore both high 
statistics and low systematics are important. 

One concern is muon stops in non-hydrogen materials. 
Detection of decay electrons from extraneous materials 
will alter the time spectrum and distort the measured 
lifetime. Note that the degrader, collimator, and target 
cell were all constructed from high-Z materials, so muon 
stops were rapidly absorbed. Also special care was taken 
to avoid the occurrence of muon stops in plastic scintilla- 
tor. Lastly, H2 or D2 of ultra-high elemental and isotopic 
purity was used. For fip experiments a small quan tity of 



D2 in H2 will result in /i-transfer (see Sec. [V.C ). For 
fid experiments, where muons are qui ckly re cycled from 
dfid molecules to fid atoms (see Sec. tV.B ), the effects 
of muon transfer to high-Z contaminants are especially 
troublesome. 

Another concern is a time distortion due to rate ef- 
fects, i.e. distortion arisin g from finite resolvin g time 
and electron pulse pile-up. Bardin et al (1981a ) found 



the main effect of pulse pile-up is an additional time com- 
ponent with an effective lifetime r/2. They corrected for 
this effect by measuring lifetimes at different rates and 
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extrapolating to zero. The correction to r was roughly 
~ 0.1 ns. Additionally time dependent and time inde- 
pendent physical backgrounds were accounted for. 

In fid experiments a correction is also necessary for fi 
capture on 3 He nuclei. The /i 3 He atoms are produced 



via dfj,d formation and dd fusion (see Sec. 1V.B). The 
3 He background, which was monitored by the detection 
of the 2.5 MeV neutrons from the dd fusion, was a 10% 
correction in the Saclay /id experiment. 

Finally we note that although the Saclay group mea- 
sured both the /i + and \i~ lifetimes they employed the 
then available world average values for r + in order to ex- 
tract the capture rates in LH 2 and LD 2 . The earlier hy- 
dro gen exper iment used r + = 2197.148±0.066 ns ( Bardin 
et a|., 1981a), though this was later updated by Martina 



198 



(19£ 2]7l984D . The later deuterium measurement (jBardin 
used t. 



et a 

discussed in detail in Sec. |V.D 



= 2197.03 ± 0.04 ns. This point is 



contamination the problem of transfer is minimized. Ad- 
ditionally a neutron background generated by the combi- 
nation of Michel bremsstrahlung and (7, n) reactions is 
observed. These photo-neutrons have the 2.2 fis lifetime 
of the n stops. Therefore studies with /i + stops, yielding 
photo neutrons without capture neutrons, are necessary 
to subtract this background. 

Also the accurate determination of the neutron detec- 
tion efficiency is a difficult problem. Calibration via the 
8.9 MeV neutrons from the -n~p — > 771 reaction is helpful, 
but careful simulations of neutron interactions in coun- 
ters, target, etc., are necessary. 

Note that the deuterium experiment using the neutron 
method is especially challenging. First the deuterium 
neutrons form a continuous distribution peaking at 1.5 
MeV. Second a large background is produced by dd fusion 
following d[id formation. 



3. Neutron method 

The 'neutron method' involves directly detecting the 
recoil neutrons from muon capture in H2 or D 2 . For 
hydiURun capLuie Ihe iieuliuiis aie iiiuiiu-enei Relic with 

E n - 



ugen cap I me Ihe 
-5.2 MeV whiL 



iieuLiuiis aie iiiuiiu-eiieigeLie wil 
for deuterium capture the neutrons 



are peaked at E n ~ 1. 5 MeV. Hydrogen d ata i s available 
on l iquid ta rgets from Blescr et al. (1962) and Rothberg 



et a 
et a 



(1973) 



19631) and on gas targets from |Alberigi-Quaranta| 
1969| ) and Bystritskii et al. (1974). Deuterium 
data is available in pure D 2 from ICargnclli et al. (19891) 



and H 2 /D 2 mixtures from Wang et al. (1965) and Bcrtin 
et a 



The 'neutron method' involves (i) counting the incom- 
ing muons and outgoing neutrons and (ii) determining 
their corresponding detection efficiencies. In a typical 
set-up the /x-beam is directed into the target vessel via 
a scintillator telescope and the neutrons are detected in 
a liquid scintillator counter array. Pulse-shape analy- 
sis enables the separation of neutrons from 7-rays and 
veto counters enable the separation of neutrons from elec- 
trons. Both electrons and gammas from /j,-decay are in- 
tense backgrounds. Typically the detection of neutrons 
is initiated about 1 fis after fi arrival. 

The low yield of capture neutrons from OMC in H 2 or 
D 2 means neutron backgrounds from fi capture in sur- 
rounding materials, fi transfer to target impurities, and 
other sources, are troublesome. By using high-Z mate- 



rials for Ihe collimator, vessel, etc., the neutron" back- 
grounds fiuiii /1 stops hi exlianeuus maleiials aie shuil- 
lived. 5 By using high-purity gas/liquid with small Z > 1 



5 Further, the expe riment s of Albcrigi-Quaranta et al. (1969 ), 
^iertin et al. (1973), and Bystritskii et al. (1974) used a counter 



arrangement in the target vessel to define the /i-stops in hydro- 
gen. 



4. Radiative capture 



RMC on H 2 has a yield per muon stop of ~ 10~ 8 and 
a continuum gamma-ray spectrum with E y < 99.2 MeV. 
The first measurement of RMC on H 2 was recently ac- 
complished at TRIUMF fljonkmans et al, 1996] ; |Wright 



et al., 199? 



The experiment detected photons from fi stops in liq- 
uid H 2 . An ultra-pure muon beam {ir/fi = 10~ 3 ) was 
directed into a liquid H 2 target via a scintillator tele- 
scope. The target flask and vacuum jacket were con- 
structed from Au and Ag in order to ensure wall stops 
in high-Z materials. The flask contained pure hydrogen 
with a D 2 contamination of about 1 ppm and a Z > 1 
contamination of less than 1 ppb. The 7-ray detector was 
a high acceptance, medium resolution, pair spectrometer 
and comprised a lead cylinder for 7-ray conversion, cylin- 
drical multiwire and drift chambers for e + e~ tracking, 
and axial B-field for momentum analysis. The photon 
detection efficiency was efl ~ 10 -2 and photon energy 
resolution was AE/E ~ 10%. 

The major difficulty in RMC on H 2 is the tiny yield. 
Consequently photon backgrounds from pion stops in liq- 
uid hydrogen, muon capture in nearby materials, and fi 
decay are dangerous. 

Pion stops in liquid H 2 undergo both charge exchange, 
ir~p — > Tr°n, and radiative capture, ir~p — > 771. The re- 
sulting photons comprise a Doppler spectrum of 55-83 
MeV and monoenergetic peak at 129 MeV, th us endan- 
gering th e region of interest for RMC on H 2 . Wright et 



al. (199? ) suppressed this background via an ultra-high 



purity muon beam, prompt photon timing cut, and the 
difference in the n/fi ranges. Additionally the authors de- 
termined the residual background from the 55-83 MeV 
photons via the residual signal from the 129 MeV 7-ray 
peak. 

The bremsstrahlung of electrons from /x-decay yields 
a continuum background with E^ < It prevents 

the measurement of the RMC spectrum below 53 MeV 
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and threatens the measurement of the 
above 53 MeV, because of the finite 
pair spectrometer 



Wright et al (199 



ibrt 



spectrometer to minimize the contr 
energy tail in the response function 
remaining background from Michel 
measured by stopping a /x + beam in 
the photon background from /i decay w 
signal from /i capture. 

In addition the backgrounds from /i si 
materials and fj, transfer to target 
mized by (i) using high-Z materials in 1 
and (ii) using high isotopic purity H2 
rial. Also ancillary measurements wen 
termine the photon backgrounds from 
and cosmic-ray interactions. 

The experiment recorded 397 ± 20 



RMC spectrum 
resolution of the 
3) designed their 
ion of the high- 
Additionally the 
brcmsstrahlung was 
quid H2, yielding 
thout the photon 



ops in extraneous 
impurities were mini- 
he target vicinity 
1 the target mate- 
employed to de- 
abcelerator sources 



ergies E 1 > 60 MeV and times t > 305 
tracting the backgrounds from ^t-decay 
(48 ± 7), Au/Ag radiative ijl capture 
sources, a total of 279 ± 26 photons f 
uid H2 was obtained. The resulting 
ratio for RMC on liquid H2, i. e. th 
ton spectrum for k > 60 MeV divided 
muon decay and capture rates, was R, 
(2.1 ±0.21) x 10~ 8 . Note this value 



ahotons with en- 
ns. After sub- 
bremsstrahlung 
±11), and other 
om RMC in liq- 
mrtial branching 
integrated pho- 
yy the sum of the 
(k > 60MeV) = 
esponds to the 



COT' 



par ticular occupancy of the uv> spin-s 



UMF experiment, i.e. with t > 365 ns 
sities. 



D. Comparison of experiment and theory 



cf 



In Table |n] we summarize the results 
measurements of muon capture on hy 
experiment we have shown the density 
relative to liquid hydrogen, the approximate 
from muon stop until counting starts 
eraged proportion of singlet, ortho an 
evant for the experiment, and the rat 
these conditions. Note that the sing! 
tios depend on the muon chemistry, 
ing parameters determining that chem stry. 
lay time At. In Table || we used a 



mation rate A T 



2.5 x 10 6 



state formation to para-state formatic 
mar and Mcn'shikov , 1999]; Ponomar 



197$ IZel'dovich and Gcrshtein , 19591) 



transition rate A 



op 



and gamma factors of 2j ortho = 1.009 
(Bakalov et al, 1982). Furthermore 



tic 



correspond to different experimental 
are not directly comparable. 

First we stress that the various expe 
tive to different combinations of the fj, 
ular states, i.e. the triplet, singlet, orth 
The situation is simplest in the H2 



Alberigi-Quaranta et al (1969) and 



ates in the TRI- 



at liquid H2 den- 



of the various 
Irogen. For each 
the target n/no, 
time delay 
At, the time av- 
I para states rel- 
obtaincd under 
let/ortho/para ra- 
ind the underly- 
and the de- 
molecular for- 
a, ratio of o rtho- 
n of 240:1 (jfalf- 



pip 



and Faifman 



an ortho-to-para 



( Bar din et al, 1981a ) 



a: id 



2l para = 1143 

rates tabulated 
clonditions and so 



iments are sensi- 
£ tomic and molec- 
d and para states, 
experiments of 



Bystritskii et al 



(1974) where muon capture is almost exclusively from 
the F = atomic state and p/ip molecule formation is a 
few p ?rcent correction. However for experiments in liquid 
H2, while capture from the ortho state is the largest, cap- 
ture from other states is significant. The precise blend of 
stated is determined by the delay At between the muon 
1 time and the counting start time. The greater 
3lay the smaller the contribution from the F = 
c state and the larger the contribution from the 
nolecular state, as the muon has additional time to 
bhe pfip molecular state and convert from ortho to 
state, in accord with the processes described in Sec. 
Thus for example the bubble chamber experiments, 
wherij At = 0, have therefore the largest contribution of 
atom capture and smallest contribution of para 
capture. In contrast the Saclay lifetime experi- 
where At ~ 2.5 /xs, had the smallest contribution 
glet atom capture and, though still dominated by 
capture, the largest contribution of capture from 
state. 6 

;he last column of this table is given the value of 
responding to the experimental rate. To get these 
numb|ers the singlet and triplet capture rates were cal- 
d and combined as appropriate for the experimcn- 
cojnditions of the individual experiment. The model 



arrivj 

the d 

atom 

para 

form 

para 

fv 



single t 
molequle 
ment 
of sin 
ortho 
the 

In 
9p co 



p ira ; 



culat' 
tal 

used vas the standard diagrammatic approach of Fearing 



coupl 
factois 
be of 
of the 
0.80 
and 
turn t 
typic, 
theor 



peariji; 

was 



te 



9p 



Here 
just 
of g P 
and 

For 
l + q 
were 
gaug 
founc 
form 
the r 



(1980 ) updated to the extent that modern values of the 
ngs, particularly g a = 1.267, were used and form 
were included. These form factors were taken to 
the form f(q 2 ) = 1 



1 tn 



q 2 < r 2 > /6 with the values 
rms radii square d < r 2 > taken as 0.59 fm 2 and 
2 for g v and g m (Mergell et al, 1996 ) respectively 



42 fm 2 for g a (Ahrens et al, 198S). At the momen- 
ransfer appropriate for OMC these form factors are 
ally 0.98 — 0.96 so that including them reduces the 
tical OMC rate by 2-4%. The constant term ap- 
g in g p , as in Eq. (||) or ([)]), was included and g p 
ameterized as 



par 



> , _ D 2m tl m N rnv,'m N g a (Q)r 1 A 



9%(Q ) = R 



(ml - q 2 ) 



9a (0) - 



(25) 



R is not intended to have physical significance, but 
a conventional way of parametrizing the variation 
from the PCAC value. At the PCAC point, R = 1 



(-0.88m 2 ) = 8.58 - 0.45 = 8.13. 



RMC form factors were included also using the 
< r 2 > /6 form and the necessary gauge terms 
generated via a minimal substitution. This gives a 
inv ariant result which includes most of the terms 
by |Adler and Dothan (19661 ). For RMC these 
actors make essentially no difference, e.g. < 1% in 
te corresponding to the TRIUMF experiment, pre- 



At " 2.5 /is is an approximate average value for the Saclay ex- 
peri nent. In analyzing the data the detailed time structure had 
to bp explicitly treated. 
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TABLE II Summary of world data for OMC and RMC on hydrogen. The columns correspond to the target density in units of 
liquid hydrogen density n a , the time delay between /i stop and start of counting, the effective, time averaged, singlet/ortho/para 
ratio corresponding to the particular experimental conditions, the capture rate corresponding to these ratios, and the value of g p 
implied using the calculation described in the text. For RMC the value given in the rate column corresponds to _R 7 (fc > QQMeV), 
the partial branching ratio obtained by integrating the photon spectrum above k — 60MeV and dividing by the sum of the fi 
decay and capture rates. Note that the ratios of molecular states in the S:0:P column depend on the parameters, e. g. Ann. 
use d, and ar e given only for those experiments that measure the neutron or gamma yield. For the Saclay experiment (Bardin 



1981a ) the relationship between the measured rate and the ortho capture rate depends on details of the experiment and 



et a, 

is given in the original paper. 



Ref. 



n/n D 


At(ns) 


S:0:P 


rate (s 1 ) 


g p (-0.88ml) 


1.0 


0.0 


0.15:0.77:0.07 


420 ± 120 


19.5 ± 11.6 


1.0 


0.0 


0.15:0.77:0.07 


428 ± 85 


18.7 ±8.2 


1.0 


0.0 


0.15:0.77:0.07 


450 ± 50 


16.4 ±4.9 


1.0 


1.0 


0.01:0.88:0.11 


515 ±85 


6.3 ±8.7 


1.0 


1.2 


0.01:0.88:0.12 


464 ± 42 


11.4 ±4.2 


0.014 


0.9 


1.00:0.00:0.00 


651 ±57 


11.0 ±3.8 


0.072 


1.4 


1.00:0.00:0.00 


686 ± 88 


8.7 ±5.7 


1.0 


2.5 




460 ± 20 


7.9 ±3.0 








435 ± 17 


10.6 ±2.7 


1.0 


0.365 


0.06:0.85:0.09 


(2.10 ±0.21) x 10~ 8 


12.4 ± 0.9 ±0. 



OMC 



Hildebrand 



Hildebrand and Doede (1962) 



Bleser et al. 



(1962j) 



Bertolini et a l. (196 2[) 
1962|) 



Rothberg et al. (1963) 



Albcrigi-Quaranta et al. (1969) 
Bystritskii et al. (1974 ) 
Bardin et al. (1981a ((original r+) 
(new r+) 

RMC 



Wright et al. (199^ ) (original theory) 
(new theory) 



12.2 ±0.9 ±0.4 



sumedly because both spacelike and timelike momentum 
transfers contribute and, with the linear approximation 
for t he form factors, tend to cancel. T he A was included 
as in |Beder and Fearing (1987j |l989[) . 

The OMC results obtained from this calculation are in 
very good agreement with other modern OMC calcula- 
tion s such as those of Ando et al. (200fj| ) and Bernard et 
al. ( 2001b] ) . It is interesting to observe however that the 
values of g p quoted in Table |TJ are 0.3-0.8 higher tha n 



those in a similar table given in |Bardin et al. (1981b ) 



This can be traced to two main effects, namely the in 
crease in g a to its modern value, 1.254— >1. 267, and the 
use of more modern form factors w hich fall somewhat les s 
rapidly with q 2 than those used in Bardin et al. (1981b). 



Both of these effects lead to a larger theoretical rate for 
a given value of g p and thus, as can be seen from Fig. |^, 
to a larger g p to fit a given experimental rate. 

A further comment is required concerning the value 
of g v obta ined from the Saclay experiment ( Bardin et 



al., 1981a| ,[b[). In analyzing their data the authors ex- 



tracted the capture rate from the difference between their 
measured value for r_ and the world average value for 
t + . They used the world average for comparison because 
their measured value for r + , although consistent with the 
world average, had a larger uncertainty. Since the pub- 
lication of Bardin et al. (1981a) the world average of t + 



has changed from 2197.15 ± 07 ns to 2197.03 ± 0.04 ns 
flParticle Data Group , 200C| ). [n determining g p for 
Table || we decided it best to update the /i + lifetime 
in extracting the (x~ capture rate, 7 which now becomes 
435 ± 17s" 1 instead of 460 ± 20s" 1 . Therefore the value 
of 5p (-0.88m2) = 10.6 ± 2.7 in Table || is one standard 
deviation larger than g p (— 0.88m^) = 7.9 ± 3.0 which 
is obtained from the original published rate with our 
theoretical calculation, and even larger than the value 
ff p (— 0.88m^,) = 7.1 ± 3.0 given in the original paper 



(Bardin et al., 1981b) 



The values of g p (—0.88m 2 l ) obtained from the five more 
recent, electronic, OMC experiments are all in general 
agreement within their errors and result in a world aver- 
age value 10.5 ± 1.8. The sin gle determination w ith the 
smallest uncertainty is that of Bardin et al. (1981~b| ) which 
gives 10.6±2.7. Both of these values are somewhat larger 
than the PCAC prediction g p (-0.88m£) = 8.13, which 
includes the constant term, though are consistent with 
it at the 1 to 1.3 standard deviation level. Note that 



7 This was also done, in a conference proceedings, h v one ol 
the m embers of the original experimental group. See Martino 



(1984 ). 
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if o ne does not update the fi + lifetime in the Bardin ei 
al. ( 1981b| ) analysis the Saclay result is g p {— 0.88m^) 



7.9 ± 3.0 and world average is g p (-0.8Sml) = 9.4 ± 1.9, 
both of which are also consistent with PCAC, but both 
still about 0.7 larger than the corresponding numbers, 
7.1 ± 3.0 and 8.7 ± 1.9 quoted in the original paper, be- 
cause of the updates to the parameters of the theory. 

However the result # p (-0.88ra£) = 12.4 ± 0.9 ± 0.4, 
or 12.2 ± 1.1 using the updated th eory, from the TRI - 
UMF radiative capture experiment jWright et at, 199i ) 
is 50% larger than, and clearly inconsistent, with PCAC. 
Note that the RMC result is quite consistent with the 
world average from OMC or the Saclay result for g p if 
one uses the updated fi + lifetime. However, with the 
older //+ lifetime that Bardin et al. (1981b| ) have used, 
the Saclay result and RMC result do not overlap within 
their uncertainties. 

We can thus summarize the situation in hydrogen as 
follows. The RMC result is several standard deviations 
larger than the PCAC prediction and clearly inconsistent 
with it. The OMC results have always been considered 
to be in agreement with PCAC, based on the results of 
the Saclay experiment. We have seen though that, as a 
result of updates to the parameters in the theory and to 
subsequent measurements of the /i + lifetime, the value 
of g p obtained from the OMC result has increased. The 
increase is only about one standard deviation, so the re- 
sult is still marginally consistent with PCAC. The central 
value however is now also high, and actually somewhat 
closer to the RMC result than to the PCAC value. 



E. Attempts to resolve the discrepancy 

1. General comments 

In the previous section we saw that the results for g p 
obtained from RMC were significantly higher than the 
PCAC prediction and in definite disagreement. The new 
interpretation of the OMC results suggest that they are 
too high also, being in good agreement with the RMC re- 
sult, but perhaps only marginally consistent with PCAC. 
Since the older interpretation of the OMC data seemed to 
be in good agreement with PCAC, there has been essen- 
tially no consideration of possible difficulties with OMC. 
However the discrepancy between the RMC results and 
PCAC has generated a lot of discussion and a number of 
attempts to find additional effects to explain it, some of 
which will be now discussed. 

A first question one asks is it logical to have a different 
g p for RMC and OMC? At a fundamental level the value 
of g p should be the same. However given the present state 
of analysis, an apparent difference could arise simply be- 
cause RMC is much more complicated than OMC and 
involves a lot of additional diagrams. Thus if something 
is left out of the RMC analysis, a fit to the data using the 
standard approach may require a different value of g p to 
compensate for the piece left out. In this view a difference 



in the values of g p extracted from OMC and RMC using 
current theory may reflect something wrong or missing in 
the theory rather than a failure of the PCAC relation of 
Eq. (||). Note however that while a difference between the 
values of g p extracted from OMC and RMC might be ra- 
tionalized this way, this would not explain any differences 
between the OMC and PCAC results since for OMC g p 
is a parameter in the most general weak amplitude. It 
can thus be determined from data largely independently 
of the details of the weak capture part, as opposed to the 
molecular and atomic part, of the theory. 

A second question to ask is what does it take to bring 
the RMC result closer to PCAC, and what does that do 
to OMC? As discussed earlier the rates for both RMC 
and OMC depend strongly on the proportion of singlet 
versus triplet components in the initial state. To get the 
RMC result to agree better with PCAC we need to in- 
crease the predicted rate for a given g p , which can be 
done by increasing the amount of triplet capture and re- 
ducing the singlet capture. Since OMC is dominated by 
the singlet rate, this reduces the OMC rate, which also 
moves the extracted g p toward the PCAC value. However 
the sensitivity of OMC and RMC to increasing triplet is 
different, and thus it becomes difficult to completely fix 
one without destroying the agreement of the other. 8 

In any case, in the subsequent paragraphs we will dis- 
cuss several aspects of the muon chemistry which change 
the average singlet /ortho/para ratio for each experiment 
and thus in principle affect the comparison of the results 
with PCAC. We will also discuss two other suggestions, 
dealing specifically with the RMC calculation which have 
been put forward as potentially resolving the discrepancy 
implied by the RMC results. 



2. Value of the ortho-para transition rate A op 

As has been noted, the experimental combination of 
singlet and triplet states is important, and this is de- 
termined by the various transition rates governing the 
chemical processes the muon undergoes between stopping 
and capture. The least certain of these rates is A op , the 
ortho-para transition rate. There is a single theoretica l 



value A op = (7.1 ± 1.2) x lO 4 * -1 ( JBakalov et al., 1982; ) 



8 To get simultaneous agreement of both the TRIUMF RMC result 
and Saclay OMC result with PCAC a larger increase in triplet 
occupancy is needed for RMC than for OMC. Since the TRIUMF 
RMC experiment detected the majority of its photons with muon 
capture times less than a few fis, whereas the Saclay OMC ex- 
periment detected the majority of its electrons with muon decay 
times more than a few /j.s, in principle such circumstances are 
possible. For example a long-lived, i.e. ~1 fis lifetime , triplet 
atom component could increase triplet occupancy for the TRI- 
UMF RMC experiment much more than for the Saclay OMC ex- 
periment. However, such a long-lived triplet atom is completely 
at odds with our current knowledge of muon chemistry in liquid 
hydrogen. 
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FIG. 6 The region in the g p - A ov plane which is al lowed by 
the TRIUMF RMC exper iment ([Wright et al. 1998| ) and the 
Sacl ay OMC experiment ( Bardin et al, 1981a ). The dashed 



line fndicates the PCAC value 
Eq, 



9^' 



8.13, as defined in 



which includes the constant term. Details of the 

Note that in deter- 



model used are described in Sec. V.D 



mining the error bands the uncertainty included in the exper- 
imental results due to the uncertainty in A op has not been 
included, as the results are being plotted against A op . 



which however is nearly a factor of two l arger than the 
experimental value ( 4.1 ± 1.4) x 10 4 s _1 ( Bardin , 1982 ; 
Bardin et al, 1981b|), 9 though a nother experiment is in 
progress ( Armstrong et al., 1995 ) . A larger value of A op 



increases the relative amount of para state, and thus in 
creases the amount of the triplet component contained in 
the initial state. Figure |^ shows the region in the q r , - A 



pla ne allowe d by the TRIUMF RMC experime nt (|Wright 



L981a) 



1998 ) and the latest OMC experiment ( Bardin et 



OTJ 



et a 
al 

For these experiments an increase in the value of A op , 
which leads to a greater triplet component in the initial 
state, increases the predicted RMC rate and decreases 
the OMC rate for a given value of g p . In the vicinity 
of the PCAC value of g p , the RMC rate increases and 
the OMC rate decreases with increasing g p . This means 
that for fixed value of the rate the extracted value of g p 
decreases with increasing A op for both OMC and RMC, 
as can be seen from the figure. However the Saclay OMC 
experiment is much more sensitive to A op than the TRI- 
UMF RMC experiment. 

The updates and improvements we have made in the 
theory, which raise the OMC band in this figure, make 
it a little easier to accommodate a larger value of A op 
and still have consistency between OMC and RMC. In 
fact now a modest increase in A op from the experimental 
value which has normally been used improves the agree- 



This group also obtained a value 7.7 ± 2.7 X 10 s 1 via a m ea- 



ment of both OMC and RMC with PCAC while keeping 
them consistent. However increasing A op enough to get 
agreement of OMC with PCAC or increasing it further 
to the theoretical value still leaves RMC several standard 
deviations above the PCAC prediction. Increasing it suf- 
ficiently to produce the PCAC value of g p from the RMC 
data would lead to a catastrophic disagreement between 
the Saclay OMC experiment and the PCAC prediction. 
Thus while one can envision a value of A op which slightly 
improves the agreement with the PCAC prediction for 
both OMC and RMC, there still appears to be no value 
which will simultaneously result in good agreement with 
PCAC for both cases. 



3. Admixtures of a J = 3/2 ortho-molecular state 



It was pointed out in the very early work of Wein 
berg (196C|) that the ortho state could also include a 



component of total spin angular momentum S = 3/2 



Subsequent theo ret ical calculations ( Bakalov et al, 1982 
Halpcrn , 1964a| ,[b|; |Wcsscl and Phillipson , 1964 ) seemed 



to indicate that this component was zero or very small. 
Nevertheless if there were such a component it would ef- 
fectively increase the amount of triplet state in the initial 
sta te. This has been s u ggeste d as a possible explanation 
by Ando et al. (200C , 2001a ). Their original numbers 



surement of the electron time spectrum. See [Martino (1982| ). State, A 



were not quite right, due to a misinterpretation of the 
appropriate initial state in the OMC experiment, but the 
idea is worth examining. Figure [7] shows the region in the 
V~9p plane allowed by the TRIUMF RMC experiment 
and the most recent OMC experiment. Here r\ is the 
fraction of the spin 3/2 state present in the ortho state. 
The OMC calcula tion here differs slightly from that of 
Ando et al. (2000| ) by virtue of the fact that form fac- 
tors have now been included. The theoretical prediction, 
rj = corresponding to no additional spin 3/2 compo- 
nent, gives the most consistent result. The updates to 
the OMC theory, which again raise the band from the 
Saclay experiment, make it possible to slightly improve 
both RMC and OMC results with a small non-zero value 
of T]. However analogous to the previous case, a value of 
?y which makes OMC agree with PCAC leaves RMC too 
high and a value which makes RMC agree leads to dra- 
matic disagreement of OMC. Thus this effect does not 
appear to resolve the discrepancy. 



4. Direct singlet - para transitions 

As shown in Fig. |l] it is in principle possible for the fip 
singlet atom to make a direct transition to the para p[ip 
molecule. This also would increase the proportion of the 
triplet component in the initial state. This transition 
is governed by the rate which is predicted to be 

~ 0.75 x 10 4 s _1 for liquid hydrogen, which is two orders 
of magnitude smaller than the transition rate to the ortho 
1 



i ortho 



x 10 6 s \ ([Faifman and Men'shikov 
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FIG. 7 The region in the g v - r\ plane which is al lowed by 
the TRIUMF RMC exper iment ( | Wright et al, 1998[ ) and the 

1981a|). 



Saclay OMC experiment (Bardin et al. 

— „fUAU 



The dashed 

line indicates the PCAC value, g^ 0/1 ° = 8.13, as defined in 
Eq. (0), which includes the constant term. Details of the 
model used are described in Sec. V.D. A 



has been taken 
as 4.1 x 10 4 s _1 and the uncertainty in its value has not been 
included in the plotted error bands. 



1999), and so should be a negligible effect. 



5. Large A resonance effects 

In a recent calculation of RMC in a Lagr angian based 
mo del, describ ed in Section V.B.3 above, Truhlik and 
Kha |nna (20"0l| ) argue that the discrepancy is at least par- 
tially resolved by a contribution coming primarily from 
the A. The first observation to be made however is that 
using their best parameters the A contribution they ob- 
tain is essenti ally the same size as the contr ibution ob- 
tained earlier ( Bcdcr and Fearing , 1987 , 1989] ) which was 
already included in the analysis of the TRIUMF RMC ex- 
periment. Other values of the parameters, but still within 
the region allowed by other experiments, can raise this 
A contribution by 3-4% in the most important part of 
the spectrum, where the discrepancy is more than 40%. 
It would appear, as they in fact note, that their primary 
results really basically agree with previous ones and so 
do not provide an explanation of the discrepancy. They 
go on to observe however that one can change one of 
the parameters corresponding to off-shell properties of 
the A and thus fit the experimental spectrum. That is 
certainly possible, as one can usually arbitrarily change 
one parameter to get a fit to another. However it re- 
quires the A contribution to be an order of magnitude 
larger than what they consider most reasonable. Fur- 
thermore one would have to reconcile such an explana- 
tion with the general result that changing off-shell prop- 
erties can not change phys i cally m easurable quantities. 
See e. g. Fearing (1998b|, |2000[); Fearing and Schcrcr 
(20fj^); (gcherer and Fearing (2001 ) and references cited 



therein. Although it is clear that there are uncertainties 
in any attempt to include a A in a calculation such as 
this, their best estimate of the A contribution is quite 
consistent with the (small) contribution found in the two 
previous independent calculations which included the A 
flBeder and Fearing , 1987|, |l989|; [Bernard et al, 2001b|). 



An enhancement obtained by varying off-shell parame- 
ters which is large enough to explain the full discrepancy 
does not seem too likely. 10 



6. Other possibilities 



Anot her suggestion was made by |Bernard et al 



(2001b), namely that the discrepancy between PCAC 
and the RMC results could be resolved by a series of 
small effects in RMC which happen to add up. Such ef- 
fects might include for example those discussed above or 
perhaps higher order loop effects which have not been 
included, radiative corrections, or variations in some of 
the parameters. Clearly one can not test this suggestion 
without calculating all such effects, which has not been 
done. We note however that in so far as the molecular 
effects we have discussed are concerned, they all serve to 
increase the amount of triplet capture, at least for current 
liquid H2 experiments. However in each of these cases, 
or in some combination of them, as they are somewhat 
equivalent, it appears that changes sufficient to make 
RMC agree, lead to drastic disagreement for OMC. 

Another particular effect they proposed was an isospin 
breaking effect. This proposal was based on the obser- 
vation that the RMC rate changed quite a bit when the 
charged pion mass in the pion propagator was replaced by 
the neutral pion mass, thus supposedly showing a sensi- 
tivity to the kinds of presumedly electromagnetic isospin 
breaking effects responsible for the pion mass splitting. 
This sensitivity clearly originates in the fact that some 
of the momentum transfers in RMC can get quite close 
to the pion pole. Thus changing the position of the pole, 
even by a few MeV, can change the results significantly. 
One should note however that the pion exchanged is in 
fact a charged one. Thus in an ideal theory which prop- 
erly incorporated all electromagnetic effects, after all the 
renormalizations, and after all the terms which lead to 
mass splittings were included, the pion mass appearing 
in the propagator would in fact be the physical charged 
pion mass, not the neutral one. Thus this particular ar- 
gument for sensitivity to isospin effects in spurious. In 
such a complete theory one would probably expect terms 
in the numerator of the amplitude proportional to the 
pion mass splitting, and likewise the neutron-proton mass 



10 There was also a suggestion by Cheon and Cheoun (1998, 1999) 
of an additional term in the Lagrangian which could contribute, 
but it was shown that such a term was already include d in stan 



1998) 



dard calculations, and so was not an additional effect (Fearing 
1998a; |5mejkal and Truhlik 
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splitting. But the scale in the denominator would likely 
be something like the pion mass, so these terms should be 
small. Nevertheless, there could be accidental enhance- 
ments, so this is a calculation which should be done. 

One might also want to consider radiative corrections 
to both OMC and RMC. Typically these might be ex- 
pected to be 0(ct/7r), and this was wha t was found by th e 
one such calculation we are aware of (Goldman , 1972). 
This is too small an effect to be relevant now, but might 
be important in the interpretation of future, high preci- 
sion, OMC experiments. 



F. Summary and outlook for hydrogen 



ties. Such a measurement might distinguish between 
two possible scenarios. In one, both OMC and RMC 
may give values of g p which are too high and inconsis- 
tent with PCAC, which would suggest that the prob- 
lem is something common to OMC and RMC, perhaps 
some difficulty with our understanding of the molecu- 
lar or atomic effects. In the other OMC might clearly 
agree with PCAC and disagree with RMC, which would 
suggest that the problem is something wrong with our 
understanding of RMC. 

It would also help to confirm some aspects of the muon 
chemistry, particularly the value of A op which is impor- 
tant in determining the initial muonic state and for which 



the ex isting experimental (Bardin , 1982 



Bardin et 



1981b) and theoretical (Bakalov et at, 1982) values differ 



It is thus possible to summarize the situation in hy- 
drogen as follows. From the theory side, all calculations 
of OMC in hydrogen are in essential agreement, as they 
must be as all couplings and form factors except g p are 
fixed by generally well accepted principles such as CVC or 
by other experiments. The several calculations of RMC, 
which in principle could differ, e.g. via the extra loop 
contributions included in the ChPT calculation relative 
to the standard diagrammatic calculation, are in fact in 
general agreement at the few percent level. 

On the experimental side, all modern OMC experi- 
ments now generally lead to values of g p which are con- 
sistent with each other within the errors. The percep- 
tion that OMC agrees with PCAC has changed somewhat 
however with updates to the theoretical calculations and 
to the n + lifetime. These updates have increased the 
world average from OMC by only about one standard 
deviation, but result in a value which is now higher than 
the PCAC value by about 1.3 standard deviations and 
which is in fact in better agreement with the RMC re- 
sult than with PCAC. The one existing RMC experiment 
gives a value of g p which is 50% larger than the PCAC 
pre diction and with uncertainties small enough that it is 



by al most a factor of two. Su ch an experiment is under 
way (Armstrong et ai, 1995) and should have results 
soon. 

Another alternative is to try to measure some other 
quantity which is especially sensitive to g p . One such 
quantity is the triplet capture rate for OMC in H2. Here 
the central problem is that the triplet rate is about thirty 
times smaller than the singlet rate. Further during the 
time evolution of the /ip system following atomic capture, 
the singlet state is always present, so the singlet capture 
dominates the OMC rate at all times. However some sug- 
gestions have surfaced for isolating the triplet capture 
from the mu ch larger single t capture. One possibility 
suggested by Deutsch (1983| ) is to measure the time de- 
pendence of the neutron polarization. Neutrons following 
capture from the singlet state have polarization -1, while 
those from the triplet state have a polarization depending 
on the weak dynamics and g p . Measuring as a function of 
time allows one to determine a relative neutron polariza- 
tion, rather than an absolute one, which is more difficult. 
Such a measurement of the neutron polarization could, 
in principle, help in identifying the triplet capture and 
isolating the g v contribution. Another suggestion (Bai 



in c jear disagreement with that prediction. Though there ley et al, 1983| ) notes that the triplet /xp atoms, unlike 



have been a number of suggestions and attempts to ex- 
plain the discrepancy indicated by the RMC result, none 
so far have been successful, and the difference remains a 
puzzle. 

So what needs to be done next to resolve this situa- 
tion? Specifically a high precision measurement of the 
OMC rate in gaseous hydrogen would be very helpful. 
By using gas rather than liquid the muon chemistry un- 
certainties are greatly reduced. Indeed a new measure- 
ment of the singlet capture rate for OMC in H2 gas at 
P ~ 10 bar using the lifetime method is currently un- 
der way at PSI (Kammel et al., 2000). The goal is a 
precision of about 1% in the singlet rate and about 6% 
in the coupling g p . Note that the experiment will use 
a novel hydrogen time projection chamber that enables 
direct monitoring of muon stops, muon transfer and /i- 
atom diffusion. 

Also useful would be a new measurement of OMC in 
liquid hydrogen with significantly improved uncertain- 



the singlet /ip atoms, will precess in a magnetic field. 
Consequently the time spectra of decay electrons or cap- 
ture neutrons will be modulated by this precession when 
muons are stopped in H2 gas at low pressures, where 
triplet atoms are relatively long-lived, thus enabling iso- 
lation of triplet capture. While interesting and worth 
reconsidering, such experiments are of course extremely 
difficult. 

For RMC, it is clear that a new experiment in liquid 
hydrogen, where triplet capture is largest, would also be 
both interesting and useful to check the TRIUMF result. 
Further a new RMC experiment in gaseous hydrogen, 
where singlet capture is largest, would be very interest- 
ing. However an RMC experiment in gaseous hydrogen 
is extraordinarily difficult since the muon stopping rate 
is much less and the singlet RMC rate is much smaller 
than the triplet RMC rate. 

Since for RMC the triplet rate is by far the largest, it 
is possible, by utilizing both a gas target, where capture 
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is primarily from the singlet state, and a liquid target, 
where there can be significant captures from the triplet 
state, to obtain as the dominant contribution to the RMC 
rate either singlet capture or triplet capture. This is dif- 
ferent from the situation for OMC since there the singlet 
rate is the largest and so in both gas and liquid experi- 
ments the singlet capture will dominate the total rate. 

In principle it is possible to measure spin degrees of 
freedom in RMC as well and a suggestion of th is type 
has also been made recently by Ando et al. (2002] ) . Such 
experiments would be extremely difficult because of the 
low rates and the necessity of either starting with a po- 
larized muonic atom or measuring the polarization of the 
outgoing photon, but they are more sensitive to g p than 
the rate and are sensitive to a different combination of 
amplitu des and so give independent information on the 
prouess ( Feaiing , 1970 y. 



perhaps less so than originally in view of the advances 
in our knowledge of the nucleon-nucleon force from other 
sources. However the sensitivity to a nn means that it 
must be well known to have any hope of using this pro- 
cess to extract g p , at least if one focuses on the region of 
low relative neutron momentum where the rate is largest. 

Furthermore the total OMC rate, which is dominated 
by the doublet rate, is only moderately sensitive to g p . 
A 50% cha nge in g v produc es only about a 10% change 
in the rate (Doi et al., 1991). 



The total rate is also affected by MEC's which have 
been calculated using specific d iagrams involving pion 
or rho exchange (Doi et al., 1990) or using a more elabo- 
ratc hard pion Lagr angian involving 7r's, p's, o/s a nd a\ ' s 



199C 



et al., 



Ivanov and Truhlik , 1979b 



Tatara 



( Adam et al. 

1990| ). Generally these MEC's are about a 10% 



correction to the total rate. 



VI. MUON CAPTURE IN DEUTERIUM 



We next want to consider OMC on deuterium, i.e. the 
reaction fi + d—>-n + n + is, which involves many of the 
same ingredients as OMC on the proton, but also some 
additional complications. 



A. Theory of muon capture in deuterium 

Since there are now two nucleons involved in the cap- 
ture process, several new ingredients appear which were 
not relevant for capture on the proton. In particular 
the process now depends on the two nucleon interac- 
tion 



The results of modern calculations now generall y agree. 
For example the doublet capture rate obtained by Tatara 
et al. (1990| ) is 398-400 s" 1 of which 31-33 s^ 1 comes 



from MEC's and where the ra nge corresponds t o different 
nucleon-nucl eon potentials. Dpi et al. (1990| ) obtained 



402 s" 1 and |Adam et al. (1990| ) found 416±7 s 



It it interesting to note that since this is a three body 
final state, just like RMC in the proton case, it is pos- 
sible to get timelike momentum transfers approaching 
and thus get near the pion pole, where in princi- 
ple the sensitivity to g p is enhanced. To our knowledge 
sensitivit y to q v in that region has not been investigated, 
although Mintz (1983) considered so me e ffects of time- 



which is reflected in the properties of the initial Reform factors and poi et al. (199Cj ) and poulard et al 



stat e dcutcron, such as the wave function , percentage D 

state, etc., and also in the final state interaction between 
the two outgoing neutrons. Furthermore now meson ex- 
change corrections (MEC's) can contribute. 

There have been many calculations of the OMC rate 
in deuterium. 11 Some of the more recent ones include 
Tatara et al. (1990|), |Adam et al (1990|), |Doi et al. (1990|, 



(1982) showed that meson exchange currents are an im- 



1991 



Morita and Morita (1992| ), [Hwang and Lin (19991 ) 



and Ando et al. (2001b| ). Originally the main emphasis 
was on using this process to extract the neutron-neutron 



portant effect there. Unfortunately this kinematic region 
occurs when the two neutrons go out back to back and 
the v has relatively low energy, which is a region which 
is strongly suppressed relative to the region where the 
neutron-neutron final state scattering is important, thus 
making experiments difficult. 

There is a very strong hyperfine effect in the capture 
rate on deuterium. The quartet rate is only 10-15 s ^ 1 



scattering length < 



In the kinematics m which the two 



199C 



( Doi et al.. 
19T4|). It has been suggested 



Ho-Kim et al., 1976: Sotona and Truhlik 



)oi et 



trurii peaks al a value several orders of magmlude larger 

than that for kinematics with large relative momentum. 
Such a measurement of this low energy neutron spec- 
trum, while extremely difficult, is still interesting, but 



1991; Morita 



nem roiis have low rclaLive momeiiLum, Lhe ractnnr spec- and ^ Qrita ; m $ fcgj^ ~ ± a 7; 1993 |) that t ' h ^t£ 



1 1 Some representative e arl ier ones include th ose of [Wang ( 1965k } 



us 



(1974 ). Dogotar et al. (1979 LE 

ano v ana 

|M gl iven~gS7^,|Ho-Kirn kTal. Q97(|4Svar.- | al n^ pi 
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Sotona and Trnhl 



[Vlintz (1973| |l983| ) 



md .Bajzcr ( 19S0|y 



varc 

oulard et al. (1985 ), and 



of quartet to doublet rate is much more sensitive to g p 
than the doublet rate and thus that a measurement of 
the quartet rate, or directly of the ratio, would give a 
better value of g p . 

In principle RMC in deuterium should also be an in- 
teresting process, giving information on g p . However one 
might expect that MEC's would be even more important 
here, since at the least the photon would have to cou- 
ple to all of the intermediate states which generate the 
MEC's in OMC. To our knowledge there have not been 
calculations of RMC in deuterium. 
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B. Experiments on muon capture in deuterium 

A summary of available data for /id total capture rates 



is g ven in Tabic [II . It lists experiments with mixed 



H2/D2 targets by Wang et al. (1965) and Bcrtin et al. 
(19^) and pure D 2 t argets by Bardin et al. (1986) and 
Cargnelli et al. (1989| ). A H 2 /D 2 target is helpful 111 rc- 
ducing the neutron background from dfid fusion. The 
Bardin et al. (1986| ) re sult w as obtained using the 'life- 
time method' (see Sec. V.C.2 ) and the other resu lts wer e 
obtained using the 'neutron method' (see Sec. V.C.3| ). 
Note the neutron measurements give partial rates for 
muon capture with neutron energies > E t hres and re- 
quire some input from theory to determine the total rate 
of [id capture 



The measurement by Wang et al. (1965 ) used a liquid 



H 2 target with a 0.32% D 2 admixture. Under these con- 
ditions the [i capture is dominantly from pfid molecules, 
since muon transfer yields /id atoms in roughly 20 ns and 
proton c aptur e yields pfid molecules in roughly 200 ns. 
See S ec. |IV.C for details . Consequently the neutron sig- 
nal in Wang et al. (1965 ) involves a superposition of cap- 



ture on protons, deuterons and 3 He, the latter being pro- 
duced via pd fusion in pfid molecules. Note that in ana- 
lyzing their results the authors assumed a statistical mix- 
ture of the /id spin-states in the pfid molecules, as would 
be expected if hyperfine depopulation by spin-flip colli- 
sions is much slower than pfid formation. Also the differ- 
ent time dependence of /i 3 He neutrons and /id neutrons 
was employed in separating their contributions. 12 Their 
final result for doublet capture was A 1 / 2 — 365 ± 96s _1 . 

The measurement by Bcrtin et al. (1973| ) used a 7.6 bar 
H 2 gas target with a 5% D 2 admixture. Under these 
conditions the /id atom formation is rapid but p/id 
molecule formation is negligible, thus avoiding the com- 
plication of contributions from protons, deuterons, and 
3 He. Strangely they found that under the assumption of 
a statistical mix of the /id spin-states their extracted rate 
for doublet capture was 1100-1500 s _1 , i.e. three times 
the theoretical value. Therefore, in order to understand 
their results, they postulated a much faster hyperfine de- 
population rate than assumed by Wang et al. (1965). 
Thus assuming a pure doublet mix of /id spin states the 
authors obtained a doublet rate of 445 ±60 s _1 . How- 
ever the subsequent studies of hyperfine depopulation in 
H 2 /D 2 systems, e.g. |Breunlich (198l"|), have failed to 



sup port this claim and therefore the [Bcrtin et al. (1972 ) 



resu lt io extremely puzzling. 



T he more recent exper iments of Cargnelli et al. (198£ ) 
and Bardin et al. (1986) were conducted in pure deu- 
terium. In these circumstances the /id + d collision rate 
is sufficient to fully depopulate the F + state. Note any 



12 The neutrons from fid capture build-up quickly with the 
nanosecond timescale of the fi transfer process. The neu- 
trons from /i 3 He capture build-up slowly with the microsecond 
timescale of the pd fusion process. 



dfid formation is followed by prompt fusion and muon 
recycling into /id atoms. The resulting r ates for doublet 
captu re were A ly / 2 = 409 ± 40 s _1 from Cargnelli et al. 



(1989 ), using a gas target and the neutron met hod, and 
A 1/2 = 470 ± 29 s" 1 from |Bardin et al. (1986| ), using a 



liquid target and the li fetim e method. Using the c alcula- 
tion of Doi et al. (1991), the |Cargnclli et al. (1989| ) result 
implies g p = 8 ± 4 and the Bardin et al. (1986[ ) result 



implies g p — 2 ±3, which are in disagreement at the level 
of about 1-1.5 0. Although one experiment used a gas 
target and the other experiment used a liquid target, in 
both cases it is believed that capture is from the doublet 
atomic state, so that muon chemistry in pure deuterium 
is unlikely to account for the different results. However 
both experiments were challenging and faced different ex- 
perimental difficulties. Therefore at present we conclude 
the two results simply indicate a rather large uncertainty 
in the value of the coupling g p extracted from the /td 
system. 



There has also been one measurement of the high en- 
ergy neutron spectrum following muon capture in deu- 
terium (Lee et al., 1987) using neutron time of flight 
methods. In this kinematic situation the neutrons come 
out back to back, so the direction is well defined, but 
there is no start signal for time of flight measurements. 
This group used a novel technique to overcome this, by 
using the detection of one of the neutrons in a counter 
placed close to the target as a start signal for the time 
of flight measurement of the other neutron in a counter 
placed further away from the target. The results sug- 
gest the e nhancement due to me son exchange corrections 
found by Goulard et al. (1982 ) and are in qualitative 
agreement with the calculations of Doi et al. (199Cj ) at 
least as far as the shape of the spectrum is concerned, 
but are not precise enough to give any real information 
on g p . In view of the potentially increased sensitivity to 
g p however, such a measurement of high energy neutrons 
is worth re-examining. 



A new higher precision measurement of the /td dou- 
blet capture rate would be useful in resolvin g the possible 



discrepancy between the current results of Bardin et al. 



(1986) and Cargnelli et al. (1989). However in deuterium, 



unlike in hydrogen, the uncertainties in calculating the 
contributions of exchange currents may ultimately limit 
the achievable accuracy in extracting g p . We also re- 
mark that in pure deuterium and hydrogen-deuterium 
mixtures the hyperfine depopulation rate is much slower 
than in pure hydrogen. Consequently measuring either 
the quartet rate or hyperfine dependence in //d capture, 
which are strongly dependent on g p , may be experimen- 
tally easier in /id. It would be worthwhile to reconsider 
this source of information on g p , as has been suggested 



by Morita and Morita (1992) and Doi et al. (1991) 
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TABLE III Summary of world data for the fid doublet capture rate. The columns correspond to the target density in units 
of liquid hydrogen density n a , target temperature, the time delay between /i stop and start of counting, the neutron energy 
threshold, the deuterium concentration, and the doublet capture rate. 



Ref. 


n/n 


T(K) 


At(fis) 


E n (MeV) 


Cd 


rate (s x ) 



Wang et al. (1961 



Bertin et al. (1973) 



Bardin et al. (1986) 



Cargnelli et al. (1989| ) 



1.0 
0.013 

1.0 
0.04 



18 
293 
18 
40 



0.8 
0.5 
2.5 
1.9 



1.4 
1.5 

1.5, 2.5 



0.0032 
0.05 
1.00 
1.00 



365 ± 96 
445 ± 60° 
470 ± 29 
409 ± 40 



" Note that the Bertin et al. (1973) doublet capture rate was obtained assuming a much faster hyperfme depopulation rate 
than is consistent with our current understanding of muon chemistry in H2/D2 mixtures. See text for details. 



VII. MUON CAPTURE IN HE 

We now want to consider both OMC and RMC on 3 He. 
Since it is possible to break up the final state there are 
several processes possible, 



fj, 



+ 3 He 



'H + v 



6 He -> A H 



[i + He -^p + n + n + v 



(26) 



(27) 



(28) 



with analogous processes involving a photon in the case 
of RMC. 

Of these, most work has been done on the triton final 
state, Eq. (p6]). Like the deuteron, since there is more 
tha n one nuclcon , the complications of exchange cur- 



- 3 H system is a spin 1/2 isodoublet just like the proton- 
neutron system. Thus one can carry over directly the cal- 
culations done for capture on the proton. The couplings 
and form factors now apply to the three-body system as a 
whole, rather than the nucleon. However they can be ob- 
tained in many cases directly from other processes such as 
electron scattering on 3 He or tritium /3-decay or from the- 
ory, i.e., by applying PC AC to the nuclear pseudoscalar 
form factor. This approach has the advantage that it 
already includes effects of MEC's, which are buried in 
the form factors, and that it can be done relativistically. 
However it is difficult to relate the nuclear form factors to 
the fundamental couplings and form factors of the weak 
current, e. g. g p . Some sort of microscopic model is nec- 
essary to make this connection. 

One of the first such EP M calculations was that of 
Kim and Primakoff (1965). The idea was la t er ap- 



rcntk must be addressed Additionally three-body forces 

are now possible. However there has been a tremendous 
amount of work done on three-body wave functions in 
recent years and these wave functions are now extremely 
welli known , Thus the wave function complications which 



plied to 3 Hc by Bedcr (1976[), Fearing (1980|) , |Klicb 



and Rood (1981), and Klieb (1982). A more modern 



and very careful calculation with particular attention to 
the form factors and the various uncertainties was car- 



plagiie interpretation nf rapture in heavier nuclei really 

do not enter here. Furthermore the rate is significantly 
higher than in capture on protons and deuterons and the 
atomic and molecular processes which make the interpre- 
tation of capture in hydrogen or deuterium difficult are 
not present here. Thus 3 He is an ideal compromise, and 
we will see that one of the best measurements of g p to 
date comes from OMC in 3 He. 



A. Theory of ordinary muon capture in 3 He 



ried out by Congleton and Fearing (1993), and their 
result s were c onfirmed by povacrts and Lucio-Martincz 
(2000| ) and by |Ho et al. (2002| ). The final result is a pre 



diction from t he EPM for the statistical capt ure rate of 
1497±21s~ 1 ( Congleton and Fearing , 1993 ), using the 
standard PCAC value for g p . 

The IA on the other hand treats the capture process 
as capture on an essentially free proton. This one nu- 
cleon interaction is then expanded in powers of 1 /m and 
used as an effective interaction in a matrix element be- 
tween nuclear wave functions. The information specific 
to the nucleus comes in only via the nuclear wave func- 
tions. Early application s of this approach app l ied to 3 He 
include Peterson (1968 ), Phillips et al. (1975 ), Klieb and 



ory of muon capture in~^He, the elementary particle 
mo del (EPM) and the impulse approximation (IA), sup- 



T [iere have been two major approaches to the the- Rood (1981 ), and Klieb (1982). Again a very carefu l 



calculation was done by Congleton and Fearing (1995 ) 



plcmcnted by explicit calculations of meson exchange 
corrections. The EPM exploits the fact that the 3 He 



using the variational wave functions of Kameyama et al._ 
(1989| ) with the result fo r the rate 1304 s -1 . This was 



essentially confirmed by Ho et al. (2002) using for the 
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wave functions momentum space Fadd< ev 
variety of potentials. They found howejve 
spread among the results using the 
tions. The fact that the I A result is ab 
the EPM is characteristic of the IA and 



solutions for a 
r about a ±3% 
wave func- 
15% lower than 
indicate s the ef- 
fect of the MEC's. |Conglcton and Truhlik (1996| ) refined 



o it 



sxplic 



these IA calculations by including e: 
based on a specific Lagrangian involvm. 
plings. The end result, which also included 
of three-body forces, was 1502 ± 32 s 
ccllcnt agreement with the EP M rest 
recently Marcucci et al. (2001) repeated 
of |Conglcton and Truhlik (1996j ) usir 



It. 



some 



of 1 



arising from the Argonnc 1*14 or vis 
three-body forces. They also made 
the weak current and reduced the size 
by fixing the A contribution to the excl 
rcctions via a fit to tritium beta deca; 
1494 ± 9 s" 1 . 

Thus to summarize, the capture rate 
is well determined by both EPM and 
tions with an accuracy of 1-3%. We not 
of accuracy radiative corrections might 
should be evaluated. 

Finally it is worth mentioning that th 
calc ulations of various spin correlatior 
gleton , 1994; Congleton and Fearing 



itly the MEC's, 
g 7r, p, 01, A cou- 
some effects 
This is in ex- 
Finally, very 
the calculation 
wave functions 
pot entials with some 
refinements to 
the uncertainty 
ange current cor- 
". They obtained 



for OMC on 3 He 
IjA+MEC calcula- 
that at this level 
be important and 

re have also been 
s in OMC (|Con^ 



and Truhlik , 1996|; Govaerts and Lucie 



Marcucci et al., 2001). Generally thes 



more sensitive to g p and less sensitive 
the rate is. There have also been 



break up channels ( Phillips et al. 
199* |Wang , 1965a^ 



B. Experiments on OMC in 3 He 

The 3 H channel in OMC on 3 He 
trino with momentum 100 MeV/c and 
energy 1.9 MeV. The triton yield per 
0.3%. Contrary to muon capture on cc 
recoil triton has a relatively high energ 



small charge, and is thcrctorc directly d etectable 



below we discuss the current 



VILB7T 

— > 3 H capture rate experiments and 
3 He — > 3 H recoil asymmetry experimer ts 



m 



1. ' 5 ie capture rate experiments 



So far several measurements of the 
have been performed. They comprise 
Aucrbach et al. ( 1967 ), Clay et al. (1G 65), Falomkon e\ 



al. (1963), Zaimidoroga et al. (1963), aid recent work of 



Ackerbauer et al. (1998) 



The basic method involves counting \ 
stops and 3 He-recoils when a beam of 
in 3 He. In Aucrbach et al. (1967) and 



1995 ; |Congleton 



Martinez , 2000 



correlations are 
t(b the MEC's than 
calculations of some 



1975: Skibinski et al. 



;lds a muon neu- 
triton recoil with 
<i stop is roughly 
mplex nuclei, the 
y and a relatively 



(19981 
Fa 



111 

liquic 
scinti 
the 



omkon et al. (1965 ) and Clay et al. (1965 ) a 3 He 
scintillator was used. One advantage of a liquid 
lation counter is the higher muon stopping rate in 
liquid helium target. However one advantage of a gas 
tion chamber is the better separation of the \x and 
gnals via its tracking capabilities. 



; recent experiment of Ackerbauer et al. (1998) em- 



iomz£ 
3 He s 
Tli 

ploye I an 3 He ionization chamber with a 15 cm 3 active 
volume and a 20 Torr gas pressure. The anode plane 
segmented to permit the tracking of the incoming 
the isolation of the recoil tritons, and the defini- 
f the fiducial volume. The anode strips were read 
iihto flash ADCs for measurement of the drift time 
the energy loss. Separate triggers were employed 
mt the incoming muons or recoil tritons with full 



cie ncy. 



was 
muon|3 
tion 
out 
and 
to co 
effii 

On 
over! 
3 He 
probl 
latin 
(ii) k 
shap 



i Hiization. Ackerbauer et 



dures 
back, 
to- 

Th 
-» 3 I 
tial p 



jrour 



over 
1410 



In Sec 
-a- 



status ot the J rie~ 
Sec. IVII.B.2I the 



lish t 

as 

safely 



well 



(1998) 



(1996 ) 



3 He 3 H rate 
;hc early work of 



he numbers of \i- 
muons is stopped 



Ackerbauer et al. 



13 For 
tain 
obt; 



a 3 He gas ionizati on ch amber was used w hereas 



experimental difficulty is muon-triton pile-up, i.e. 
pping of the prompt ionization and the delayed 

(1998| ) addressed this 



m by both (i) ignoring pile-up events and extrapo- 
to t = the non pile-up /u- 3 H time spectrum, and 
?eping pile-up events and using the anode pulse- 
(■s for computing ^- 3 H time differences. The con- 
sistency of (i) and (ii) was helpful in demonstrating the 
correct treatment of pulse pile-up. 

Another experimental difficulty is background pro- 
cesses. The 1.9 MeV triton energy-loss peak is super- 
imposed on backgrounds that include 3 H break-up chan- 
nels, thermal neutron c apture, and /z-decay electrons. 



Ackeibauer et al. (1998) employed various fitting proce- 



in order to separate the triton peak and continuum 
round and estimate the uncertainties. Their signal- 
background ratio was roughly 20:1. 

1998| ) experiment yielded a 3 He 



Ackerbauer et 



rate, corresponding to a statistical hypcrfinc ini- 
Dpulation, of 1496.0 ± 4.0 s _1 , with the largest sys- 
tematic uncertainty being the subtraction of the back- 
d. The ±0.2% accuracy is a ten-fold improvement 
he earlier experiments of Falomkon et al. (1965] ) 



t 140 s" 1 , |Clay et al. (1965|) 1467 ± 67 s' 1 , and 



|Aucrbach et al. (19671 ) 1505±40s 1 . In addition [Acker- 



bauci et al. (1998| ) utilized their time spectrum to estab 



lat transfer from the F_ to the F + hyperfine state 
as capture from the 2S metastable state can be 
neglected. 13 



Taken together the measurement of Ackerbauer et al 



and the calculation of Congleton and Truhlik 
give g p (— 0.954m^) = 8.53 ± 1.54 which implies 
3p(-(|l.88m*) = 8.77 ± 1.58 = (1.02 ± 0A8)g^ CAC , or 
(1.08 ± 0.19)3^° if we use the value of g^ CAC which 
includes the constant term. These results are in nice 



the hyperfine transition rate Ackerbauer et al. (199^ ) ob- 
id (0.006 ± 0.008) fis' 1 and for the 2S state lifetime they 
ined < 50 ns. 
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agreement with the PCAC prediction for the coupling 
g p . Here the dominant uncertainty originates from the 
theoretical calculation and not the experimental data. 
Clearly the key question is the exact size of the theoret- 
ical uncertainty. We have quot ed the theoretical error of 
Conglcton and Truhlik (199€ ) which originates mainly 
from th e A contribution to th e exchange currents. Subse- 
quently Marcucci et al. (2001 ) have suggested this uncer- 
tainty can be reduced by fixing some of the parameters 
using results from exchange current effects in tritium beta 
decay. However Ho et al. (2002| ) found significant theo- 



retical uncertainty arising from different choices of the 
nuclear wave functions, though some of this may be due 
to the fact that three-body forces were not included and 
tuned to give the same binding energies for all choices. 
Further efforts on establishing and reducing the theoret- 
ical uncertainty in calculating the rate for 3 He — > 3 H 
capture are definitely worthwhile. 

Lastly we note that there have been some recent mea- 



sur e i nients of paiLial caplme rates and ciiuigv spcclia 
for the breakup channels following muon capture in 3 IIc 

Jummings et al, 1992; Kuhn et al, 1994). Comparisons 



with a recent calculation^ Bkibinski et al., 1999) indicate 
that a good treatment of final state effects is necessary 
to describe the data. However at the present stage these 
channels give no information on g p . 



2. 3 He recoil asymmetry experiments 



As discussed for nuclei in general in Sec. IX. B. 2 below 



the recoil asymmetry in muon capture can be very sen- 
sitive to the coupling g p . Unfortunately production of a 
measurable 3 H asymmetry requires production of a large 
/i _3 He polarization, which has proven to be difficult. 

Early attempts to polarize /i~ 3 He, using a polarized 
beam on a unpolarized target and an unpolarized beam 
on a polarized targ et, were unsuccessful. For example, 
Souder et al. (1975) obtained a /i~ polarization of only 
~6% w hen stopping polarize d muons in unpolarized 3 He 
gas and Newbury et al. (1991 ) obtained a fi~ polarization 
of only ~7% when stopping unpolarized muons in polar- 
ize d 3 He gas. 14 The b reakthrough was made at LAMPF 
by Barton et al. (1993 ) by repolarizing the /i~ 3 He system 



with a laser-pumped Rb vapor after the atomic cascade. 

How does repolarization work? Following the /i _3 He 
atomic cascade a (^ _3 He) + positive ion is produced. 
However the (/x _3 He) + ion is short lived, forming ei- 
ther (/z _3 He) +3 He molecules via collisions with sur- 
rounding He atoms or (/x _3 He) + e~ atoms after collisions 
with electron-donor impurities. The repo larization tech- 
nique developed by Barton et al. (19"93| ) employs both 
spin transfer in the dissociation of the (/i~ 3 He) +3 He 



molecules 



Rb t +He{p-He) + -> Rb + + He + {p-He) + e~ \ 



(29) 



and spin transfer in the collisions of the (/z 3 He) + e 
atoms 

Rb t +{fi-He)+e- |-> Rb J, +{^He)+e- \ (30) 

where the arrows denote the polarization state of trans- 
ferred electrons. Following the dissociation process of 
Eq. ( p9| ) and collision process of Eq. ( pO] ) the polariza- 
tion of the electron i s shared with the m uon via their 
spin-spin interaction. Barton et al. (1993| ) discovered at 
high Rb| densities both processes contribute to muon re- 
polarization on the time scale of the muon lifetime, and 
produced a /x~ polarization in 3 He gas of 26%. 

Recently the first measurement of the recoil asy mmetry 
in the 3 He — > 3 H transition was accomplished by Souder 



et al. (1998) at TRIUMF. The experiment utilized a 3 He 



14 In both cases the polarizatio ns were smaller than expe cted from 
atomic cascade calculations ( Reifenrother et al., 1987). 



ionization chamber to stop the incoming muons, repo- 
larize the /i _3 He system, and track the triton recoils. 
For experimental d etails see Bogorad et al. (1997) and 
3oudcr et al. (1998). The chamber had an instrumented 
volume of 140 cm 3 and a gas pressure of 8 bar. Anode 
segmentation assisted the definition of the fiducial vol- 
ume and pulse-shape read-out assisted the separation of 
the /x/ 3 He pulses. Running the chamber at high volt- 
age, high temperature, and high Rb density was a great 
achievement. 

Souder et al. (1998) employed an interesting approach 
in determining cos 9, the direction between the 3 H recoil 
and /x 3 He polarization. First the magnitude of cos 9 was 
determined from the anode signal width. If the recoil 
travels perpendicular to the anode plane the drift time 
range and anode pulse width are large. If the recoil trav- 
els parallel to the anode plane the drift time range and 
anode pulse width are small. Second the sign of cos0 
was determined from the anode pulse shape. For a re- 
coil traveling towards the anode the early-arriving Bragg 
peak yields a fast rise time pulse. For a recoil traveling 
away from the anode the late-arriving Bragg peak yields 
a slow rise time pulse. 

The /i 3 He polarization was determined by detecting 
the /i-decay electrons and using the well known correla- 
tion between the electron momentum and the muon spin 
in the /i — + evv decay. The electrons were detected us- 
ing a telescope consisting of wire chambers and plastic 
scintillators, which permitted ray tracing into the ioniza- 
tion chamber. The asymmetry in the electron counts, on 
reversing the laser pol arization, gave P^. 

Souder et al. (1998) achieved a muon polarization of 
P M = 30 ±4%. By normalizing the observed triton asym- 
metry with measured muon polarization the authors ob- 
tained a preliminary value for the vector asymmetry for 
3 He -> 3 H capture of A v = 0.63±0.09(siai.)±°;^(si/si.). 
The largest uncertainty was the systematic error in the 
experimental determination of the muon polarization. 
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Taking the preliminary result of Bouder et al. (1998 ) 
and the model calculation of |Conglcton and Truhlik 
(199^) we obtain, at q 2 = -0.954m£, gp / g PCAC = 
0.40±gf§, which implies that g p (-0.88m*) = 3.3±£f. 
This result is consistent with PCAC and the /x 3 He cap- 
ture data, but the experimental uncertainties are quite 
large. 



C. Radiative muon capture in He 

1. Theoretical calculations 

Much less has been done on RMC in 3 He than on 
OMC, though there have been calculations of the exclu- 
sive rate both in IA and EPM. Calculations in the EPM 

who how- 



were made by Hwang and Primakoff (1978 ) 
ever made an incorrect assu mption, and th u s obtained 



2. Experimental results 



In 3 He RMC the 3 H-channel has a strong peaking at 
the recoil energy corresponding to the kinematic limit 1.9 
MeV. In contrast the energy spectra of charged particles 
from break-up channels are a broad continua reaching to 
53 MeV and therefore are straightforwardly distinguished 
from 3 He — » 3 H. 

While the 3 He— > 3 H transition following 3 He RMC and 
the p— >n transition following X H RMC are spin-isospin 
analogs, the 3 He — ► 3 H process is experimentally easier. 
First the 3 He rate is much larger and consequently the 
background difficulties are greatly reduced. Second the 
/i 3 He atom circumvents the various chemical processes 
that confound muonic hydrogen. Of course for 3 He RMC 
the presence of exchange currents and the contribution 
of break-up channels have to be considered. 



inco rrect resu lt s, and also by Fearing (1980 ), Klieb and 



Roold (19811), Klieb (19821) and Ho et al. (20021). For The » channel in 3 He RMC has recently been mea- 



RMC however the EPM seems to be less reliable than for 
OMC. The nuclear form factors are more rapidly varying 
than the nucleon form factors so that t he gauge terms in- 
volving derivatives of the form factors ( Adler and Dothan 
lgj|; phris tillin and Scrvadio , 1977t |Ho et al, 2002| 



Klieb , 1985) are somewhat larger. F urthermore, as em- 
phasized by Klieb and Rood (1984a b|), some important 
terms involving intermediate state excitations are missed 
in the EPM. 



sured by |Wright et al. (2000| ) at TRIUMF. In the experi 
ment a [i~ beam was stopped in liquid 3 He and gamma- 
triton coincidences recorded. The photons were detected 
using a pair spectrometer and the tritons were detected 
via their scintillation light in the liquid helium. Pulse- 
shape read out permitted the discrimination of the in- 
coming muon from the recoil triton. The detection effi- 
ciency was calibrated via 7r stops in liquid 3 He, i.e. with 
the well-known branching ratios for 7r _3 He — > 



T he RMC rate h as also been c alculat ed in IA by Klieb 



Rood (198t ), [Klieb (1982] ), and 



Ho et al. 



(2002|) 



r-3 



He 



r o3 



7 3 H and 



H. 



and 

The earlier calculations used wave functions from the 
Reid soft core potential and made a number of approxi- 
mations in evaluat ing the nuclear matrix elements. The 
later calculation of Ho et al. (2002) improved on a number 
of these approximations, kept some additional terms, and 
used a variety of nuclear wave functions obtained from 
modern momentum space Faddeev calculations. The end 
results for the ph oton spectrum however were about the 
same as those of |Klieb (1982| ), though the photon po- 



larization was somewhat different, apparently because of 
the higher order terms kept. 

Typical results for the RMC statistical capture rate 
for photons with energies greater than 57 MeV are 0.1 7 
s- 1 in the IA and 0.25 s" 1 in the EPM (|Ho et al., 20021) . 
Just as for OMC, the IA result is significantly below that 
obtained in the EPM, and presumedly the difference is 
again at least partially explained by MEC's, which will 



be nhuch more complicated here than they were tor UlVlC, 
and have not yet been calculated in detail. In view how- 
ever of the problems with the EPM applied to 3 He, as 
noted above, it is not possible to simply ascribe the dif- 
ference just to MEC's. Such corrections will have to be 
calculated explicitly to finally obtain a reliable rate. 

In principle the various breakup channels would give 
additional information, but to our knowledge there have 
not been any calculations of such channels for RMC. 



In the 3 He RMC experiment the major photon back- 
grounds originate from muon decay in the target ma- 
terial, muon capture in the nearby materials, and pion 
contamination in the muon beam. However energy cuts, 
i.e. E 1 > 57 MeV, and timing cuts, i.e. t 1 > 440 ns, to- 
gether with the gamma-triton coincidence requirement, 
were effective in discriminating the 3 He signal from back- 
ground processes. Note the potential background from 
random coincidences of photons from RMC and tritons 
from OMC was negligible. 

Only preliminary res ults are presently av ailable for the 
3 He RMC experiment flWright et al, 20001 ). The results 
indicate that the shape of the measured 3 He RMC energy 
spectrum is consistent with the theoretical prediction of 
Klieb and Rood (198l| ), but that the overall magnitude 
is 20-30% smaller than the impulse appro ximation prc- 
diction. To get agreement in magnitude Wright et al 
(2000P stated that one needs g p — 3.4, which is much 



smaller than the PCAC value. 

We stress that the 3 He^ 3 H radiative capture rate is 
potentially a very valuable data-point for the determina- 
tion of g p . We therefore urge the publication of the final 
result from the TRIUMF experiment and the undertak- 
ing of a modern impulse approximation plus exchange 
current calculation for the process. 
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VIII. OTHER FEW BODY PROCESSES 

There have also been attempts to extract g p from other 
processes, notably pion electroproduction. To do this 
Choi et al. (1993) measured the near threshold 7r + elec- 
troproduction at several momentum transfers. The con- 
nection to q v was made via a low energy theorem de- 
scribed by Vainshtein and Zakharov (1972 ) and reviewed 
by Drcchscl and Tiator (1992). The theorem is valid in 
the limit of vanishing pion mass, so it had to be assumed 
that hi gher order correc tions were in fact small. Never- 
theless Choi et al. (1993) obtained results which they in- 
terpreted as g p {t) over a range of four-momentum trans- 
fer squared 0.07 < t < 0.18 GeV 2 . These results agreed 
quite well with the predicted pion pole dominance and 
with the PCAC value near t = 0. 

It is clear that such an approach to 'measuring' g p is 
quite different than that described for OMC and RMC 
above. In fact there has been a quite heated disagree- 
ment over the question of whether or not pion electro- 
pro duction does give new i nformation about cither g v or 
g a ( Bernard et al., 2001c|; Guichon , 2001; Habcrzcttl , 
200% |2001| ; fTruhlik , 200% . 

In OMC g p is defined as a parameter in the most gen- 
eral form of the amplitude, so that a measurement of 
the rate gives directly information about g p . In RMC 
this same amplitude appears directly, albeit with one leg 
off shell and with additional gauge terms, some of which 
probe the interior structure of the amplitude. To the ex- 
tent these are small, one again in RMC has an almost 
direct measurement of g p . 

In other processes however, such as pion electropro- 
duction, one needs to make a number of assumptions in 
order to make the connection to g p . For example, the 
low energy theorem approach is based on chiral symme- 
try, which implies the usual PCAC relation between g a 
and g p , which in turn allows one to eliminate one or the 
other from the expressions for the amplitude. Alterna- 
tively in a Born approximation diagram approach, which 
is often used, a pion pole diagram is included. This is the 
diagram which is assumed to be responsible for g p . Thus 
in such cases is one measuring g p or just the importance 
of a pion exchange diagram? 

It thus appears that there is a degree of circularity in 
trying to use processes such as pion electroproduction to 
'measure' g p . To make the interpretation, one has to as- 
sume something nearly equivalent to the PCAC relation 
between g a and g p , which is what is being tested. In 
our view then processes such as pion electroproduction 
are very interesting, and very worth measuring, but pri- 
marily as consistency checks on our understanding of the 
underlying theory, rather than as independent measure- 
ments of g p . 

Finally one should note that in the context of ChPT 
the calculations of OMC give a relat ion between g v and 



could be considered as a way of 'measuring' g pi but only 
within the context of ChPT. 



IX. EXCLUSIVE OMC ON COMPLEX NUCLEI 

Determining g p from exclusive OMC on complex nuclei 
is not easy. First, the reaction products are a ~100 MeV 
neutrino and a ~0.1 MeV recoil, and therefore the j3- 
rays or 7-rays from the recoil's decay must be detected. 
Second, the effects of g p in OMC are small and subtle, 
and therefore measurements of spin observables are usu- 
ally required. Last, the observables are functions of both 
the weak couplings and the nuclear wave functions, and 
disentangling these contributions is unavoidably model 
dependent. 

In allowed transitions on Ji = targets the Op tran- 
sition 12 C(0+,0) -> 12 B(l+,0) and ls-Od transition 
28 Si(0+,0) -> 28 A1(1+,2201) have attracted the most 
attention. In order to isolate the coupling g p , the 12 C 
experiments have measured recoil polarizations via the 
12 B /3-decay and the 28 Si experiments have measured 
7-ray correlations via the 28 Al 7-decay. In allowed tran- 
sitions on Ji ^ targets the Op nucleus 11 B and ls-Od 
nucleus 23 Na have been studied. These experiments mea- 
sured the hyperfine effect in 3/2 ± — > 1/2^ transitions to 
extract the coupling g p . Last, the capture rate of the 
16 O(0+,0) ^ 16 N(0", 120) transition has been measured 
by several groups; since the + — » 0~ spin sequence offers 
special sensitivity to g p . 

The mat erial on exclusive OMC is organized as follows. 
In section IX. B we describ e the physical observables in 
partial transitions. In Sec. [X.C we discuss the dynam- 
ical conten t of the physical observables in OMC and in 
Sec. IX. D we describe the specific manifestations of the 



coupling g p in OMC. The experimental work on exclu- 
sive transitions is covered in Sec. [X.E. In Sees. IX. F 



and IX. G respectively we review the theoretical frame- 
work and nucl ear mo dels for exclusive OMC on complex 
nuclei. In Sec. IX. H we discuss the results for 
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some of the low ene r gy constants (Ando et at, 2000 
Bernard et at, 1998; Fearing et al., 1997^ Thus any 



process which contains these same low energy constants 



One goal is a unified discussion and comparison of the 
different experiments and calculations. Another goal is a 
critical assessment of the current 'best values' and model 
uncertainties in the extraction of the coupling g p . 

We don't discuss the related topic of total OMC rates 
on complex nuclei. Although a large body of accurate 
data is nowadays available on total OMC rates, because 
the total rate is weakly dependent on g p and the model 
uncertainties in total rates are large, extracting g p isn't 
feasible. We point out though that total OMC rates are 
employed in analyzing the total rate of radiative capture 
on complex nuclei (see Sec. |x|). 



A. Free couplings versus effective couplings 

We first wish to comment on the meaning of the cou- 
pling g p that's extracted from partial transitions on com- 
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plex nuclei, and specifically the issue of free nucleon 
coupling constants versus effective nucleon coupling con- 
stants. 

Suppose one performed an exact calculation of the 
initial-final wave functions and the weak transition, that 
is, incorporating all components that contribute signifi- 
cantly to the wave functions and all currents that con- 
tribute significantly to the transition. Naturally it would 
require a large number of basis states, a realistic interac- 
tion between nucleons, and the contribution of nucleonic 
excitations, such as the A resonance, and exchange cur- 
rents, such as pion exchange. However such a calculation 
would involve the free nucleon values of the weak cou- 
plings g a , g p , etc. 

However, in practice only for few-body systems are 
such highly accurate calculations feasible. Rather in 
complex nuclei a common practice is to consider only a 
handful of a.ctive nucleons and a, truncated model spa.ee, 



targets we shall denote the capture rate by A and for 
Ji > targets we shall denote the hyperfine rates by A± . 
In addition on J,; > targets it's useful to define the 
statistical rate 



As = (4 ±i 7)A + + ( T A r)A- 



'2J t 



V 



■h + r 



(32) 



where (Ji+l)/(2Jj+l) and Jj/(2Ji+l) are the statistical 
populations of the F± hyperfine states. 



2. Recoil asymmetries 

Generally the direction of emission of recoils after cap- 
ture is anisotropic about the (X— atom orientation axis. 15 
For Ji = targets the ang u lar distribution for recoil emis- 
sion i s (Bcrnabcu , 1975|; |Cicchanowicz and Ozicwicz 



and I to omit both exchange currents and nucleonic exci- 
tations. In such calculations the free coupling constants 
are replaced by effective coupling constants, with their 
renormalization accounting for the overall effects of the 
missing components of the nuclear model. The use of ef- 
fective electric charges and an effective axial coupling are 
well known examples of this approach to the problem. 

In this article our interest is in the coupling g p of the 
free nucleon. We therefore have tried to use the most so- 
phisticated models available, in order to extract the cou- 
pling g p . In almost every case we discuss there are full- 
shell model calculations with perturbation theory com- 
putations of core polarization effects, i.e. the effects of 
truncating the model space, and exchange current effects, 
e.g. the effects of pion exchange and A excitations. In 
this way we believe we come as close as possible to ex- 
tracting the free nucleon value for the induced coupling 
9v 



1984 [Mukhopadhyay , 1977j ~ 



— = — [1 

dn 4tt l 



a Pi(P M -z>)] 



(33) 



where P M is the muon polarization, v is the neutrino di- 
rection, Pi is the £ = 1 Legcndre polynomial, and a is 
the asymmetry coefficient. Note that the recoil asymme- 
try is a pseudoscalar quantity, and therefore an example 
of parity violation in muon capture. 

For Ji > targets the angular distribution for recoil 
emission is generally more complicated. For example a 
F = 1 /i atom may possess both a vector polarization, 
i.e. non-statistical populations of the mp = +1 and 
rap = — 1 sub-states, and a tensor polarization, i.e. non- 
statistical populations of the \mp\ = 1 and mp = 
sub-states. Consequently for F = 1 capture the angu- 
lar distribution of recoil emission is (piechanowicz and 



B. Physical observables 

In muon capture the initial state is the IS ground state 
of the muonic atom. Capture then yields a muon neutrino 
with momentum v and recoil nucleus with momentum 
— v, i.e. 



Oziewicz , 1984 ; Conglcton and Fearing , 1993; Galindo 
and Pascual , 1968| J 



dA A r , , -» , „ . a 
- = -[l + ^ |P,|P l( P, 



A t |P M | P 2 (P M -z>)] 
(34) 



[A,Z] — ► [A,Z-1] 



(31) 



For Ji — targets the total spin of the /z-atom is 
F = 1/2, and for Ji > targets the possible spins of 
the /j,-atom are F± — Ji ± 1/2. On Ji = targets the 
necessary variables for describing capture are the /z-atom 
orientation, the recoil orientation and the neutrino direc- 
tion. On Ji > targets the capture process is further 
dependent on the P-state. 



where A v is the vector asymmetry coefficient and A t is 
the tensor asymmetry coefficient. Note for F > 1 atoms 
even higher ranks of /it-atom orientation are possible in 
principle. 

We stress that P M in Eq. (||) and Eq. @ is the muon's 
residual polarization in the /u-atom ground state. Un- 
fortunately in all targets the muon nucleus spin orbit 
interaction leads to substantial depolarization during 
the atomic cascade. For example for J; = targets 
the residual polarization is typically 17%. Further in 
J =/= targets the muon-nucleus spin-spin interaction 



1. Capture rates 



The easiest observable to experimentally determine is 
the capture rate of the partial transition. For Ji = 



15 An exception is the isotropic distribution of the recoil emission 
from a F = /i-atom. 
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leads to additional depolarization during the atomic cas- 
cade. The degree of depolarization is dependent on J, 
and F. For further details see |Favart et al. (1970 ]j_ Ham- 

~~|, 11977| ' 



bro |md Mukhopadhyay (1975| , ^L977|) ; |Kuno et al. (1987|) 



Mukhopadhyay (1977). 



3. Recoil orientations 



recoil asymmetry in the capture process. For unpolarized 
muons a 7-1/ directional correlation is possible and for po- 
larized muons a 7-1/-// triple correlation is also possible. 
In discussing the correlations we shall denote the spin 
sequence in Eq. (^) by J; — > Jf — ► j where J^, J/ and 
j are the angular momenta of the three nuclear states. 
Note that below we consider only unpo larized targets and 



Jj_ < 1 recoils. For further details see Ciechanowicz and 



Generally after muon capture the recoil nucleus is ori- 
ented, where orientation along both the i/-momentum 
axis and the /x-spin axis is possible. For a Jf — 1/2 recoil 
the orientation is described by its recoil polarization P, 
where P = {p +1 / 2 ~ P-1/2)/ J^Pm with p m denoting the 
population of each magnetic sub-state m. For a Jf = 1 
recoil the orientation is described by its polarization P 
and its alignment A, where A = (p+\ +p-\ — 2po)/^2pm- 
For Jf > 1 recoils even higher ranks of orientation are 
possible in principle. Conventionally the recoil orienta- 
tion along the ^-momentum axis is designated the longi- 
tudinal polarization and longitudinal alignment A^, 
and the recoil orientation along the /i-spin axis is des- 
ignated the average polarization P av and average align- 
ment A av . The recoil polarizations, i.e. Pl and P av , are 
pseudoscalar quantities and additional examples of parity 
violation in m uon capture. For further details see |Bern- 



Oziewicz (1984) 



abe u (1975 ); Devanathan et al. (1972 ); Mukhopadhyay 
(197 7|) ; [Bubramanian et al. (1976| ) 



Additionally triple correlations of the recoil orientation 
with the /x-spin direction and the i/-momentum direc- 
tion are possible. Of experimental importance are the 
so-called forward hemisphere and backward hemisphere 
polarizations. Defining the hemispheres relative to the 
/i-spin axis, Pf is the polarization along the /x-axis for 
the 'forward hemisphere' recoils, and Pb is the polar- 
ization along the /Lt-axis for the 'backward hemisphere' 
recoils. Note that P^ and Pb are related to Pi and P av 
via 



P F = \{P a 



Pb — -^{Pav - -^Pl) 



(35) 



(36) 



which demonstrates their dependence on both the ^(-spin 
direction and the ^-momentum direction. 



4. Gamma-ray correlations 

Gamma-ray directional correlations with the v- 
momentum axis and the /z-spin axis in the sequence 



[Z,A] 



\Z-1,A]* 



[z-Mr+7 

(37) 



are additional observables in partial transitions. Such 
correlations originate from the recoil orientation and the 



For unpolarized muons P^ 
relation is given by 16 



the 7-1/ directional cor- 



W = \ + a 2 P 2 {k ■ v) 



(38) 



where k is the 7-ray direction, v is the neutrino direction, 
P 2 is the i = 2 Legendre polynomial, and a 2 is the j-v 
correlation coefficient. Note that this directional corre- 
lation a 2 is a consequence of the longitudinal alignment 
Al of the recoil nucleus in the capture process. 

For polarized muons P M 7^ the j—u—p, triple correla- 
tion is given by 



W = 1 + (a+ §ci)P M • v k ■ v 
+ (a 2 +b 2 P^ ■ v k ■ v)P 2 (k ■ v) 



(39) 



and contains three distinct correlation terms involving 
(a + |ci), a 2 and b 2 . Note in Eq. ( gg ) that (i) k en- 
ters only quadratically because of parity conservation in 
7-decays, (ii) Pp enters only linearly because the muon 
is spin-1/2, and (iii) v enters in powers of two or less 
for Jf < 1. In Eq. ([39]) the correlation involving a is a 
consequence of the recoil asymmetry and the correlation 
involving a 2 is a consequence of the longitudinal align- 
ment. The remaining terms, i.e. c\ and 62, originate 
from the triple correlation of the recoil alignment with 



If either unpolarized muons, i.e. P M 



0, 



P M and v. 

or perpendicular geometry, i.e. P M • k — 0, is employed 
the only non-vanishing correlation is a 2 , and Eq. ( |39| ) 
becomes Eq. (|38|). 

It's important to recognize that c±, a 2 and b 2 are func- 
tions of both the /i-capture process and the 7-decay pro- 
cess. For example the coefficient a 2 may be written in the 
form of a product AlB 2 , where the alignment Al is gov- 
erned by the \x capture process and the parameter B 2 is 
governed by the 7-decay process. A handy compilation 
of the coefficients B 2 for various Jf — ► j spin-parities, 
multipolarities, and mixing ratios is given in Table 1 of 
Ciechanowicz and Oziewicz (1984] ). Note that in cer- 
tain cases, e.g. for Ml emission in a 1 + — > + decay, 



16 The theory nf -v-rav correlations was developed hy Popov and 
co-wor kers in Buk at and Popov (1964); FSukhy o stov and Popo 1 



7|Parthasi 



(19678,44 l97fH - 3zi£3d£Z and Popov Popo v GSffi j 



See alsoH'arthasarathy and Sridhar (1978 ) and |Cicchanowicz and 
Oziewicz (1984[ ). 
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£?2 is large, but in other cases, e.g. for Ml emission in 
a 1 + — > 2 + decay, B 2 is small. This makes the former 
1 + — > + case more favorable, and the latter 1 + — > 2 + 
case less favorable, for 7-ray correlation experiments. 



C. Helicity representation 

Herein we consider the overall dynamical content of 
exclusive muon capture. Utilizing the helicity represen- 
tation we discuss the constraints imposed on [i capture 
by ^-handedness and T-invariance. We shall denote the 
corresponding helicity amplitudes by T^, where A is the 
recoil helicity and F is the /x-atom hyperfine state. We 
use the superscript '+' for F = Ji + 1/2 capture, the su- 
perscript '— ' for F — Ji — 1/2 capture, and no superscript 
for Ji — atoms. Assuming the absence of T violation 
in fi capture the helicity amplitudes are real numbers. 

Further details on the helicity represen t ation in the 
(fx, v) reactio n are given in Bernabeu (1975), [Mukhopad- ' 



hya f (1977j ) and Ciechanowicz and Oziewicz (1984) 



1. Capture on zero-spin targets 

For Ji = targets the helicity representation depicts 
muon capture as the two-body decay of a spin-1/2 (j,- 
atom into a left-handed neutrino and a spin- J/ recoil. 
Choosing the z-axis along the v-axis, the definite neu- 
trino helicity of = —1/2 means the allowable recoil 
helicities are Xf = 0, +1. The corresponding helicity 
amplitudes, denoted by To and T±, are the underlying 
dynamical variables in fi capture on Ji = targets. 

In Table ^ we compile explicit formulas for various 
physical observables in — * Jf transitions in terms of 
T , Ti and their ratio X = ^2T Q /T X . The 7-ray cor- 
relation coefficients also involve the quantity B2 that is 
governed by the 7-decay. With the exception of A, the 
reaction dynamics are completely determined by X , and 
con sequently there exist numerous relations between ob- 



TABLE IV Helicity decomposition of physical observables in 
— ► 1 transitions. We give the capture rate A, recoil asym- 
metry a, longitudinal polarization and alignment Pl and Al, 
average polarization and alignment P av and A av , and 7-ray 
correlation coefficients d2, b 2 and (a+ fci). The helicity am- 
plitudes for recoil helicities of A/ = 0, +1 are denoted by To 
and Ti and X = %/2To/Ti. The parameter B2 involved in the 
7-ray correlations depends on the spin-parities of the initial- 
final states and the multipolarity of the gamma-radiation in 
the 7-decay. 



Observable 



Helicity decomposition 



Pl 
A l 
p 

J av 
A a v 



t>2 



\T 1 \ 2 (2 + X 2 ) 
(X 2 -2)/(X 2 + 2) 

-2/(2 + X 2 ) 
2(1 — A 2 )/(2 + X 2 ) 
2/3(l + 2A)/(2 + A 2 ) 


B 2 V2{2 + 2X - X 2 )/(2 + X 2 ) 

B 2 V2(l-X 2 )/(2 + X 2 ) 
B 2 V2(l -2X + X 2 )/(2 + X 2 ) 



lar momentum coupling the total number of recoil helic- 
ity states Xf and contributing helicity amplitudes Tf is 



serv ablcs in — > Jf transitions, for example sec Bernabeu Oziewicz , 1984; Mukhopadhyay , 1977) 



(197|5|) and [Mukhopadhyay (1977| ). Clearly although dif- 



thc lesser of cither 2 J f + 1 or 2T + 1 (|Cicchanowicz and 



ferent observables offer alternative possibilities for exper- 
imental measurements, the dynamical content of — > Jf 
transitions is somewhat limited. 

Note that the 0^0 sequence is a special case. Com- 
pared to Jf > 0, where the Xf = 0, +1 helicity states are 
allowable, for Jf = 0, only the Xf = helicity state is 
possible. Consequently a single amplitude To is the sole 
dynamical quantity for physical observables in — » 
transitions. 



2. Capture on nonzero-spin targets 

For Ji > targets the helicity representation depicts 
capture as a two-body decay of the F = Ji ±1/2 di- 
atom. Based on definite neutrino handedness and angu- 



For concreteness let's consider the example of a 1/2 — > 
1/2 transition, where the /it-atom spin is either F = 
or F = 1. For F — 1 both the Xf = -1/2, +1/2 re- 
coil helicity states are populated, but for F = only 
the Xf = +1/2 recoil helicity state is populated because 
of the single magnetic sub-state of the F — ^t-atom. 
Consequently one helicity amplitude governs the F = 
capture, two helicity amplitudes govern the F = 1 cap- 
ture, and three dynamical variables underlie 1/2 — > 1/2 
transitions. By comparison in 3/2 — > 1/2 transitions a 
total of four helicity amplitudes contribute, two for F + 
capture and two for T_ capture. 

In summary for Ji > targets many independent 
quantities are experimentally accessible. In principle the 
increased number of variables allows an increased number 
of cross checks on model calculations for J, > transi- 
tions. 
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D. Induced currents 

Most often in partial transitions the leading contribu- 
tions originate from the axial coupling g a . This subsec- 
tion concerns 'where to find the coupling g p T 

To answer this question we shall examine capture in 
(i) the q/M — » limit, where q is the three-momentum 
transfer and M is the nucleon ma ss, and (ii) in the Fujii 



Primakoff (FP) approximation (Fujii and Primakoff 
195!fl). In the q/M — ► limit the effects of g p are ab 
sent. However in the FP approximation, which keeps 
q/M terms involving allowed operators, but drops l- 
forbidden and gradient operators, the leading effects of 
induced currents are present. Therefore comparison of 
the capture process in the q/M — > limit and the FP 
approximation is helpful in understanding the manifes- 
tation of g p in capture. 

At this point it is helpful to recall the Fujii- Primakoff 



effective Hamiltonian for muon capture (Fujii and Pri- 
malJoff , I959| [Primakoff , 19591) . Its form is 



T+ - — y~ - J2 T * ( Gv 1 ' Xi + Ga a ' CTi 

8=1 

+G P cr u tn-fi) S(r - n) (40) 

where 1, li, cr and er^ are the 2x2 unit and spin matrices 
and r and r j are the spatial coordinates of the lepton and 
the i th nucleon respectively, is the //-momentum unit 
vector, and r + , tT convert the muon into a neutrino and 
a proton into a neutron. Last, Gy, Ga and Gp are the 
so-called Fujii-Primakoff effective couplings 



G\ 



g (1+ JL) 

y y 2M' 



(41) 



the ^-momentum axis. However this limitation does not 
apply to the Ga term and the operator (1 — cr ■ C>) cr ■ er^. 
Consequently in Eq. the Gp term admits only a 
longitudinal coupling of the weak currents whereas the 
Ga term admits also a transverse coupling of the weak 
currents. This makes correlations, orientations, etc., sen- 
sitive to g p /g a - For further details see prenacs' (1985| ). 

To further illustrate the g p sensitivity we consider the 
example of a + — > 1 + transition. In the q/M — > limit 
the two helicity amplitudes, i.e. Tq and Ti, are both 
determined by the product of the coupling g a and the 
allowed Gamow-Teller operator, and X = 1. However 
in the FP approximation, the amplitude To is propor- 
tional to (Ga — Gp) and involves a longitudinal coupling 
of the weak currents, and the amplitude T\ is propor- 
tional to G a and involves a transverse coupling of the 
weak currents, and X — (Ga ~ Gp)/Ga — 0.59 for 
q = 100 MeV/c. Therefore the asymmetries, correla- 
tions and orientations, which are governed by X, per- 
mit the determination of g p /g a - For further details see 
Mukhopadhyay (1977| ). 

In Table |V| we compile FP approximation expressions 
for physical observables in + — > 1 + transitions in terms 
of Ga and Gp. Note in the q/M — > limit the recoil 
polarizations are Pl — —2/3 and P av = +2/3 and the 
recoil alignments are Al = and A av = 0, thus indicat- 
ing that the recoil is 'highly polarized' but 'not aligned'. 
Most strikingly, in the FP approximation the longitudinal 
alignment Al = +0.55 is very large and highly sensitive 
to g p /g a - This occurs because the alignment measures 
the difference in population of the Ay = recoil sub- 
state, which is populated by Gp, and the A/ = 1 recoil 
sub-state, which is not populated by Gp. It therefore 
represents a golden observable for the spin structure of 
the induced coupling g p . 



Ga = -{g a + 2j7(9v + 9m)) 



(42) 



2. Hyperfine dependences 



G P = 



' 2M ^ 9p 



g a + 9v + g m ) 



(43) 



Note that, when q/M — ► 0, the coupling Gp vanishes and 
Gy and Ga are determined by g v and g a respectively. 
Furthermore the induced currents are order q/M, and g p 
appears in Gp only. For the canonical values of the weak 
couplings, as given in Table |lj (see Sec. 0), G v = 1.00, 
G A = -1.27 and G P = 0.00 in the q/M -> limit, and 
G v = 1-03, G A = -1-52 and G P = -0.62 in the FP 
approximation and with q = 100 MeV/c. 



1. Asymmetries, orientations and correlations 

To understand the g p sensitivity of orientations, cor- 
relations, etc., we compare the terms with Ga and Gp 
in Eq. ([l0"|). Observe the operator corresponding to the 
Gp term, i.e. (1 — cr • u) cr ■ cri ■ 0, cannot change the 
spin projection of either the lepton or the nucleon along 



The hyperfine dependence A + /A_ of AJ T = Jf — J, = 
±1 + transitions in muon capture on Ji ^ targets is also 
sensitive to the coupling g p . Specifically in the multipole 
expansion of the Fujii-Primakoff Hamiltonian the G a~ 
term makes allowed contributions to neutrino waves with 
total angular momentum = l/2 + only, whereas the 
Gp-term makes allowed contributions to neutrino waves 
with total angular momentum j n = 3/2 + also. This 
difference is because of hatu in Eq. ((40|). In A = +1+ 
transitions, neutrinos with j = 1/2^ may be emitted 
in F + capture, but neutrinos with j = 3/2 + must be 
emitted in i+ capture. This makes A + /A_ >> 1 and a 
strong function of the ratio g p /g a - However in AJ*" 

I transitions, neutrinos with j = 1 /2 + may be emitted 
in F_ capture, but neutrinos with j = 3/2 + must be 
emitted in F + capture. This makes A + /A_ << 1 and a 
strong function of the ratio g p /g a - 

For example, we consider the important case of 
3/2+ — > 1/2+ transitions. In the q/M limit, cap- 
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TABLE V Expressions and values of observables for + — ► 1 + 
transitions in the q/M — * limit and the Fujii-Primakoff ap- 
proximation. Ga and Gp are the Fujii-Primakoff effective 
constants defined in Eqs. and (^3|) and the common de- 
nominator is r = (3G A + G P — 2GaGp). We tabulate the 7- 
ray directional correlations for the + — ► 1 + — > + sequence. 



Observable 


FPA 


q/M -f 


FPA 




Eq. 


value 


value 


a 


(3Gi + Gp - 2G A Gp)/r 


-0.33 


-0.70 


Pl 


-2Gi/r 


-0.67 


-0.85 


Ai 


(-2G 2 p +4G A G P )/r 


0.00 


-0.55 


av 


1 c\/~i2 4 /~i /~1 \ /"n 
(2G A — gLrALrpj/l 


+0.(3 < 


+0.61 













OL + |C1 


(3Gi - G P )/r 


1.0 


1.21 


a 2 


(-Gp + 2GpGA)/r 


0.00 


0.28 


62 


G 2 p/r 


0.00 


0.07 



ture from the F— = 1 hypcrfinc state is governed by g a 
but capture from the F + = 2 hyperfine state is forbid- 
den. Therefore A + /A_ = 0. However in the FP approxi- 
mation, the two F + helicity amplitudes are proportional 
to Gp while the two F- helicity amplitudes are propor- 
tional to either (AGa — Gp) or (2Ga — Gp). Consequently 
the hyperfine dependence A + /A_ is highly dependent on 

9p/9a- 

The FP approximation expressions for various observ- 
ables in 3/2+ — ► l/2 + transitions are given in Table [v| 
In passing we mention that the asymmetries, correlations 
and orientations in F„ capture for 3/2+ — > 1/2+ transi- 
tions are also sensitive to the coupling g p , for the same 
reasons as described for the + — > 1 + transitions in the 
This is not true for the case of 



preceding Sec. IX.D.l 



asymmetries, correlations and orientations for F+ cap- 
ture in 3/2+ — » 1/2+ transitions, because only j 77 = 3/2+ 
neutrino emission is possible. 



Related FP 
AJ 71 = ±1+ 



approximation 
transitions are 



Mukhopadhyay 

(pair 



expressions for other 
given for example by 



TABLE VI Equations and values for selected observables in 
3/2+ -> 1/2+ transitions in the q/M -> limit and Fujii- 
Primakoff approximation. The common denominator is V — 
(24Gi - 16GaGp + 3Gp). 



Observable 


FPA 
Eq. 


q/M -» 
value 


FPA 
value 


A+/A_ 


3Gp/r 


0.00 


0.027 


a~ 


-i(4GA-Gp) 2 /r 


-0.22 


-0.24 


a+ 


2 
5 


+0.60 


+0.60 


Pl 


(8Gi - G 2 P )/r 


+0.33 


+0.44 


pi 


1 

5 


+0.20 


+0.20 



3. Capture rates 

In general in fi capture the rate A has a fairly weak 
dependence on the coupling g p . For example in + — > 
1 + transitions the contribution of g p is typically 10%. 
However an exception is the capture rate of the first- 
forbidden + 0~ transition. 

To understand the sensitivity it is instructive to as- 
sume the dominance of the 1=1 retarded Gamow- Teller 
operator in the first-forbidden + — > CP transition, i.e. 
ignoring the contribution from the axial charge operator. 
For details see Towner and Khanna (1981). The capture 
rate in + 



(1977) and Ciechanowicz and Oziewicz 



0~ transitions is then governed by the cou- 
pling constant combination {Ga — Gp), which is strongly 
dependent on g p . For example in going from g p = to 
g p = 6.7g a the rate is increased by roughly 50%. In short, 
the quantum numbers of the A J 71 " = 0~ multipolc are ef- 
fective in isolating the longitudinal contributions to weak 
currents, such as g p . 



E. Experimental studies of partial transitions 

This subsection concerns the experimental work on 
exclusive OMC. We discuss two allowed transitions on 
Ji = targets, 12 C(0+,0) -» 12 B(l+,0) and 28 Si(0+,0) 
— > 28 A1(0 + ,2201), two allowed transitions on Jj + 1 tar- 
gets, n B(3/2-,0) -> n Be(l/2-,320) and 23 Na(3/2+,0) 
-> 23 Ne(l/2 + ,3458), and the first-forbidden transition 
16 O(0+,0) -> 16 N(0", 120). The experiments include 
measurements of capture rates, recoil orientations, 7-ray 
correlations, and hyperfine dependences. 
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1. 12 C(0+,0) -> 12 B(l+,0) 

The 12 C(0+,0) -> 12 B(l+,0) reaction is an allowed 
Gamow- Teller transition between the spin-0, isoscalar 
12 C ground state and the spin-1 isovector 12 B ground 
state. T he transition wa s first observed in cosmic-ray 
data by Godfrey (1953 ) via the identification of the 



backgrounds. Transverse-field precession of the /i spin is 
used to measure the muon polarization P M . 

One challenge is the experimental determination of a 
small /3-ray asymmetry (~3%) with a reasonable accu- 
racy (±10%). Consequently false asymmetries, such as 
geometrical and inst rumental effects, must be minimized 
and then measured. Possoz et al. (1977, 1974) used a po- 



Godfrey— Tiomino cycle, i.e. fx capture on C followed 
by /3-decay of 12 B. 

The first investigation of the recoil polarization in 
the 12 C(0 + ,0) ^ 1 2 B(l + ,0) transition was conducted by 
Love et al. (1959) in order to measure the fj,~ helicity 



in 7r decay. The application of polarization measure- 
me nts to induced currents was pioneered by [Possoz et 



larization preserving target material (graphite) and a po- 
larization destroyi ng target (polyethylene) to determine 
false asymmetries. Kuno et al. (1984 , 1986| ) used a novel 
magnetic resonance technique to periodically destroy the 
muon polarization. 

In summary the resulting values of th e average po- 



al. (119771. Il974|) at Saclav and extended by [Rocsch ei (1974| ) 
B a 



larization were P av = 0.38 ± 0.07 from Possoz et al. 
P av = 0.452 ± 0. 042 from possoz et al. (1977|) 



al. |l981a| ,Tb|) and |Truttman et al. (1979|) at PSI. Us- p av 0.462 ± 0.053 from |Kuno et al ( 1984| , |1986|) . We 



ing ingenious techniques these researchers have amassed 
impressive data on several polarization observables in 
12 C (0+,0) -^ 12 B(l+,0) capture. More recently |Kuno ei 
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al ( |L984 11986| ) at the BOOM facility have polished some 
techniques and re-measured the polarization P av . 

The procedure for measuring the polarization P av of 
12 B recoils from 12 C capture is straightforward in princi- 
ple. First one makes polarized /z _12 C atoms by stopping 
polarized muons in carbon-containing material. Next 
the polarized /i 12 C atoms produce the polarized 12 B re- 
coils via muon capture. Then one measures the 12 B (3- 
ray asymmetry to determine the recoil polarization P av - 
Note that the method relies on the known correlation of 
the ejected /3-rays with the 12 B orientation. Also note 
that this method gives the average polarization P av , i.e. 
along the /Lt-axis, not the longitudinal polarization P^, 
i.e. along the v-axis. 

Unfortunately the measurement of P av is difficult in 
practice. First the muon polarization is largely destroyed 
via the spin-orbit interaction in the atomic cascade. Sec- 
ond the 12 B polarization is easily destroyed by the spin- 
spin interaction in the host material. Third only about 
1% of n sto ps in 12 C undergo th e 12 C(0+, 0) -^ 12 B(1+, ) 



transition (Maier et al, 1964; Reynolds et al, 1963). 
Consequently, the beta-ray asymmetries are small, the 
target choices are limited, and backgrounds are trouble- 
some. Fortunately the short 29 ms lifetime and high 
15 MeV end point for 12 B beta-decay are ideally suited 
to a measurement of the polarization. 

Average polarization measurements have been per- 
for mcd by Possoz et al. (1977 , 1974) at Saclay and Kuno 
et a [ (19847p"986| ) at BOOTVL The basic set-up comprises 
a beam telescope for detecting incoming muons, beta- 
ray counters for detecting 12 B decays, and Miche l coun- 
ters for detecting /x-decays. Possoz et al. (1977 , 1974 ) 
found a >0.3 kG longitudinal B-field and a graphite 
target were sufficient to preserve the 12 B polarization 
for t > 29 ms. Note that the distinctive lifetime/end 
point are used to identify the 12 B /3-rays and their for- 
ward/backward count rates are used to determine the 
12 B /3-asymmetry. Pulsed-beam operation allows beta- 
ray detection under beam-off conditions and reduces the 



discuss the extraction of g p from P av in Sec 

Additionally the forward (Pf) and backward (P_b) 
polarizations, defined in Eqs. ( |35"| ) and j3^), for the 
12 C (0 + ,0) — > 12 B(l + ,0) tr ansiti on have been measure d 



by [Truttman et al. (1979| ) and [Roesch et al. (1981b[ ) 
These experiments are masterpieces of ingenuity and 
technique. 

The experiments employed a novel target with recoil 
direction sensitivity. The targets comprised a multi-layer 
sandwich of 'triple-foils' with each triple-foil comprising 
a carbon target foil (C), a polarization preserving foil 
(P), and a polarization destroying foil (D). Arranging 
the P— foil upstream and D-foil downstream of the car- 
bon foil, i.e. DCP, permits selective retention of the for- 
ward hemisphere recoil polarization. Arranging the D- 
foil upstream and P-foil downstream of the carbon foil, 
i.e. PCD, permits selective retention of the backward 
hemisphere recoil polarization. Thereby /3-ray asymme- 
try measurements for the DCP configuration yield Pf 
and the PCD configuration yield Pb- Additionally a 
PCP sandwich enables the measurement of the average 
polarization and a DCD sandwich enables the measure- 
ment of the false asymmetries. 

One advantage of measuring P f/b is a larger recoil po- 
larization and a larger beta-ray asymmetry. 17 Another 
advantage is that combining both Pf and Pb to deduce 
Pav/PL where 



av _ 2 r F 



Pl 



Pf-Pb 



(44) 



reduces the sensitivity to false asymmetries and system- 
atic uncertainties. Note that the measurements of Pf 
and Pb are achieved by simply rotating the multi-layer 
target by 180 degrees. 

However a major challenge in measuring Pf/b is the 
small quantity of the 12 C material in the multi-layer tar- 



17 Recall from Sec. 1X.B.3 that Pp/s have contributions from both 
Pav, which is decreased by fj, depolarization, and from P^, which 
is unchanged by fj, depolarization. 
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get, since the carbon foils must be thin enough for recoils 
to emerge and the P /D foils must be thick enough for re- 
coils to stop. Consequently any 12 B background from 
stops in nearby carbon is es pecially perilous. The refore 
Truttman et al. (1979) and Roesch et al. (1981b) took 
great care to avoid using any carbon-containing materi- 
als in the neighborhood of the experiment. Based on 
these data the authors obtained P av /PL = —0.516 ± 
0.041. 



In Table |VII we summarize the various measurements 

12 B(l+,0). Note 



of recoil polarizations in 12 C(0 + ,0) 
we quote all results in terms of the helicity amplitude 
ratio X so as to facilitate their comparison. The experi- 
mental results are mutually consistent. 

We stress the experimental results for P av and P F/ / B 
are for muon capture to all bound states in the 12 B nu- 
cleus. Therefore a correction is necessary to obtain the 
interesting ground state polarization from the observed 
bound state pola rization. Measurements of the total rate 
to bound states ( Maier et al., 1964 ; Reynolds et al, 1965 ) 
and the individual rates to e xcited states (Budyashov ei 
L970|; piffon et al, 198l|; |Miller et al., 1972a|; |Roesch 



et a 



, 1981a) are employed to determine this correction. 



Unfortunately the 7-ray measurements of fi capture 
to individual 12 B excited states are complicated by near 
equal energies of several Doppler broadened 7-rays. 18 



Furt hermore the latest rcaulta of Rocach et al. (1981a 



and I earlie r roou lto of - Budyaohov ct al. (197C )y- Miller m 



al. Ql972a| ) and piffon et al. (198l| ) are in disagreement 
and consequently there exists some uncertainty in the 
correction for the capture to the 12 B excited stat es. For 
a d etailed d iscussion of the experimental d ata see [Meas- 
day (2001 ). Using the 7-ray yield s from Giffon et al. 



(195 lj) and the model calculations of Fukui et al. (1987) 
the world average for X in 12 C(0+,0) ->■ i2 B(l +, 0)is 



0.23 ±0.06. H owever using the 7-ray yields from Roesch 
et a l. (1981a ) and the model calculations of [Fukui ei 
al. 1987| ), the world average for X in 12 C(0+,0) -> 
12 B(1 + . 0) is 0.20 ± 0.06. Fortunately the corrections are 
not too large. 



Be and 23 Na 



! Ne 



Historically studies of the hyperfine effect in nuclear 
muon capture reaction were important in demonstrating 
the V — A structure of weak interactions. These exper- 

iments wer e performed in m uonic i9 f by [^nn gan et al 



(19^1]) and |Winston (1963| ). The first application of the 



hyperfine effect to the induced coupling g p was conducted 
by [Deutsch et al. (1968 ) in muonic n B. Recently Wiaux 
et a . (2002) at PSI have improved the data on n B, and 



Johnson et al. (1996| ) at TRIUMF have extended the data 



IX.B.l that the IS ground state of 



to^Na. 

Recall from Sec. 
a Ji 7^ /i-atom is a hyperfine doublet with a spin 
F± = Ji ± 1/2. These states are split by the spin-spin 
interaction of the muon-nucleus magnetic moments. For 
a positive nuclear magnetic moment the F- state is the 
true atomic ground state and for a negative nuclear mag- 
netic moment the F + state is the true atomic ground 
state. 

The two hyperfine state are initially populated with 
statistical weights, i.e. n + = (Ji + l)/(2Jj + 1) and 
n_ = Jj/(2Jj + l). 19 Thereafter hyperfine transitions 
from the upper F-state to the lower F-state occur by 
Mi-Auger emission, changing the relative occupancies 
of the hyperfine states. The rate of hyperfine con- 
version by Auger emission is governed by (i) the wave 
function overlap of the electron and the ^t-atom and (ii) 
the relative sizes of the electron binding and the hyperfine 
splitting. The wave function overlap leads to a system- 
atic increase in A^ with Z whereas the electron binding 
leads to sudden decreases in A/j at Z ~ 6, where K- 
shell emission is forbidden, and Z ~ 18, where L-shell 
emission is forbidden. A detailed account of hyperfine 
conversion by Auger emission is given by Winston (1963) 
and experimental determinations of conversion rates have 



generally confirmed the calculated rates (Gorringe et al. 



1993J; |Measday , 200 1| ; [Suzuki et al., 1987| ). Most im 



portantly for investigating g p , there exist a handful of 
/i-atoms with comparable hyperfine transition and muon 
disappearance ra tes, e.g. n B, 19 F, 23 Na and 35 C1. 

iciegcli (iyt>yj was first to recognize that the hyper- 
fine conversion during the /i-atom lifetime would alter 
the muon occupancy of the F-states, and modify the 
time spectrum of the capture products, and thus permit 
the determination of the hyperfine dependence A + /A_. 
Formulas for the time evolution of capture products in 
the presence of hyperfine conversion have been published 
by several authors. For the case of a positive nuclear 
magnetic moment, e.g. n B and 23 Na, 



dN 
~dt 



oc 



-A D U(n , Ah \ 

A+ A ft \„-(A h +AA D )tl 



e 

+ (A± _ _ 



-AAj 



(45) 



where n + , n and A+, A_ are the initial populations and 
the capture rates of the F+, F_ hyperfine states, A^ is 
the F- state disappearance rate, AA_d = A^ — A^ is the 
hyperfine disappearance increment, and A^ the hyperfine 
conversion rate. In almost all cases of interest A/j >> 
A Ad and therefore the factor Ah /(Ah ± AAd) is close to 
unity. 



18 Specifically (i) the 947 keV (2620 1674) and 953 keV (953 -> 
0) gamma-ray lines and (ii) the 1668 keV (2620 — > 953) and 
1674 keV (1674 — > 0) gamma-ray lines 



The exception is the use nf a polarized target and a polarized 



beam (Hambro and Mukhopadhyay 
197?| ). 



1977 



Mukhopadhyay 
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TABLE VII Compilation of results from the recoil polarization experiments for the transition 12 C(0 + ,0) — > 12 B(l + ,0). The 
results are presented in terms of the dynamical parameter X = v2 Tp/Tx. They incorpor ate the corrections fo r capture to 12 B 
bound excited states using the 7-ray data from Giffon et al. (1981) in column three and (Roesch et al, 1981a) in column four. 
The experimental observable is listed in column two. 



Ref. 



obs. 



X using 



Siffon et al. (1981) 



X using 



Roesch et al. (1981a 



Possoz et al. (1974) 



Possoz et al. (1977) 



Roesch et al. (1981b) 



Kuno et al. (1984| ) 
FPA 



p 

i av 
Pav 

Pf/b 
p 

J av 



0.10 ±0.11 
0.22 ±0.07 
0.27 ±0.07 
0.23^ 
0.59 



1-0.10 
-0.08 



0.08 ±0.11 
0.20 ±0.07 
0.24 ±0.07 
0.21^ 
0.59 



, +0.10 

-0.08 



As discussed in Sec. IX. D. 2 the hyperfine dependence 

±1 + spin- 
[3) demon- 



A + /A_ of partial transitions with AJ X : 
sequences is especially sensitive to g p . Eq. 
strates clearly that the hyperfine dependence A + /A_ in 
such capture is encoded in the time spectrum of the re- 
action products 



T fre ll B(3/2~,0)-> ll Be(l/2~, 320) reaction is a 
AJ W = 1 + allowed transition from the n B ground state 
to the 320 keV n B first-excited state. The 320 keV, 1/2" 
11 Be excited state decays by 7-emission (r = 115 ±10 fs) 
to t he n Bc ground state, and the 11 Bc l/2 + ground 



pcriment. 

Both experiments fit their 320 keV 7-ray time spectra 
to the time dependence of Eq. ( f45f ) in order to extract 
A + /A_. The instrumental time resolutions were deter- 
mined via muonic x-rays and the continuum backg rounds 
were subtracted via neighboring energy windows. Wiaux 



et al. (2002) analyzed both the 320 keV 7-ray time 



spectrum and the Michel time spectrum in order to re- 
duce the correlations between A^,, A^ and A + /A_. Fo r 



n B(3/2-, 0) -> n Be(l/2~, 320) the |Deutsch et al. (196E 



experi ment yielded A + /A_ < 0.26 and the Wiaux et al 



staf f decays by /3-emission (r 

n B ground state. For reference, the 11 B disappear- 



13.81 ± 0.08 s) to the (2002| ) experiment yielded A + /A_ = 0.028±0.021. These 



anc e rate is A D = (0.4787 ± 0.0008) x 10 6 s" 1 (Wiaux 



et < \l, 2002 ) and the n B hyperfine transition rate is 
A h = (0.181 ± 0.016) x 10 6 s" 1 ( |Wiaux et al, 2002| ). 

One nice feature in the experimental study of the hy- 
perfine effect in the n B(3/2~, 0)-> n Be(l/2-, 320) tran- 
sition is the level structure of the 11 Be nucleus. Since 
the only particle-stable states are the l/2 + ground state 
and the l/2~ excited state this minimizes any concerns of 



experiments clearly demonstrate the strong dependence 
of the capture process on the hyperfine state. 

In allowed capture on 23 Na nuclei a large fraction of 
Gamow-Teller st rength is to several levels in the energy 
region 1-4 MeV flSiebels et al, 1995| ). In particular two 
A,r = -1+ reactions, ^Na(3/2+, 0H 23 Ne(l/2 + , 1017) 
and 23 Na(3/2+,0)^ 23 Ne(l/2+,3458), exhaust a sizable 
fraction of GT strength. For reference, the 23 Na disap- 
pearance ra te is A^, = (0.831 ± 0.002) x 10 6 s" 1 ( [Suzuki 



cas qadc feeding from muon capture to higher lying Be 



states. However one difficulty is the very low rate for 
the 11 B(3/2-,0)^ 11 Be(l/2-,320) transition, with only 
about 0.2% of n stops in n B undergoing this reaction. 
Two 11 B(3/2~,0)-> 11 Be(l/2~,320) mea surements, 



et al, 1987 ) and the 23 Na hyperfine conversion rate is 
A^ = (15.5 ± 1.1) x 10 6 s" 1 ( porringc et al, 1994| ) 



the pioneering work by Dcutsch et al (1968 ) at CERN 
and the recent work by Wiaux et al (2002) at PSI, have 



been performed. Both counted incoming muons in a plas- 
tic scintillator beam telescope, detected outgoing 7-rays 
in a high resolution Ge detector, and used a natural iso- 
topic abundance boron target. The 320 keV 7-ray signal - 
to-noise ratio was 1:2 in the earlier Deutsch et al (1966 ) 
experiment and 5:1 in the later Wiaux et al (2002[ ) ex- 



Compared to n B, in 23 Na (i) the capture rates are 
considerably larger and (ii) the hyperfine rate and disap- 
pearance rate are more readily distinguished, thus mak- 
ing the measurement more straightforward. However the 
large number of 23 Ne states, and greater fragmentation 
of Gamow-Teller strength, means cascade feeding from 
higher-lying levels to lower-lying levels is a worry. 

The study of the hyperfine ef fect in the /i 23 Na atoms 
was conducted at T RIUMF ( porringc et al, 1991 ; 
Johnson et al, 1996|). Incident muons were counted 



in a plastic scintillator telescope and stopped in a 
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sodium metal target. Emerging 7-rays were detected in 
high purity Ge detectors with surrounding Nal Comp- 
ton suppressors. For 23 Na(3/2+, 0)^ 23 Ne(l/2+, 1017) 
and 23 Na(3/2+,0)^ 23 Ne(l/2+,3458) the measurement 
yielded A+/A_ = 0.18 ± 0.03 and A+/A_ < 0.19 and 
revealed a very large hyperfine effect. Note that in ana- 
lyzi ng their data the authors had to account for the direct 



production and the indirect production of the 1017 kcV 



7~r 4ys. For more details on the interpretation of their 
data see pohnson et al. (1996 ) . 



nately measurement-to-measurement discrepancies of a 
factor of two are apparent for 16 O(0+,0) -> 16 N(0~, 120). 
Probably one should reject the early measurements with 
Nal detectors due to poor resolution and thus uncertain 
backgrounds. Further the peculiar time dependence of 
the interesting 120 keV 7-ray is probably the origin of 
disagr eement between the e arly studies of Dcutsch et al. 
(1969) and later studies of Kane et al. (1973 ) and |Gm" 



chon et al. (1979). In the later exp eriments special eff orts 
was made to tackle this difficulty. Kane et al. (1973) em- 



3. 16 O(0+,0) -» 16 N((T,i20) 

The 16 O(0+,0) -^ 16 N(0~, 120) transition is a first for- 
bidden transition from the J? — + , 16 ground state to 



, 16 N metastable state. As discussed in Sec. 
IX.D.3I the capture rates of AJ 71 " = 0~ transitions are 



the JjT 



ploy ed a continuous beam without a //-stop time gate 
and Guichon et al. (1979) employed a pulsed beam with 
a long /i-stop time gate. These experiments, using dif- 
ferent timing techniques and yielding consistent results, 
give a mean value for the 16 O(0 + ,0) — > 16 N(0~,0) rate 
of A = 1520 ± 100s- 1 . 



TABLE VIII Muon capture rates in 16 to the (0 _ , 120) and 



axial current, and hence to a„. The 16 N level structure Cohen et al. (1965). 


Astbury et al. (1961 


), 


Dcutsch et al. 


com 


prises four particle-stable bound states: the (2 , 0) (19690, 


Kane et al. (1973) and 


Guichon et al. (1979) in units 



excited states. Note that the 16 N(0", 120) state decays, 16 N (3 ,298) excited state. 



with a half-life ti/2 — 5.25/is, both by 7-ray emission to 
the 16 N(2~, 0) ground state and by /3-ray emission to the 
16 O(0 + , 0) ground state. Also note that the dominant de- 
cay of the 1~, 16 N state is via 7-ray cascade through the 
0", 16 N state. 

One experimental difficulty is the low rate of the 
16 O(0+,0) -> 16 N(0~,120) transition. This difficulty 
is compounded by the low energy of the de-excitation 
7-ray, resulting in large backgrounds from Michel 
bremsstrahlung and Compton scattering. Furthermore 
because of the meta-stability of the 16 N(0~,120) state 
the 120 keV 7-ray time spectrum is a convolution of the 
1.8/xs /i 16 lifetime and the 5.2/xs 16 N(0~,120) lifetime. 
Consequently some care is needed in applying /i-stop 
timing gates as the 120 keV 16 N 7-ray and other 16 N 
7-rays are impacted differently. Last the feeding of the 
16 N(0", 120) level from the 16 N(1~,398) level also com- 
plicates the extraction of A. 

The first experimental studi es of 16 — > 16 N p a rtial 
transitions were conducted by Cohen et al. (1963| ), Co- 



Ref. 



Cohen et al. (1963[ ) 



Astbury et al. (1964) 



Deutsch et al. (1969[ ) 



Kane et al. (197;) 



Guichon et al. (1979) 



1.1 ±0.2 
1.6 ±0.2 



n ss+ - 15 



1.56 ±0.03 
1.50 ±0.11 



1.7 ±0.1 
1.4 ±0.2 



■1 oc+0.36 
1.80_ 17 



1.31 ±0.11 
1.27 ±0.09 



hen et al. (1964) and Astbury et al. (1964) in the early 
sixties. The experiments detected the dc excitation 7- 
rays from 16 N excited states using Nal detectors. Un- 
fortunately, because of the limited resolution of the Nal 
detectors, they suffered from poor signal-to-noise and 
unidentified background lines. Consequently from the 
late sixties to the late seventies a further series of 7- 
ray experiments wer e conducted using Ge d etectors by 
Deutsch et al. (1969| ), Kane et al. (1973 ) and Guichon ei 



al. ( |1979| ). 



The experimental results for /i capt ure r ates to 16 N 
bound states are summarized in Table VIII . All exper- 
iments indicate substantial capture to the 2~ ground 
statfc, some capture to the 0~ 
and negligible capture to the 3~ 



4. 28 Si(0+,0) -> 28 AI(1+,2201) 

The study of g p by measurement of 7-ray angular cor- 
relations in \x~ [A,Z] — * v [A,Z-1]* — > 7 [A,Z-1]** transi- 
tions was originally proposed by Popov and co-workers. 20 
The authors examined 7-ray correlations in allowed and 
forbidden transitions on J = and J =/= targets, and 
emphasized the g p sensitivity of the + — > 1 + — * + spin 
sequence. 

An experimental method for 7-ray correlation mea- 
surements was subsequently suggested by Grenacs et al. 



and 1 excited states, 
excited state. Unfortu- 



iBukat and Podov C 19641): 


Bukhvostov and Podov (1967a 


C|); fDzicwicz and Popov 


(1965)); |Popov (1963 


)■ 



3G 



(1968). The method exploits the Doppler energy shift 
of nuclear 7-ray from in-flight decay. To illustrate the 
method we consider a recoil nucleus with a 7-decay life- 
time denoted by r and a stopping time denoted by ts (the 
stopping time is typically ~ 0.5 ps in medium- weight nu- 
clei). If r << ts the recoil is in motion as it decays and 
consequently the 7-ray energy is Doppler shifted by 



AE E~E Q 
—— = — — — = pcoso 
E n E n 



(46) 



where E a is the 7-ray energy in the recoil reference frame, 
E is the 7-ray energy in the laboratory reference frame, 
(3 = v/c is the velocity of the recoil in the laboratory, and 



is the angle between the 7-ray and the recoil momen- 



tum' As a consequence of Eq. (|46|) the energy spectra 
of nuclear 7-rays from muon capture are Doppler broad- 
ened when r << Tg. Further the exact lineshape of the 
Doppler spectrum is a function of the correlations be- 
tween th e 7-ray, recoil and /i -spin directions. As recog- 
nized by Grenacs et al. (1968 ), this permits a determina- 



tion of the 7-ray correlation coefficients by measurement 
of the 7-ray Doppler energy spectrum. 

Two experimental configurations for correlation mea- 
surements have special significance. In the first arrange- 
ment, the j—v configuration, the experiment is conducted 
with either unpolarized muons or perpendicular geome- 
try, so that P^-k — where k is the 7-ray direction. From 
Eq. ( 39 ) this configuration yields sensitivity to the corre- 
lation coefficient a 2 only. In the second arrangement, the 
"f-v-fj, configuration, the experiment is conducted with 
both polarized muons and non-perpendicular geometry, 



so that P^-k 7^ 0. From Eq. (39) this configuration yields 
sensitivity to the correlation coefficients a + |ci and 62 

also. Note that, due to different powers of k in Eq. (39), 



the Doppler lineshape arising from a 2 is symmetric about 
E whereas the Doppler lineshapes arising from a + |ci 
and 62 are asymmetric about E a . 

Measurements of 7-1/ and j—v—fj, correlations have pros 
and cons. One disadvantage of 7— 1/— )U correlation mea- 
surements is the small fi polarization in the //-atom 
ground state. However one advantage is that the Doppler 
lineshape may be manipulated by varying the /^-spin di- 
rection or 7-detector position, which is helpful in separat- 
ing the Doppler signal from the continuum background. 
Furthermore distortion of the 7-ray lineshape due to 
slowing-down of the recoil nucleus is more straightfor- 
wardly separated from the asymmetric effects of j-v—n 
correlations than from the symmetric effects of 7-^ cor- 
relations. 

Two n~ [A,Z] -> v [A,Z-1]* -> 7 [A,Z-1]** transitions 
have attracted the most attention, 



2^ 



Si(0+,0) 



2."v 



Al(l + ,2201) 



2x, 



Al(0 + ,973) (47) 



that the + — » 1 + — > + sequence involves a pure Ml 7- 
decay whereas the + — * 1 + — * 2 + sequence involves a 
mixed E2/M 1 7-decay with mixing ratio S(E2/M1) = 
0.37 ±0.11 flKudoyarov et ai, 19981) . Additionally the 
2201 keV state lifetime is 59 ± 6 fs and slowing-down ef- 
fects are non-negligible. Although at first glance the two 
sequences involve six correlations, i.e. (a+|ci) 
and b\ 229 for the 1229 keV 7-ray and (a+§ci) 
and b 2 170 for the 2170 keV 7-ray, these coefficients are 
related to a single underlying dynamical param eter, th e 



\1229 „1229 

; "2 
2170 „2170 



helicity amplitude ratio X, as described in Sec. IX.C.l 



The ground-breaking work on the 7-ray correlations in 
these spin se quences was conducted at SREL by Miller 



et al. (1972b). They measured the 7 ray lineshapes from 
fi stops in both natural Si and enriched 28 SiC*2 targets. 
Unfortunately the statistics were limited and the authors 
were forced to assume the absence of distortions of the 
Doppler lineshape due to slowing-down of the recoil ion. 
More recently an improved measurement of the 02 co- 



efficient was conducted at TRIUMF by Moftah et al. 



(1997). The experiment was performed in perpendicu- 



lar geometry, i.e. with ■ k = 0, and utilized a coin- 
cidence technique with Compton suppression to improve 
the 1229 keV 7-ray signal-to-noise. In analyzing the 
data the authors treated the 02 coefficients for 1229 keV 
and 2170 keV 7-rays as independent, since the 2170 keV 
multipolarity was not measured at the time. However 
they included recoil slowing-down effects in the fit of the 
lineshapes. 

Also recently a new measurement of all the c oefficients 
has been performed at the Dubna phasotron by Brudanin 



et al. (1995| ) and [Briancon et al. (2000| ). To identify both 



(a + Ici) and 62 the authors recorded the Doppler spec- 



tra for different values of P^ • k. In Brudanin et al. (199£ ) 
they used forward/backward positioned Ge detectors to 
vary P^ ■ k and in |Briancon et al. (2000| ) they used /i 
spin precession to vary P M • k. In analyzing their data the 
authors enforced the dynamical relations between cor- 
relation coefficients and fit the 2201 keV lifetime and 
2170 keV mixing ratio. 

Note that a concern in all experiments is the produc- 
tion of the 28 A1(1+,2201) state by either (i) (fj,, nv) or 
(/i, nnz/) capture on 29 Si or 30 Si isotopes, or (ii) 
capture to higher-lying 28 Si levels. Such contributions 
would distort the lineshapes an d impact the extrac tion 
of (a+§ci), a 2 and 6 2 . Note that [killer et al. (1972a| ) ob- 
tained constraints o n co ntributions from 29 Si(/ i, nu) and 
Moftah et al. (1997[ ) and |Briangon et al. (2000[ ) obtained 
limits on cascade feeding from higher levels. However a 
small contribution from (i) or (ii) is not excluded. 

In Table IX we summarize the various measurements of 
7-ray correlations for 28 Si(0 + ,0) 



28 A1(1+,2201). The 
experimental configurations are denoted by P^ ■ k = 



2>N 



Si(0+ 0) 



28 A1(1 + , 2201) ► 28 Al(2 + ,37). (48) or P M • fc ^ and indicate the absence or presence of 



They involve a common GT transition from the 28 Si(0 + , 
0) ground state to the 28 A1(1+, 2201 ) excited state. Note 



sensitivity to (a+ |ci) and b 2 . To assist the comparison 
of experiments we quote all results in terms of the helicity 
amplitude ratio X. The 'world data' weighted mean is 



37 



X = 0.554 ±0.042. 



TABLE IX Compilation of recent results from the j— 
ray correlation experiments for the transition 28 Si(0 + ,0) — > 
28 A1(1 + , 2201). The results are presented in terms of the dy- 
namical parameter X = v2Trj /T\ . The analyzed gamma-rays 
are listed in column two and the experimental technique is 
listed in column three. FPA is the Fujii-Primakoff approxi- 
mation value. 



Ref. 


7-ray 


7-ray X = v^To/Ti 




trans. 


corr. 



Moftah et al. (1997| ) 1229 



V' 0.454l°;i 2 



Brudanin et al. (1995| ) 1229, 2171 'j-v-fi' 0.543 ± 0.052 
Briangon et al. (2000| ) 1229, 2171 'j-v-fj,' 0.566 ± 0.045 
FPA 0.59 



F. Theoretical framework for exclusive OMC 

Herein we describe the theoretical treatment of physi- 
cal observables in (//, v) reactions. Our main goals are to 
outline the steps and assumptions in calculating the ob- 
servables and provide some details o f work on 11 B, 12 C, 



16 0, 23 Na and 28 Si. In Sec. |IX.F.1| we describethe op- 



erators that contribute to \i capture. In Sec. IX. F. 2 we 



discuss the application of the impulse approximation and 
in Sec. IX. F. 3 we discuss the effects due to exchange cur- 
rents . The detailed discussion of nuclear models is left to 
Sec. tX.G. 



Cj — Mj involves a longitudinal coupling to the lepton 
field and is dependent on g p , while q~ji—q~p a 9 i nvo l V es a 
transverse coupling to the lepton field and is independent 
of g p . 

+ — > 0~ transitions, e.g. 16 0, involve a unique J 71 = 
0~ multipole, a single electroweak weak operator, Lq — 
Mq , and two multipole operators, Mqi ■ o" and M cr ■ V. 
The Moi ■ cr operator is the I = 1 retarded Gamow- 
Teller operator originating from the space component of 
the axial current. The M^cr ■ V operator is the axial 
charge operator originating from the time component of 
the axial current. Note that the contribution of g v in 
+ — > 0~ transitions is via Mqi ■ <J - 

+ — * 1 + transitions, e.g. 12 C and 28 Si, involve a 
unique J n = 1 + multipole, two electroweak weak op- 
erators, L\ - M\ and Tf 5 ~ T" lag , and four multipole 
operators, Mio ■ cr, M\2 ■ cr, M\cr ■ V and Afn • V. 
The M io • cr operator is the allowed GT operator. The 
remaining contributions include the axial current's time- 
component, i.e. Mia ■ V, and second-forbidden correc- 
tions, i.e. M12 ■ cr. Note that the leading contribution 
of g p in + — * 1 + transitions is via Miq ■ cr. 

For transitions on Jj =/= targets a range of multi- 
poles are involved, i.e. \Ji — Jf \ to (J^ + Jf). For exam- 
ple, a l/2 + — > l/2 + transition, e.g. 1 H or 3 He, involves 
,r = 0+,l+ multipoles and a 3/2+ -> 1/2+ transition, 
e.g. 1X B or 23 Na, involves J T = l + ,2 + multipoles. For 
1/2+ — > 1/2+ transitions the Lq — Mo operator yields 
an additional contribution from the allowed Fermi op- 
erator. For 3/2+ — > 1/2+ transitions the — M2 and 
^ei rpmagb operators yield additional contributions from 
the ^-forbidden multipole operators M2 and M22 • cr and 
the gradient multipole operators M 2 i • V and M23 ■ V. 
However the leading contribution of g p in l/2 + — > 1/2+ 
transitions and 3/2 + — > l/2 + transitions is still via 
Mio ■ cr. Note that the multipoles with J x = + and 
= 2 + are independent of g p . 



1. Multipole operators 

The model calculation of partial transitions on com- 
plex nuclei is generally conducted via a multipole expan- 
sion of the Fujii-Primakoff effective Hamiltonian of Eq. 
(Eol) . We refer the reader who is interested in the details 



of the formalism to the articles by Primakoff (1959), 


Morita and Fujii (196C), 


Luyten et al. (196S), Walecka 


(197 
Ozi( 


5), Mukhopadhyay (1977 


) and Ciechanowicz and 


wicz (1984]). Herein we briefly describe the specific 



operators that contribute to transitions of interest for 
g p , i.e. the transitions + — > 0~, + — > 1 + and 3/2 ± — > 
l/2 ± . We employ the notation of Walecka (1975[ ) for 
the contributing electrowe ak operators, (Cj — A4j and 



Tf l —T™ a9 ), and notation of Donnelly and Haxton (1979 



1980] ) for the contributing multipole operators, Mj(qx) 
M JL (qx)-cr, M JL (qx)-V and Mj(qx) cr V . Note that 



2. Impulse approximation 

In principle the weak amplitudes in muon capture have 
one-, two- and many-body contributions. However in 
practice the starting point for most calculations is to ap- 
proximate the weak nuclear amplitude by a summation of 
A one-body amplitudes, i.e. the impulse approximation. 
This amounts to ignoring the effects of pion exchange 
currents, A-hole excitations, etc. 

Assuming a one-body form for nuclear currents, 
the required multi-particle weak matrix elements 
< J/||0 ,7 ||Ji > between an initial state |Jj > and final 
state I Jf > may be written in terms of single-particle 
weak matrix elements < a'HO^Ha > between single- 



particle states labeled by \a > and \a' > as (Donnelly 



38 



and Haxton , 1979) 
< J f \\0 J \\J t >=J2C(J,a,a',Jf,Ji) < a'\\0 J \\a > 

a, a' 

m- 



allowed Fermi operator this approach is not helpful. In- 
stead the modifying effects of exchange currents must be 
addressed by explicitly evaluating a specific set of Feyn- 
man graphs. An example of such an approach is the 
N7rpAi phcnomcnol ogical Lagrangian model of Ivanov 



and Truhlik (1979a ) . Note that the Feynman graphs and 



where the states a = [n,j,£] and a' = [n',j' 
capture the operator O 7 



is M 



J 



C 



j or T] 



and in fi 



The coefficients C(J,a,a' , J/, Jj) are called one-body 
transition densities and determine the contributions of 
each single-particle matrix element < a'||0 J ||a! > to 
the multi-particle matrix element < J/||0 J ||Ji >. Note 
that the one-body transition densities are determined by 
the nuclear structure while the single-particle matrix el- 
ements are functions of the weak couplings. Eq. flls| ) 
therefore represents a convenient separation of nuclear 
structure and weak dynamics. 21 We stress that in prac- 
tice in applying Eq. (|49j) the summation is truncated to 
a finite number of the single particle transitions. 



3. Exchange currents 

At some level the impulse approximation will break 
down, and consequently the evaluation of contributions 
from pion exchange, delta excitation, etc., is important in 
extracting g p . One approach to computing exchange cur- 
rents is to use low energy theorems to constrain soft pion 
contributions. Another approach involves enumerating a 
plausible set of Feynman diagrams that incorporate it's, 
A's, etc. We shall not attempt to cover in detail the 
broad topic of exchange currents, but rather we summa- 
rize their application to muon capture on complex nuclei. 

The work of Kubodera et al. (1978| ) was pivotal in es- 
tablishing the importance of soft pion exchange currents 
in various electroweak processes. The authors observed 
that soft pions produce large effects in the time compo- 
nent of the axial current and the space component of the 
vector current. Further they recognized that arguments 
based on chiral symmetry fix the size of these effects, 
yielding a powerful tool in determining the corrections 
from 2-body currents. In particular for muon capture 
the soft-pion corrections to axial charge operators, while 
substantial, are well determined. Such Mjcr ■ V opera- 
tors compete with the leading contributions of g p in both 
allowed and + — » CP transitions. For further detai ls see 
Guichon et al. (1978); |Guichon and Samour (1979 ) and 



Towner (1986^ 



Unfortunately, for the space component of the axial 
current and the time component of the vector current the 
constraints dictated by chiral symmetry are ineffective in 
determining the contributions of 2-body currents. For 
example for the allowed Gamow- Teller operator and the 



corresponding operators are obviously not unique and the 
couplings constants and form factors are frequently not 
well known. Consequently the calculation of 2-body cor- 
rections to GT matrix elements and Fermi matrix ele- 
ments have significant uncertainties. In particular the 
A's contribution is poorly determined. For examples of 
applications to muon capture see |Adam et al. (1990| ). 



G. Nuclear models for partial transitions 

Next we consider the specific structure and nuclear 
models for partial transitions on A= 11,12,16,23 and 
28 nuclei. Our focus is on the elements of the models 
that impact the determination of g p . 

Nowadays full-space shell model calculations are rou- 
tinely performed for 4 < A < 16, i.e. Op shell nu- 
clei, and 16 < A < 40, i.e. ls-Od shell nuclei, and 
the empirical determination of effective interactions from 
least-squares-fits to nuclear data is well established. 
For example, the parameters of the Op shell interaction 
were obtained by |Cohcn and Kurath (196"?] ) by fitting 
4 < A < 16 level energies and the parameters of the 
ls-Od interaction were obtained by Wildenthal (1984) by 



fitting 16 < A < 40 level energies. In addition semi- 
empirical interactions, which incorporate assumptions for 
the p articular form of the effective interaction, are avail- 
able ( Brown et at, 1988|) . Both empirical and semi- 



empirical interactions are capable of reproducing many 
features and phenomena in Op and ls-Od nuclei. 

Note that an alternative approach is the microscopic 
derivation of the effective NN interaction from the free 
NN interaction. The derivation involves a power-series 
relating the effective interaction and free interaction, 
but unfortunately the question of convergence is tricky. 
In general such interactions either give less satisfactory 
model-data agreement than empirical interactions, e.g. 
the interaction of Kuo and Brown (196(f ), or need some 
ad-hoc tuning of model parameters, e.g. the interaction 



of Haugc and Maripuu (1973). However the compar- 



ison of results from empirical interactions and realistic 
interactions is helpful in understanding and evaluating 
the model uncertainties. 



1. n B(3/2-,0) 



L Be(l/2-,320) 



St rictly the radial form of the nuclear wave functions also enters 



Several calculations of capture rate and hyperfine de- 
pendences have been perfor med for the n B(3/2~,0) — > 
n Be(l/2~ 3 20) t r ansition (|Bcrnabcu 1971|; [Koshigiri 
et al., jjjgg, |1984| ; [Kuz'min et al, 1994| ftuzuki , 1997Q 



the computation of the single— particle matrix elements. 



Bernabeu (1971) first discussed the relatively high sensi- 
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tivity to g p and relatively low sens itivity to nuclear struc - 



ture of A + /A_. More recently Kuz'min et al. (1994 ) 



assessed th e effects of diffe rent Op-shell effective interac- 
tions while [Buzuki (1997 ) assessed the effects of a 11 Be 
neutron halo. 

The simplest picture of n B(3/2~,0) -> 
n Be(l/2~, 320) consists of a 0s*i 2 QpJ/ 2 initial state, 

0s^/ 2 0pg i 2 0p\ i 2 final state, and 0p 3 /2 — ► Opi/2 single 
particle transition. However in reality the 11 B(3/2~,0) 
initial state has substantial contributions from both 



o4/2 



°P3/2 



0Pi/2 an d ® s i/2®P3/2®Pi/2 configurations, 
and interference between the 0pi/ 2 — > QP3/2 single 
particle transition and the 0p 3 /2 ~ * Ofi/2 single particle 
transition is important. For example the capture rate 
is grossly over-estimated in a simple 0p 3 / 2 — > 0pi/ 2 
picture. 

Full Op-shell model calculations for n B(3/2~,0) -> 
n Be(l/2~, 320) with well established effective interac- 
tions confirm the leading transition is 0p 3 /2 — > Qpi/2j 
with a substantial contribution from 0pi/ 2 —> 0p 3 / 2 and 
a significant contribution from 0p 3 /2 — > 0p 3 /2- Since the 
Gamow-Teller matrix elements for 0p 3 / 2 - ► Qpi/2 an d 
QPi/2 ~~ * 0j»3/2 have opposite signs they interfere destruc- 

and A_. Typ- 



The simplest picture of 12 C(0+,0) -> 12 B(l+,0) con- 
sists of a 0sy 2 0p^, 2 initial state, 0s^, 2 0j>g , 2 0pL 2 final 
state, and a 0p 3 / 2 - ► 0piy 2 single particle transition. 
Like the 11 B g.s., the 12 C g.s. has substantial contri- 
butions from both 0sy 2 0py 2 0p^ 2 and 0sf^ 2 0py 2 0p^ 2 
configurations, and the interference of amplitudes from 
OP3/2 ~~ * 0p 3 / 2 transitions and Opx/2 ~~ * OP3/2 transitons 
is important. Again the capture rate is grossly over- 
estimated by a simple 0p 3 / 2 — + 0pi/ 2 picture. 

Note that the overall pattern of the transition densi- 
ties obtained from full Op-shell model calculations with 
well established effective interactions is quite similar 
for n B(3/2-,0) -> n Be(l/2-,320) and 12 C(0+,0) -> 
12 B(l + ,0). Specifically for 12 C the largest density is 
QP3/2 ~^ 0pi/2, the next-to-largest density is 0pi/2 — * 
0P3/2> and the contributions from 0p 3 / 2 — > 0p 3 / 2 and 



Op 



1/2 



tively and dramatically decrease both A + 



ical model-to-model variations in the 0p 3 / 2 — > 0pi/ 2 and 
0pi/2 - * OP3/2 densities are roughly ±10%. Note that 
the remaining 0pi/ 2 — > 0pi/ 2 and 0p 3 / 2 — > 0p 3 / 2 den- 
sities show larger model-to-model variations, however 
their contributions to capture are smaller. 

A unique complication for the n B(3/2~,0) — > 
n Be(l/2-, 320) transition is the n Be neutron halo. The 
halo is interesting in the context of nuclear structure 
studies b ut worrisome in the context of extracting g p . 
Recently ISuznki (19971) has assessed the impact, of the 



0pi/2 are small. As in 11 B, in 12 C the interfer- 
ence of 0p 3 / 2 — > 0pi/2 with 0pi/2 — > 0p 3 / 2 is important 
in reducing the Gamow-Teller matrix element. 

A nice feature of 12 C(0+,0) -> 12 B(l+,0) is that 
related data on other electroweak processes are avail- 
able, e.g. 12 B /3-decay, 12 C(e,e') scattering and 
12 C(1 + ,15.1 MeV) 7-decay, and are helpful in testing 
the model calculations. For example the one-body tran- 
sition densities have been extracte d fro m these data by 
Doyle and Mukhopadhyay (19951 ) an d [Haxton (197^ ) 
and generally support the model calculations. 



3. 23 Na(3/2+,0) ^ 23 Ne(l/2+, 3458) 

Calculations of capture rates and hyperfine dependen- 



cies m 



5 Na 



2.3 



Ne t ransitions have been p erfor med by 



Johnson et al. (1996| ), Koshigiri et al. (1997] ) and Biisko 



iieuJron halo on the capture rate and its hvpcrfinc dcpcn- ncn et al. (1998j)TThe calculations have been conducted 



dence. Suzuki (1997) reported it tends to reduce both 



A + and A_ but barely changes A + /A_ 



2. 12 C(0+,0) 



2 B(l+,0) 



A large number of model calculation are available for 
12 C(0+,0) -> 12 B(l+,0). Early studie s of (j, capture rates 



and /3- decay rates were performed by Flamard and Ford 
(19 ^91), lO^onncll et al. (1972[), Mukhopadhyay and Mag 



farhne (1971 jTlfrnmelc and Mukhopadhyay (1975 ), and 



Mukhopadhyay and Martorell (1978). More comprehen- 
sive studies of rates and polarizations in fi capture were 
pub lished by Morit a and co-workers, 22 a nd by Subra- 
marian et al. (1976), Roscnfcldcr (1979| ), Ciechanowicz 
( 1 9 S 1|) , Hayes and Towner (2000] ) and Auerbach and 
Bro vn (2002). They include investigations of two-body 



currents and core polarization effects. 



in the full ls-Od model space using the Wildcnthal empir- 
ical interaction ( Wildcnthal, 198 4| ) and the K uo-Brown 
realistic interaction ( JKuo and Brown , 1966 ). A study 
of 2-body currents a nd core polarization was made by 
Koshigiri et al. (1997] ). 

The 23 Na ground state is J 71 " = 3/2+ , indicating the 
simplest picture of a single unpaired proton in a 0g?5/2 
orbital is wrong. Rather the A = 20-24 mass region 
is well known for examples of light deformed nuclei and 
rotational band spectra. Consequently the spherical shell 
model representation of these nuclei is quite complex, 
with A = 23 wave functions having small occupancies for 
the (d 5 / 2 ) 7 configuration and large occupancies of the 
lSx/2, ld 3 / 2 orbitals. 

For concreteness we describe the 23 Na(3/2 + ,0) — > 
23 Ne(l/2 + , 3458) transition, which is the strongest tran- 
sition in the /x _23 Na experiment. Full ls-Od shell model 
calculations with well tested effective interactions show 
0rf 3 / 2 is the strongest single particle transition 



Od 



5/2 



22 See 


Ami et al. C 1981 


I: Pukui et al. (1983a i. 


1987): 


Kobavashi et 


al. ( 


1978); Koshigiri et al. (1985); Morita et al. (1994J). 



and 0d 3 /2 — > 0d 3 / 2 is the next strongest single particle 
transition. Other contributions are typically 10-20% of 
0^5/2 — * 0d 3 / 2 - The variations of the one-body transition 



densities between models are < 10% for Od 



5/2 



Od 



3/2 



40 



densities and < 25% for — > 0<i3/2 densities. 



4. 2S Si(0+,0) -> 28 AI(1+,2201) 

Several authors have calculated the rates and correla- 
tions for 28 Si( 0+,0) -> 28 A1(1+,220 1). The first efforts 



wer e made by Ciechanowicz (1976 ) and Parthasarathy 
and Sridhar (1978| , 1981| ), but employed relatively crude 



nuc lear wave functions, t he lp-lh wave function of Don- 
nell y and Walker (1970) and the truncated ls-Od wave 



functio n of |McGrory and Wildenthal (197l|). More re- 



cently Kuz'min and Tetereva (2000[). [Kuzmin et al 




siiskonen et al. (1999), and Ciechanowicz et al. 



admixtures. They found the effects of 2p-2h configura- 
tions on the axial charge matrix element Mqct ■ V and the 
retarded Gamow- Teller matrix element M\q ■ a were op- 
posite in sign, and consequently the near-cancellation of 
model uncertainties in A u / A g broke down. Subsequently 
Towner and Khanna (1981 ) evaluated the effects of 2p- 
2h configurations in ie N(0 , 120) <-> 16 O(0+, 0) transi- 
tions with more extensive configurations and various ef- 
fective interactions. They concluded that 2p-2h config- 
urations decrease the /3-decay rate by factors of 2 to 4, 
decrease the \i capture rate by factors of 1.5 to 2, and 
increase A^/Ap by factors of 1.5 to 2, thus confirming 
the m odel dependence observed by puichon and Samour 
Xl979| ). 



have conducted full ls-Od calculations |5iiskoncn| Recently jHaxton and Johnson (199C| ) and [Warburton 



1999|) and |Uicch anowicz et al. (1998|) also con 
sidered both 2-body currents and core polarization in 
28 Si(0+ 0) 28 A1(1+,2201). 
In the simplest picture 28 Si(0 



,0) 



consists of a full 0dg 2 2 initial state, J 71 " = 1 + (Od 
0d^ 2 ) final state, and 0d 5 / 



28 A1(1+,2201) 

5/2' 

2 — > 0d 3 / 2 single particle tran- 
sition. However the model calculations show the simple 
picture is insufficient and configurations with several par- 
ticles in lsi/2-0d3/2 orbitals are important. 

A special remark on the one-body transition densities 
in the 28 Si(0+,0) -> 28 A1(1+, 2201) transition is worth- 
while. Unlike the previous examples of n B, 12 C and 
23 Na, in 28 Si the variations f rom model to model are 
large, e.g. the densities from Kuo and Brown (196€ ) 



and Wildenthal (1984) are quite different. Also no par- 
ticu lar single particle transition is dominant, e.g. the 



Or/ 



inte 'action of Wildenthal (1984) shows 0^5/2 



5/2 



Od 



Od 



3/2; 



5/2j 



and lsj/2 — > Od^ /o with similar magni- 



tudes. Therefore, as discussed by Kuz'min and Tetereva 
(20Cj^), the model calculations are especially sensitive to 
interference effects. 



et al. (1994 ) have performed full 4 Klo (3 Hlo) calculation 
for 16 O ( le N) low-lying levels. These calculations nicely 
reproduce the excitation energies of the 16 O isoscalar pos- 
itive parity states and the 16 N isovector negative parity 
states, and indicate significant 4p-4h probabilities in 16 0. 
In both works the authors stress the large destructive in- 
terference between MqCT ■ V and Mio ■ er in /3-decay, 
and therefore substantial model dependences in Ap and 
A^/Ap. Consequently the authors argue that A p , not 
A^/ ' Ag, is preferable for extracting the coupling g p . 

The foregoing discussions of 16 O(0+,0) «-> 16 N(0~, 120) 
transitions show a worrisome sensitivity to the multi- 
particle wave functions. In addition the matrix element 
Mq(t ■ V is modified considerably by 2-body currents 
from soft -pion exchange. The se eff ects are discussed in 
detail by Guichon et al. (1978| ) and Towner and Khanna 



(1981). The calculations indicate that they increase the 
/3-decay rate by a factor of about four and increase the /i 
capture rate by a factor of about two. The contribution 
of MqCT ■ V is larger in /3-decay than /x capture. This 
further complicates the extraction of g p from capture on 
16 0. 



5. ie O(0+,0) -^ 16 N((T,120) 



Because the nucleus 16 O is doubly magic and the tran- 
sition 16 O(0+, 0) -> 16 N(0~, 0) is first forbidden this case 
is special. The simple model for 16 O(0+, 0) -> 16 N(0~, 0) 
comprises a Op closed shell initial state, and a lp-lh 
J n = 0~ final state, and involves lpi/2 — > 2s!/ 2 and IP3/2 
— > lc?3/2 single particle transitions. Several authors have 
computed the rates of fx capture (A M ) and /3-decay (A/3) 
for 16 O(0+ 0) <-> 16 N(0~,0) wit hin this scheme. For ex- 
ample see |Guichon et al. (1978 ) and references therein. 
They found the rates to be strongly dependent on the 
0p-i/9~ 1 -0d-i/9 admixture in the 16 N wave function, but 



H. Coupling g p from partial transitions 

I. Recommended values of g p /g a 

In Table [x] we list our recommended values for g p /g a 
from exclusive OMC on complex nuclei. By 'recom- 
mended values for g p /ga we mean our assessment of the 
best values from the current world experimental data- 
set and the most complete model calculations. In quot- 



the ratio A u / Ag to be near-model independent. 


Possoz et al. (1974), possoz et al. (1977 


)• 


Roesch et al. 


Si 


nee the work of 


Brown and Green (1966) we know 


(1981b|) and Kuno et al. (19840 to obtain 


a world aver- 



the above picture is not complete, and that 2p-2h con- 
figurations in the 16 O ground state are important. The 
effects of 2p-2h configurat ions on /i capture and /3-decay 
were first investigated by puichon and Samour (1979|) 
who considered (2s 1/2 ) 2 (lpi/2)~ 2 and (ld 3 / 2 ) 2 (lp3/2) 



age of the helicity amplitude ratio X = 0.23 ± 0.06 for 
12 C(0 + ,0) ^ 12 B(l + ,0) and comb i ned the experimenta l 
resul ts of iBrudanin et al. (1995j ), |Moftah et al. (1997| ) 
and Briancon et al. (2000Q toobtain a world average 
of the helicity amplitude ratio X — 0.55 ± 0.04 for 
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28 Si(0+,0)^ 28 Al(l+,2201). 23 For t he hyperfinc deperi - 
dence on n B we employ the results of Wiaux et al. (2002 ) 



and for th e hyperfine dependen ce on 23 Na we employ the 
results of Johnson et al. (1996 ). For the capture rate of 
the 16 O(0+,0)^ 16 N(0 ~,120Uransition we a veraged the 
exp erime ntal results of Kane et al. (1973 ) and Guichon ei 
al. ( |L979D . Note that in Table |Xj the quoted errors include 
only experimental uncertainties. Also note we quote the 
results in Table [x] in terms of g p /g a not g p . In most 
cases the measured quantities are recoil polarizations, 7- 
ray correlations or hyperfine dependences, and therefore 
are governed by ratios of nuclear matrix elements and of 
weak coupling constants. Consequently quoting g p /g a is 
more natural and more appropriate. 

In order to extract the coupling g p /g a from experi- 
mental data a model is necessary. Our model choices, 
and arguments for selecting them, are given below. 



TABLE X Recommended values of g p /g a from ordinary cap- 
ture on complex nuclei. The quoted errors are experimental 
uncertainties, and do not include model uncertainties. 



transition 



u B(3/2-,0) -> n Be(l/2-,320) 
12 C(0+,0) -> 12 B(l+,0) 
16 O(0+,0) -» 16 N((T,12u) 
23 Na(3/2+,0) ^ 23 Ne(l/2+,3458) 



For n B we took the results of Kuz'min et al. (1994) 
yielding g p /g a 



2Sr 



Si(0+,0) ^ 2S A1(1+,2201) 



9.8 ±1.8 

6.0 ±0.4 
6.6± 2 1 



i-olri 



These authors used the full Op 
space with Cohen-Kurath interaction but omitted the 
effects of core polarizatio n and exchange c urrents. Note 
the earlier calculation of Bernabeu (19711) yields a sim - 

ilar value of g p /g a = 4±|. According to (Suzuki (1997| ) cu rrents from soft pions and A excitati ons. The results 
the effects of the n Rc neutron halo on the hyperfine de- of [Johnson et al. (1996|) are similar to [Koshigiri et 



pen dlence are small. 



(1997) 



F or 12 C we used the model calcu lations of Fukui ei 
al. (1987) and excited state yields of Giffon et al. (1981) 
to obtain g p /g a — 9.8 ± 1.8. These authors used the full 
Op space with Hauge-Maripuu interaction and accounted 
for com polajdzatkm and exchange cuxmnta The mrp 1 



' po- 



larisation effects were computed to second-order in per- 
turbation theory and the exchange currents were com- 
with enntributinns from pair currents, pionic cur= 



puted 



For 28 Si we took the results of Siiskoncn et al. (199S ) 
yielding g P /g a = l-0±i'2- The authors used the full ls- 
Od space with the Brown- Wildenthal interaction. They 
included core polarization corrections but omitted ex- 
change current corrections. Note that Ciechanowicz ei 
al. (1998| ) found the effects of soft -pion exchange were 
negligible. For comparison, the earlier calculations of 
Ciechanowicz (1976] ) and jParthasarathy and Sridhar 



rentk and A excitations Earlier calculations with more 

rudimentary wave functions and less sophisticated treat- 
ments of core polarization and exchange currents gave 
similar results. 

For 16 a number of determinations of g p /g a are pub- 



1981) gave values for g p /g a of 3.3 ± 1.0 and 1.5±" ^ re- 



spectively. 



2. Model sensitivities of g v j 



lish ed. We took th e value g p /g a — 6.0 ± 0.4 from War 
burton et al. (1994), which incorporates both 4p-4h ib O 
configurations and 3p-3h 16 N configurations and repro- 
duces the 16 N /3-decay rate. The 4p-4h calculation of 
Haxton and Johnson (199C ) yields a similar value of 
g p /ga = 5-7. We note, however, that the earlier calcu- 
lations which omit 4p-4h configurations have generally 
preferred higher values for g p /g a - 

For 23 Na we took the results of Koshigiri et al 



The interesting observables in allowed transitions are 
governed by ratios of nuclear matrix elements. More 
specifically, in + — > 1 + transitions the observables are 



completely determined by, and in 3/2- 



1/2* tran- 



(1997) 



yielding g p /g a = 6.6±|'|. The authors used the full ls- 
Od space with the Brown- Wildenthal interaction, and ac- 
counted for core polarization to first-order and exchange 



sitions the observables are strongly dependent on, the 
helicity amplitude ratio X = \/2To/Ti. 24 Therefore un- 
derstanding the model uncertainties in computing X is 
central to tracing the model dependences in recoil polar- 
izations, 7-ray correlations and hyp erfine dependences. 
For example see Junker et al. (2000| ). 

Tables jx| and XII show the relative contributions of 
M 10 • <t, etc., to If- Ml (or T ) and Tf 5 - M™ ag (or 
Xi) for the relevant transitions. The tables show that 
Tf' 5 — M" lag is entirely dominated by M\o ■ er. However, 



23 In the literature both the helicity amplitude ratio, denoted X, 
and the neutrino-wave amplitude ratio, denoted x, have been 
used in this context. Although both X and x are suitable for 
representing the dynamical content of AJ W = 1+ transitions, 
they are different, i.e. X = (-2x + V2)/(x + y/2). Th erefore it 
is important not to confuse the two variables. See Sec. [X.F.I 



24 Note for A.J 1 ' = 1 + multipoles the helicity amplitude ratio 
V2T /T! and multipole amplitude ratio V2{L\ - M^)/(Tf B - 



M™ aB ) are identical 
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while Mio ■ a is the leading piece in the L\ — Ml term, 
the contributions from Mia ■ V of about 30-40% and 
M12 ■ a of up to 25% are important. Clearly the leading 
source of model dependence in computing X is therefore 
uncertainties in the ratios ot M\<j ■ V/iW io • a and M 12 • 
a/M 10 • a, i.e. corrections arising trom axial charge and 
2" -tor bidden cttects. 



TABLE XI Comparison of the corrections from the terms in- 
volving M\o ■ V, M12 ■ cr, and M11 ■ V to the multipole 
amplitude of L\ — M\ for A = 11, 12, 2 3 and 28 and several 



TABLE XII Comparison of the corrections from the terms 
involving M\cr ■ V, M12 • cr, and M11 ■ V to the multipole 
amplitude of Tf 5 -T™ ag for A = 11, 12, 23 and 28 and several 



effective interactions denoted CKPOT ( Cohen and Kurath 



1967), PKUO (Kuo and Brown , 196 



USD ( [Wildenthal 
196 



1984| ) and KUOSD ( [Kuo and Brown , 19661 ). The values in 
columns 3-5 correspond to the percentage change in the ma- 
trix elements as the correction terms are successively included. 
Note that Micr • V doesn't contribute to Tf 5 - T 1 ma9 . 



effective interac tions denoted CKPOT (Cohen and Kurath 
1967), PKUO (|Kuo and Brown , 1966|), U SD ( jwiidcnthal 
1984) and KUOSD ( |Kuo and Brown , 1966| ). The values in 
columns 3-5 correspond to the percentage change in the ma- 
trix elements as the correction terms are successively included. 
Note that Mn • V doesn't contribute to L\ - Mj\ 



A Int. 


Lf-Mf Mi cr 


V 


M12 • cr 


Mn ■ V 




corr. 


(%) 


corr. ( % ) 


corr. ( % ) 



11 


CKPOT 


-0.103 


-30.6 


-2.6 


0.0 


11 


PKUO 


-0.094 


-32.4 


0.6 


0.0 


12 


CKPOT 


0.107 


-34.5 


-9.4 


0.0 


12 


PKUO 


-0.071 


-41.6 


-18.4 


0.0 


23 


USD 


0.090 


-31.8 


4.1 


0.0 


23 


KUOSD 


-0.071 


-44.4 


13.1 


0.0 


28 


USD 


0.070 


-37.2 


-27.5 


0.0 


28 


KUOSD 


-0.071 


-54.5 


-18.4 


0.0 



A 


Int. 


J i - J i 


Micr 
corr. 


V 

(%) 


M12 ■ cr 

corr. ( % ) 


Mn ■ V 

corr. ( % ) 


11 


CKPOT 


-0.243 


0.0 




0.9 


0.6 


11 


PKUO 


-0.222 


0.0 




-0.2 


1.7 


12 


CKPOT 


0.249 


0.0 




3.1 


-0.4 


12 


PKUO 


-0.164 


0.0 




5.4 


2.1 


23 


USD 


0.213 


0.0 




-1.4 


-2.2 


23 


KUOSD 


-0.170 


0.0 




-3.6 


-4.8 


28 


USD 


0.171 


0.0 




8.6 


4.7 


28 


KUOSD 


0.106 


0.0 




4.2 


0.7 



and M12 • cr/Mio 

have been studie d by Fukui et al. (1987) for 



ncn et al. 
"S^Na! 



Such core polarization effects 
12 C, giisko 



For 



Uncertainties in one-body transition densities are im- 
portant sources of model dependencies in Aficr • V/7Vfio • 
cr and -M12 • cr/M 10 ■ a. For example let's consider 
the transitions on boron and carbon which involve the 
interference of a 0p 3 / 2 - ► OP1/2 single particle transi- 
tion and a 0pi/ 2 — * 0p 3 / 2 single particle transition. Un- 
der the interchange of initial and final states the matrix 
element .Mio ■ cr changes sign but the matrix clement 



(1999b for *Si, and |Koshigiri et al. (1997| ) for 
12 C |Fukui et al. (1987jT found downward renor- 
malizations of 13% for Mio ■ cr, 33% for Mjj ; cr, and 
22% for M x a- -V, and for 28 Si |Siiskonen et al. (1999) ) 



Mi c ■ V docs not. Consequently Mia ■ V/M w ■ a is Pincdo and Vogcl (1998f T 



found downward renormalizations of 11% for Mio ■ cr, 
30% for M 12 • cr, and 39% for M x a ■ V. The small ef- 
fect of core polarization on M\§ ■ a is because the model 
spaces are 'complete spaces' for this operator. Note that 
the renormalization of Mi§ ■ a has strong support from 
experimental data on allowed /3-decay and (p,n)/(n,p) 
reactions, and the renormalization of Myi • a has some 
support from experime ntal data on secon d for bidden 0- 
decay. For example see Warburton (1992 ) and Martinez- 



quite sensitive to the Opi/2 —>■ QP3/2 admixture in these 
A = 11, 12 transitions. Fortunately the destructive inter- 
fere nce of Opi/2 ^ OP3/2 amplitudes is also reflected in, 



and |thus calibrated by, the capture rates of 11 B(3/2 , 0) 
-> 11 Be(l/2-,320) and 12 C(0+,0) -> i2 B(l+,0). 

A special comment is worthwhile for 28 Si(0 + ,0) — > 
28 A1(1 + , 2201). The calculations indicate this transition 
involves the destructive interference of numerous single 
particle t ransitio ns with comparable one-body densities. 
See Sec. IX. G. 4 for details. Consequently for 28 Si the 
calculation of X may be especially sensitive to the model 
uncertainties. 

The inevitable truncation of model spaces, e.g. Op or- 
bitals for A = 11, 12 and ls-Od orbitals for A = 23, 28, is 
another source of model dependence in Micr • V/TWio • a 



The contributions arising from exch ange currents in al - 
lowed tr ansitions have been studie d by Fukui et al. (1987 ) 
for 12 C, Ciechanowicz et al. (1998) for 28 Si, and Koshigiri 
et al. (1997| ) for ^Na. For i2 C [Fukui et al. (1987[ ) found 
corrections of —4% to Mio ■ cr, +10% to M 12 • cr, and 
+41% to Aficr • V. Recall the large renormalization of 
the axial charge operator Mi<r • V arises from large soft- 
pion contributions in the axial current's time component. 
This renormalization is supported by experimental data 
on hrst-forbidden /?-decay. 

The assumed form for the radial dependence of the 
nuclear wave functions is another source of model depen- 
dence, e.g. at the surface of the nucleus the difference 
in harmonic oscillator and Wood-Saxon wave functions 
are large. Such effects were investigated for 12 C, 23 Na 
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and 28 Si by [Kortelainen et al. (200"(i| ). They found the 
sensitivity of X to the radial form of the nuclear wave 
function was typically 5% or less. 

Lastly we consider the determination of g p /g a from the 
capture rate of the 16 O(0+,0) -> 16 N((T, 120) transition. 
Not e that this requires knowing the absolute values of 



M [ <t and M\<j- V, i.e. not ratios like Micr-'V /M w -cr. 



Also the capture rate is highly sensitive to the 2p-2h 
4p-4h structure of the 16 O ground state and the 2-body 
contributions d ue to soft-pion exchang e. C onsequently, 
as d iscussed by Warburton et al. (1994 ) and Haxton and 
Joh nson (1990 ), the extraction of g p /g a from 16 O( + , 0) 
— > l b N(0~, 120) is a formidable challenge. Indeed War- 
burt on et al. (19*94] ) have cautioned their result for g p /g a 
is highly model dependent. 



3. Conclusions and outlook for g p /g a 

The results in Table |x| are largely consistent with 
PCAC and g p /g a — 6.42. Specifically the values ob- 
tained from 11 B, 16 and 23 Na all support the predic- 
tion of PCAC. The situation for 12 C is more borderline, 
with theory and experiment just under 2a apart, but not 
enough to cause concern. In contrast however for 28 Si 
the experiment determination and theoretical prediction 
are in obvious disagreement. 

Ignoring the puzzle of 28 Si for now, we believe the the- 
oretical uncertainties in extracting g p /g a from these data 
are likely about ±2 or so. For example for recoil polar- 
izations and 7-ray correlations the largest contribution 
to model uncertainties arises via corrections at the 30- 
40% level from the M\<r ■ V operator. Although the 
matrix elements for this gradient operator are rather dif- 
ficult to compute accurately, we note even a 50% un- 
certainty in M\cr ■ V will produce only a 10-20% un- 
certainty in extracting g p /g a - A similar situation arises 
in extracting the coupling g p /g a from hyperfine depen- 
dences, although here the additional contributions from 
A J = 2 + multipoles may increase somewhat the model 
uncertainty. However the model uncertainty is probably 
larger for 16 O(0+,0) -> 16 N(0~, 120), since absolute val- 
ues and not ratios of matrix elements are needed, and 
those required are quite difficult to accurately compute. 
See Sec. |IX.G.5| for details. 

The value of g p /g a = 1-0±};| fr° m 28 Si is rather puz- 
zling. Interestingly the experimental results for 12 C of 
X = 0.26 ± 0.06 and 28 Si of X = 0.55 ± 0.04 are 
quite different, while the theoretical predictions for 12 C 
and 28 Si are not, i.e. the calculated corrections arising 
from axial charge and 2 nd -forbidden operators are sim- 
ilar in 12 C and 28 Si. Inspection of the one-body tran- 
sition densities for 12 C and 28 Si does indicate a differ- 
ence in the two cases. Whereas for 12 C the transition 
is mainly Ops/2 — > Opi/2 with some Opi/2 — > OP3/2, f° r 
28 Si the transition involves numerous single particle tran- 
sitions of similar magnitudes. Therefore sensitivity to 
the nuclear model for the 28 Si(0+,0) ^ 28 A1(1+, 2201) 



transition is likely larger. In addition we observe that 
the agreement of theory and experiment is good for the 
12 C(0+,0) ^ 12 B(l+,0) capture rate but poor for the 
28 Si(0+,0) -+ 28 A1(1+,2201) capture rate, indicating the 
GT matrix element is reproduced well for 12 C but repro- 
duced poorly for 28 Si. For further details see |Gorringe ei 
al. (1999Q 



Finally what new experimental and theoretical work 
is worthwhile? For partial transitions on complex nu- 
clei a limiting factor is the small number of the available 
transitions, i.e. four allowed transitions and one first- 
forbidden transition. Given the unavoidable sensitivity 
to nuclear structure a larger data-set of partial transi- 
tions would be helpful. With more data an improved un- 
derstanding of contributions from M\cr ■ V and M12 ■ cr 
is presumably possible. Unfortunately of course the ex- 
perimental techniques for measuring recoil polarizations, 
gamma-ray correlations and hyperfine dependences are 
often limited to special cases. However new experimen- 
tal studies of muon capture on 32 S and 35 CI are under 
way. Additionally a measurement of the hyperfine de- 
pendence in the 6 Li(l + ,0) — * 6 He(0 + ,0) reaction, where 
accurate nuclear wavefunctions and related nuclear data 
are available, is definitely interesting. 



X. INCLUSIVE RMC ON COMPLEX NUCLEI 

Several factors have motivated investigations of inclu- 
sive RMC on complex nuclei. First the radiative rate 
on complex nuclei is highly sensitive to g p . Second 
the branching ratios for nuclear RMC are comparatively 
large, e.g. the 12 C rate is about 100 times the 1 H rate 
and the 40 Ca rate is about 1000 times the 1 H rate. Third 
the early theoretical studies implied the ratio of radiative 
capture to ordinary capture was only mildly model de- 
pendent. 

We note comprehensive reviews that cover nu clear 
RMC were published by [Mukhopadhyay (1977 ) and 
Gmitro and Truol (1987). These authors have discussed 



in detail the formalism and methods for RMC calcula- 
tions on complex nuclei. Herein we simply outline the 
major approaches, referring the reader to these reviews 
for more detail, and focus mainly on recent developm ents. 
We update the status of model calculations in Sec. X.A 



and experimental data in Sec. X.B. We discuss the inter- 



pretation of the RMC branching ratio data in Sec. X.C 
and describe the situa tion for the RMC photon asymme- 
try data in Sec. X.D. 



A. Theory of nuclear RMC 

The t heory of RMC in nuclei has a long history, dating 
back to Rood and Tolhock (1965 ). Although they were 



not the first to consider the problem, they carefully laid 
out the approach which has become the standard in sub- 
sequent calculations. The standard approach to RMC 
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in nuclei is to develop a non-relativistic Hamiltonian for 
the process which is then evaluated in impulse approxi- 
mation (IA) between nuclear states. There are thus two 
mai n ingredients, the Hamiltonian and nuclear structure. 



sum rules. This approach wa s considered for example by 
Sloboda and Fearing (1978) and more r ecently applied 
to OMC by Navarro and Krivine (1986). It was devel- 
oped further and us ed more recently for RMC by Roig 



M the usual approach the Hamiltonian is derived from and Navarro (199C| ). While still somewhat phenomeno 



the same five diagrams, Fig. H(a)-(e), used to describe 
capture on the nucleoli. This Hamiltonian, which orig- 
inates as a relativistic amplitude, is then expanded in 
powers of the nucleon momentum either directly, or us- 
ing a Foldy-Wouthuysen procedure. The leading order 
is independent of the nucleon momentum and the linear 
order correction is typically 10-20% and is neglected in 
ma ny calculations. The correction to the non-relativistic 



logical, and still depending to some extent on the closure 
approximation, this approach is much less sensitive to 
average excitation energies than the closure approxima- 
tion. 

In principle the best method is to calculate explicitly 
transitions to all possible excited states and sum the re- 
sults. This takes the energy dependence into account 



prop erly, which can be important (Fearing and Welsh 



Har jiiltoman from the terms quadratic in the nucleon mo- 1992j ) . One must always truncate the sum somewhere 



mentum was considered by 31oboda and Fearing (198C) 
and found to be very small. 

This single particle Hamiltonian is then summed over 
all nucleons and matrix elements of the result are taken 

hpt-syppn mirlpar Qt.at.pc; Tn parly ralfiilat.inns vpry rrnHn 



however, which in principle introduces errors. Modern 
shell model codes are good enough however to produce 
enough of the excited states so that for at least some 
nuclei the important transitions can be calculated and 



summed. Most recent calcula t ions (Eramzhyan et al., 



nncllear states wprp hspH , p. g a. simple Sla.tpr Hptermina.nt 
of harmonic oscillator wave functions. In recent calcula- 
tions however more modern shell model wave functions 
have been used which are derived using realistic interac- 
tions. 

The muon deposits a lot of energy in the nucleus and 
thus there can be many levels excited in the final nucleus. 
However all measurements so far are of the inclusive rate, 
and so a technique must be developed to sum over all final 
states. In early calculations closure was used. However 
this introduces a new parameter k max corresponding to 
the average maximum photon energy, or equivalently the 
average nuclear excitation. Unfortunately the rate is just 
as sensitive to k max as to g p . Furthermore, as pointed out 
by Christillin (198l| ) , the closure sum includes many ex- 
cited states which are not allowed in the radiative process 
due to energy conservation. 

There have been several approaches to attempt to by- 
pass this problem. One approach involves obtaini ng the 
spe ctrum of excited s tates from some other source. Foldy 
and 



1998; Gmitro et al, 1990, 1986, 1991) have used this ap- 



Walocka (196d| ) did thio for OMC, by recognizing 



that) much of the strength wont to the giant dipolo reso 



proach. 

One can also attempt to modify or improve the ba- 
sic Hamiltonian. Recall th at the standard Hamiltonian 
( [Rood and Tolhoek , 1961 ) is a single particle operator 
coming from the same basic five diagrams, Fig. ^(a)-(e), 
used for RMC on the proton. In principle one should 
include various meson exchange corrections, which lead 
to two body operators, in the same fas hion as has been 
done for OMC in light nuclei (See Sec. |XT\3|). To our 
knowledge this has not been done, at least recently, for 
heavy nuclei. 

An alternative approach is to look at these meson ex- 
change corrections as effects of the intermediate pion 
rescattering in the nuclear medium. This leads to a 
renormalization of the effective couplings in the nuclear 
medium. This approach has a long an involved history, 
and applies to a number of processes involving axial cur- 
rent matrix elements. It is a bit outside the scope of this 
review however. The reader interested in pursuing this 
further can c ons ider the recent pap ers of Kirchbach and 
Riska (1994|) or|Kolbe et al. (2000 ) or older studies such 



nance state and that one could get this strength from 
empirical photoabsorption cross sections. 2 5 This idea 
was generalize d to RMC by Fearing (1966 ) and further 
generalized by Christillin (1981) who added additional 
phenomenological components to account for transitions 
to quadrupole states. This approach lessened the depen- 
dence on k max but did not totally eliminate it, as there 
were still matrix elements, e. g. the leading p/m correc- 
tions which could not be obtained this way. 

Another approach involved sum rules, in which various 
matrix elements were evaluated using energy weighted 



25 Later calculations extended this to consider SU4 breaking and 
the excitatio n of the spin dipolc s tate by the axial current. For 



example see Cannata et al. (1970) 



as that of Akhmedov et al. (1985). 

Still another attempt to modify the basic Hamiltonian 
has been proposed by Gmitro et al. (1986). The idea 
here is to use current conservation to evaluate parts of 
the matrix element. This is analogous to the Siegert the- 
orem approach which has been used for other processes, 
and originates in the observation that the nuclear ma- 
trix element of the impulse approximation Hamiltonian, 
unlike the nucleon matrix element, does not correspond 
to a conserved current. An attempt is made to fix this 
by expanding the photon field and using the continuity 
equation to eliminate parts of the three-vector current 
in favor of the charge distribution. This leads to some 
different terms and to what the authors call a modified 
impulse approximation (MIA). 

Calculations using this MIA suggest that it is ex- 



tremely important (Eramzhyan et al., 1998; Gmitro et 
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al, 1990, |!986j |1991| ). It reduces the RMC rate by a 
factor of two or more and thus increases the value of g p 
needed to fit a given experimental result rather dramat- 
ically. The approach also seems to produce rates which 
are much less sensitive to g p than the IA. It does how- 
ever suppress usual impulse approximation results so that 
they are in better agreement with phenomenological ap- 
proaches. However there are some caveats. The very 
fact that what is effectively enforcing gauge invariance 
makes such a huge difference is worrisome. The results 
also depend, though not strongly, on an arbitrary choice 
among various ways to make the original expansion of the 
photon field. It also appears that one is using the con- 
tinuity equation on only part of the current, still leaving 
some reference to the three vector current in the prob- 
lem. Clearly this needs to be looked at more carefully 
to determine if it is indeed correct since it makes such a 



larg^ difference in the final results. 



Finally we should mention one other modern calcula- 
tion (Fearing and Welsh , 1992) which uses a relativistic 
Fermi gas model together with relativistic mean field the- 
ory to examine the A and Z dependence of the RMC and 
OMC rates. The Fermi gas model was used in very early 
calculations, but is clearly too crude a model to give good 
results for specific details for individual nuclei. It does 
however allow one to elucidate general trends and select 
out features which are important. This will be discussed 
further in the next section. 



B. Measurement of nuclear RMC 



The first observation of radiative capture was made at 
CEIjiN by |Conforto et al. (1962|) — Muons were stopped in 



Fe alnH photons wprp r\p±pctpA hy gamma-ray r-onvnrsion 

in a iron sheet/spark chamber sandwich and energy de- 
position in a large volume Nal crystal. After subtraction 
of backgrounds a total of five photons from RMC on Fe 
were identified and yielded a Fe RMC branching ratio of 
roughly 1CP 4 . During the following ten years a few more 



studies of nuclear RMC were conducted by Convcrsi ei 
al. ( 1964D, |Chu et al. (1965] ) and Roscnstcin and Ham- 
mer nan (1973[ ). However the measurements remained 



extremely difficult because of the low signal rate and the 
high background rates. In particular under-counted neu- 
tron backgrounds most likely corrupted these early ex- 
periments. 

More recently the availability of higher quality muon 
beams and higher performance pair spectrometers has 
dramatically improved the experimental situation. The 



new era for nuclear RMC was pioneered by Hart et a/. 
(1910) for RMC on Ca. Later experimental programs at 
PSI in the late 1980's and TRIUMF in the early 1990's 
have produced an extensive body of nuclear RMC data. 
Today's beams and detection systems allow data collec- 
tion on Z > 12 targets at count rates of ~ 1000 RMC 
photons/day with nearly background-free conditions. 
In principle the method for determining the branch- 



ing ratio for nuclear RMC is straightforward. It sim- 
ply requires counting the incoming muons and outgo- 
ing photons and determining their dete ction e fficiencies. 
However as discussed in detail in Sec. V.C.4 the small 



branching ratio means troublesome 7-ray backgrounds 
arise from fi decay in the target material, radiative fj, cap- 
ture in the neighboring materials, and the pion contam- 
ination in the muon beam. Additionally neutron back- 
grounds from ordinary capture in the target material and 
other sources on the accelerator site are dangerous if the 
n/7 discrimination is not highly effective. 

Two basic types of pair spectrometers have been widely 
employed in the RMC studies of recent years. In one ap- 
proach, as used by Frischknccht et al. (1988), the pho- 
tons are converted in a relatively thin passive converter 
and the e + e _ pair is tracked in a multiwire ch amber ar- 
rangemen t. In another approach, as used by Dobeli ei 



al. (1988 ) , the photons are converted in a relatively thick 



active converter and the e + e _ pair is measured by com- 
binations of Cerenkov detectors and Nal crystals. The 
former approach offers good resolution, typically 1-2%, 
but at the expense of a low efficiency, typically 10~ 5 . 
The latter approach offers high efficiency, typically 0.5%, 
but at the expense of a poor resolution, typically 20%. 
Such experimental set-ups have yielded excellent n/7 dis- 
crimination and provided nearly background-free RMC 
spectra. 

Most recently a novel large solid- angle pair spectrom - 
eter was developed at TRIUMF by |Wright et al. (1992|) 



for RMC on 1 H. The spectrometer offers both a relatively 
high detection efficiency and a relatively good energy res- 
olution. It has permitted quick and straightforward mea- 
surements of nuclear RMC for numerou s targets, and sig- 



nificantly extended the RMC data-set (Armstrong et al. 



1992; Bcrgbusch et al., 1999; Gorringe et al, 1998) 



C. Interpretation of nuclear RMC 

The world data-set 26 for nuclear RMC, consisting of 
targe ts from carbon to bismuth, is summarized in Table 
XIII| in which we tabulate the ratio i? 7 of the E 1 > 57 



MeV partial radiative rate to the total ordinary rate. In 
addition we plot the data versus atomic number Z in 
Fig. H and versus neutron excess a — (A — 2Z)/Z in 
Fig. H. The figures illustrate some intriguing trends, 
i.e. that i? 7 is observed to decrease from ~ 2 to ~ 0.6 
with increasing Z and increasing a. 27 Note the overall 
trend is somewhat smoother with neutron excess than 
atomic number. Specifically, the isotope effect in the 
mass 58 > 60 > 62 Ni isotopes and the odd/even-A effect in the 



26 We omjt 



he early experiments of ponvcrsi e.1 
al. (1964), |Chu et al. (196j), and 



(1975|) due to large neutron backgrounds. 



Rosenstein and Hammerman 



z ' Note that in this section we will quote all values of ij 7 in units 
of 10" 5 . 
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FIG. 8 The world data for versus Z on A > 3 nuclei. 
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is often difficult to quantify the effect of the various ap- 
proximations and assumptions and so to set limits on the 
uncertainty in the theoretical predictions. 



The phenomenological calculations of Christillin 



(1981), based on the Giant Dipole Resonance model 



'caring , 1966 ; Foldy and Walecka , 1964 ) with parame- 
ters determined from electromagnetic data but including 
effects of a quadrupole resonance with parameters fit to 
the OMC rate, gave R 7 = 2.4 using the canonical value 
g p /g a — 6.7. The calculation of Roig and Navarro (199C ) 
used fairly realistic RPA wave functions but with the phe- 
nomenological aspects of a sum rule calculation to obtain 
= 1.87. The microscopic shell model calculation of 



Gmitro et al. (198E) used simple lp-lh wave functions, 



with some corrections to obtain i? 7 = 4.25 in the stan- 
dard IA and i? 7 = 2.28 in the MIA. Thus the MIA makes 
a large difference, and brings the results in closer agree- 
ment with the more phenomenological approaches. Note 
that all of these calculations, except the standard IA ap- 
proach, apparently neglected the velocity or p/m terms 
which could be up to 10-20%. 

On a superficial level all of these results, except per- 
haps for the standard IA result, are reasonably consistent 
with the experimental result. However if one turns the 
question around and asks what limits are set on g p j g a 
things become much less clear. A theoretical value of i? 7 
less than the experiment implies a larger value of g p /g a 
is needed to fit the data. T hus the sum rule calculation 
of Roig and Navarro (1990 ) requires g p /g a ~ 8. On the 



other hand the phenomenological c alculations of Chris- 
tillin (198l] ) and the MIA results of |Gmitro et al. (1986j ) 
require g p /g a ~ 4 — 5 and the IA result implies that g p /g a 
is much smaller. Thus one has to conclude that it is un- 
reasonable to claim that g p /g a ~ 6.7 is well-established 
by RMC on 40 Ca. 

For 16 O a number of measurements and calculations 
are also available. Unfortunately measurement of RMC 
on 16 O is more difficult, since the 7-ray yield per \i~ stop 
is - 2 x 10~ 5 for 40 Ca and - 0.4 x 10~ 5 for 16 0, and most 
probably the earlier experiments have under-subtracted 
contamination from 7-ray b ackgrounds. Therefore we 
employ the recent result from Bergbusch et al. (1999 ) of 



FIG. 9 The world data for i? 7 versus a 
A > 3 nuclei. 



Al-Si, Ca-Ti pairs fit the a-dependence but not the Z- 
dep endence. Below we discuss the fl-y-data in the context 



= (A — 2Z)/Z on ^7 = 1-67±0.18. Like 40 Ca, both phenomenological cal- 
culations and microscopic calculations are available, but 
similarly the sensitivities of results to approximations are 
difficult to estimate. Taking g p /g a — 6.7 the phenom eno- 
logical calculation of Christillin and Gmitro (1985) gave 



R.~, = 2.08, the shell model calculation of Gmitro et al. 



of t he determination of the coupling g v . 



For u Ca the value of i? 7 is well-established by nu- 
mero us exp eriments. The average of those included in 
Table pCHI| is R 1 = 2.08 ±0.11. A number of calculations 
are available including the phenomenological models of 
Fearing (1966) and Christillin (1981), the mi croscopic 
shell model calculations of IGmitro et al. ( 198q) and the 



(1986|) gave i? 7 = 1.61 in MIA and i? 7 _ 3.10 in IA 



and sum rule calculation of Roig and Navarro (199C) 
gave i? 7 = 1.73. Again, with the exception of the IA 
result, these a re reasonably consisten t with the experi- 
mental result ([Bergbusch et al, 1999|) i? 7 = 1.67 ± 0.18 
with the sum rule and MIA approaches implying r oughly 



the canonical value of g p /g a and the approach of Chris 



fortunately as discussed earlier the calculation of inclu- 
sive RMC on complex nuclei is notoriously difficult. It 



sum rule calculations of Roig and Navarro (1990|). Un- tillin and Gmitro (1985 ) requiring a somewhat smaller 



value. The IA result would require a significantly smaller 
value to fit the data. However, like RMC on 40 Ca, due 
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TABLE XIII Summary of the world data on the quantity 7? 7 
on A > 3 targets, where i? 7 is the ratio of the RMC rate with 
E-y > 57 MeV to the OMC rate in units of 10" 5 . In order to 
assist comparisons in most cases the quoted results are for the 
closure approximation spectra shape where the corresponding 
value of the parameter k max is given in column four. The 



quar tity a is the neutron excess {A — 2Z) /Z. 



Target 



R 7 ( xlO" 5 ) k max ( MeV ) 



Bergbusch et al. (1999) 



l 6 o 


0.000 


1.67 


± 0.18 


88.4 


± 


2.3 


?IA1 


0.077 


1.43 


± 0.11 


90.1 


± 


1.8 


?ISi 


0.000 


2.09 


± 0.20 


89.4 


± 


1.8 


22 -LI 


0.173 


1.30 


± 0.12 


89.2 


± 


2.0 


natry 
40 ^ r 


0.280 


1.31 


± 0.15 


89.2 


± 


3.4 


? 7 ai Ag 


0.296 


1.12 


± 0.13 


89.0 


± 


3.2 



Gorringe et al. (1998|) 



Ni 






071 


1 


48 


± 





08 


92.0 


± 


2.0 





143 


1 


39 


± 





09 


90.0 


± 


2.0 





214 


1 


05 


± 





06 


89.0 


± 


2.0 



Armstrong et al. (1992) 



14»1 
20^ a 

J 2 a 'Mo 

natcr 
50 H 

rtat- 
82 
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0.077 1.43 ± 0.13 
0.000 1.93 ± 0.18 
0.000 2.09 ± 0.19 
0.283 1.11 ± 0.11 
0.3/4 0.98 ± U.Uy 



90.0 ± 2.0 
92.0 ± 2.0 
93.0 ± 2.0 
90.0 ± 2.0 
87.U ± 2.0 



Armstrong et al. (1991) 

12,"! 

6 ^ 

l 6 o 



Frischknecht et al. (1988) 



»0 



Dobeli et al. (1986) 



2 C 
3 



0.527 0.60 ± 0.07 — 84.0 ± 3.0 

0.000 1.98 ± 0.20 
0.000 2.18 ± 0.20 
0.000 2.04 ± 0.14 



0.000 3.80 ± 0.40 

0.000 2.70 ± 1.80 
0.000 2.44 ± 0.47 



89.9 ± 5.0 



13 1 

ISCa 
ae'Fe 

165 T 

67 

209 1 

83 



Ho 
Bl_ 



0.077 1.83 ± 0.26 
0.000 2.30 ± 0.21 
0.146 1.71 ± 0.17 
0.463 0.75 ± 0.13 
CL518 nfi? + nns 



88.8 ± 1.8 
92.5 ± 0.7 
90.2 ± 1.1 
84.1 ± 5.1 
ss 9 + n f\ 



to the approximations in the calculations and the differ- 
ences among the results of different calculations the case 
for g p /g a ~ 6.7 in RMC on 16 is not firmly established. 

There are also both theoretical and experimental re- 
sults for 12 C. The most recent result, and the one with 



Frischknecht et al. (1985} 



Hart et al. (1977) 



tgCa 



0.000 1.92 ± 0.20 90.8 ± 0.9 



0.000 2.11 ± 0.14 86.5 ± 1.9 



by far th e smallest uncertainty is that of Armstrong ei 



al. (1991 



1.98 ± 0.20. 



This is to be compared 
with the sum rule result (Roig and Navarro , 1990) of 



Ry = 1.42 
result 



and the I A result, i? 7 



R 



1.48 flGmitro et al. 



= 3.60 , and the MIA 
1990| ). In this case 



both MIA and sum rule results are too low, implying 
9pI ' 9a ~ 10 — 13, to fit the data whereas the IA is too 
high, implying a very low value of g p /g a - 

Finally we should mention t he calculations for 
58 ' 60 ' 62 Ni of Eramzhyan et al. (1998| ) carried out in a mi- 
croscopic model using the quasiparticle RPA. Again the 
IA results for i? 7 are much higher than those in MIA, 
but even the M IA results arc signific antly higher than 
the experiment (Gorringe et al., 1998). 

Thus by now there have been a number of experiments 
and enough calculations that we can make comparisons 
of i? 7 for several different nuclei. The situation can only 
be described as confused. The standard impulse approx- 
imation calculations are consistently too high, implying 
a value of g p /g a much smaller than the canonical value. 
The MIA start out too low for 12 C, implying a large 
value of g p /g a , and rise with increasing A to become sig- 
nificantly too high for the Ni isotopes, implying there 
a small value of g v /g a . T he ph cnomcnological calcula- 
tions of Christillin (1981 ) and Christillin and Gmitro 



(1985) for 15 O and 4U Ca are both too high and the sum 
rule calculations show no consistent pattern. It appears 
clear that we are not yet at a stage where theoretical un- 
certainties are sufficiently under control to consistently 
reproduce experimental results. 

In addition to discussing the detailed calculations for 
specific nuclei it is worthwhile to consider the systematics 
of RMC data versus atomic number and neutron excess. 
Such systematics were e xamined in the non-re lativistic 
Fermi-gas calculation of Christillin et al. (198C ) and the 
rclativ istic Fermi-gas calculation of Fearing and Welsh 
(1992). These models clearly oversimplify the nuclear 



structure, for example omitting the important effects of 
giant resonances in muon capture, or the effects of shell 
closures. However they demonstrate a number of inter- 
esting dependences of muon capture on Z and a. 



For concretenes s we consider the calculation of Fearing 
and Welsh (1992] ) . These authors have carefully studied 



the dependence of inclusive OMC and inclusive RMC on 
the input parameters and the model assumptions. They 
stress that the OMC rate, the RMC rate, and their ra- 
tio, are highly sensitive to phase-space effects, i.e. things 
that alter the available energy for the neutrino and the 
photon. For example including the /x~ atomic binding 
energy, which increases from ~ 0.1 MeV in light nuclei 
to ~ 10 MeV in heavy nuclei, decreases the OMC rate by 
a factor of two and the RMC rate by a factor of eight for 
the heaviest nuclei. Other parameters, e.g. for Fermi-gas 
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models the Coulomb energy, symmetry energy, etc., also 
have large effects on the available energy and therefore 
the rates. Consequently the authors caution that reli- 
ably extracting the coupling g p from inclusive RMC on 
complex nuclei is difficult. 



D espite suc h concerns the calculation of [Fearing and 
Wei ;h (1992] ) does reproduce the overall dependence of 



Rj data with Z and a. This model-data agreement how- 
ever suggests no reason to invoke a large-scale medium- 
modification of g p /g a to explain the Z-dependence of i? 7 
data. 

Earlier suggestions for an A-dependence of the cou- 
pling g Vl as discus sed by Gmitro and Truol (1987) and 
Dobeli et al. (1988), were largely grounded in large values 
of g P /g a obtained from shell model calculations on light 
nuclei and small values of g p /g a obtained from Fermi- 
gas calculations on heavy nuclei. Given the different sys- 
tematics of the various models and the difficulty of all 
models in consistently fitting the data, such differences 
do not provide real evidence for an A-dependence of the 
coupling g p . 

In summary, although the RMC rate on complex nuclei 
is undoubtedly g p -dependent, the problem of separating 
the coupling constant and nuclear structure is extraordi- 



nari ly tricky. While the overall features of nuclear RMC 
are generally accommodated by g p /g a = 6.7 there is a 
model uncertainty of at least 50%. Thus RMC in heavy 
nuclei is not yet competitive with RMC or OMC in hy- 
drogen or light nuclei as a source of information on g p /g a - 
For inclusive RMC on complex nuclei we find the 
progress on model calculations is lagging the progress on 
experimental data, and a breakthrough on the theoreti- 
cal side is needed before this particular tool for studying 
the coupling is competitive. However, exclusive RMC on 
complex nuclei, which is so far unmeasured, might offer 
interesting possibilities. 

D. Comments on photon asymmetries in nuclear RMC 

In RMC the gamma-ray direction and muon polariza- 
tionl are correlated according to 



N(6) = N{1 ±aP At cos( 



(50) 



where P M is the muon polarization, a is the asymmetry 
coefficient, and 9 is the angle between the /z-spin axis 
and the 7-ray momentum axis. The left handedness of 
the neutrino leads to a ~ +1, and g v is m anifest as de- 
partures of a from unity ( Fearing , 1975 ). The specific 
depende nce of a on g v /g a for RMC on 40 Ca was com- 



puted by|Rood and Tolhoek (1965[), |Gmitro et al. (198l[) 



Christillin (198l| ) and pmitro et al. (1987j ). 

To measure the asymmetry one stops the incoming 
muons in a suitable target, precesses the muon spin in 
a magnetic field, and measures the time spectrum of the 
outgoing photons. The resulting time spectrum consists 
of an exponential decay with a sinusoidal modulation, 
where the amplitude of the sine-wave is governed by the 



TABLE XIV Summary of the world data on the photon 
asymmetry in radiative capture on calcium. 



Ref. 



asymmetry a 



Hart et al. (1977) 



Dobeli et al. (1986) 



Virtue et al. (199C| ) 
world average 



0.90 ±0.50 
0.90 ± 0.43 



-, 09+0.54 
J -° £ — 0.47 



1.02 ±0.25 



product P p a. Note that the time spectrum of the Michel 
electrons is employed to measure P u and isolate a. In 



the first measurement by di Leila et al. (1971) the au- 
thors used a single Nal crystal for 7-ray detection. In 
the later meas urements by |Hart et al. (1977 ), Dobeli et 
al. (1986 ) and Virtue et al. (1990| ) the authors used a 
separate 7-ray converter and Nal calorimeter. A 40 Ca 
target was used in each experiment. 

The experimental difficulties originate from the tiny 
RMC branching ratio and the small fi residual polariza- 
tion. Consequently the backgrounds are severe. They 
include 7-rays from pion capture in the target, 7- 
rays from Michel bremsstrahlung in the target, neutrons 
from muon capture in the target, and n/7 backgrounds 
from cosmic-ray and accelerator sources. Note that a 
prompt cut reduces the 7-rays following 7r-capture and 
a E > 57 MeV cut reduces the bremsstrahlung follow- 
ing /j,-decay. The neutron background is reduced via 
the converter-calorimeter set-up and cosmic-ray back- 
ground is reduced via combined passive and active shield- 
ing. 



In Table XIV we summa rize the results of Hart et al. 
(1977| ), pdbcli et al. (1986| ) and|Virtuc et al. (199C| ). We 



omit the earliest results of |di Leila et al. (1971| ), which 
suffered severe neutron backgrounds. The various exper- 
iments are mutually consistent and yield a world average 
value a = 1.02 ±0.25. 

Unfortunately the effect of g p on a is (i) not large and 
(ii) model dependent. For exa mple with g D /g a ~ 6.7 
the phenomenological model of Christillin (1981) give s 
a ~ 0.80 and the shell model of Gmitro et al. (1987 ) 
gives a ~ 0.90. Conservatively the measurements yield 
g p /g a < 20. 



XI. SUMMARY 

The determination of g p is important for several rea- 
sons. First, whereas the values of the proton's other 
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weak couplings are nowadays well determined, the in- 
duced pseudoscalar coupling is still poorly determined. 
Second, based on chiral symmetry arguments a solid the- 
oretical prediction for g p with 2-3% accuracy is avail- 
able. Third, the prediction is founded on some elemen- 
tary properties of the strong interaction, and determining 
the coupling is therefore an important test of quantum 
chromodynamics at low energies. 



Since the classic review of Mukhopadhyay (1977) an 
extensive body of experimental data on the coupling g p 
has been accumulated. This work spans the elementary 
processes of ordinary muon capture and radiative muon 
capture on hydrogen, muon capture on few-body sys- 
tems, and exclusive OMC and inclusive RMC on complex 
nuclei. The experimental approaches have ranged from 
ultra-high precision measurements to ultra-rare process 
mea surements, and include some novel studies of spin 



phe nomena in complex nuclei. The disentangling of the 
coupling constant from the physical observables has in- 
volved diverse fields from ^t-molecular chemistry to tra- 
ditional nuclear structure and exchange current effects. 

One would expect muon capture on hydrogen to be the 
most straightforward and most easily interpreted of the 
muon capture reactions. This is the situation from the 
nuclear perspective, but unfortunately there exist atomic 
and molecular complications which also must be thor- 
oughly understood. For OMC there are a several older 
experiments tog ether with the most recent and precise 
measurement of Bardin et al. (1981a) whereas for RMC 



there is only the TRIUMF experiment (Wright et al. 
199c]]). To extract g p from these results we updated the 
theoretical calculations to include the current best values 



of t le other weak couplings and their q dependences 



Ad ditional ly we updated the capture rate of Bardin et 
al. jl981a ) for the present world-average value of the 
positive muon lifetime. We found, for the standard val- 
ues of the n chemistry parameters and specifically with 
values of g„ = 12.2 ±1.1 from 



A op = 4.1 x 10 



New work on fi capture in liquid H2 and gaseous H2 
is in progress. At TRIUMF an investigation of the \x 
chemistry in liquid hydrogen has taken data, and at PSI 
a measurement of fi~ lifetime in gaseous hydrogen is now 
under way. Hopefully these experiments will help to clar- 
ify the situation on g p . Looking further ahead, perhaps 
both new facilities, e.g. an intense muon source at a neu- 
trino factory, and new techniques, e.g. neutron polariza- 
tions or hyperfine effects, will permit a determination of 
g p to an accuracy of 2-3% or better. 

The situation for muon capture on pure deuterium is 
inconclusive. Recall that in deuterium the g p sensitivity 
is smaller than in hydrogen, the neutron experiments are 
harder, and nuclear models are needed. At present the 
world data for mu on capture in pure deut erium comprises 
the experiment of Cargnclli et al. (1989), which appears 
consistent with PCAC, and the experiment of p3ardin ei 



(1986 ), which may be inc onsistent with PCA C. Also 



the results on fid capture by Bertin et al. (1973) using a 



H2/D2 target are even more puzzling. 

Further studies of fid capture are clearly worthwhile. 
Obviously resolving the possible discrepancy between 
the doublet rates obtai ned by Cargnclli et al. (198S| ) 
and Bardin et al. (1986| ) is important. Additionally we 
note the chemistry of muons in pure H2 and pure D2 is 
quite different in several ways, including a slower rate 
for hyperfine depopulation and near absence of muonic 
molecules in pure D2. Such features may permit the 
study of g p via alternative approaches, such the hyperfine 
dependence of the capture reaction. 

In /^+ 3 Hc — > 3 H±z; capture, the recent pr ecision mea- 
suremcnt of the statistical capture rate by |Ackerbauer 
et al. (1998[ ) and ground-bre aking measurement of the 



the TRIUMF RMC experiment ([Wright et al. 1998|) 



10.6±2.7 from the Saclay OMC experiment (|Bardin 



recoil angular correlation by |5oudcr et al. (1998 ), were 
major achievements. Preliminary data on the radiative 
captu re rate in the 3 He — > 3 H channel are also avail- 
able ( Wright et al, 200C ). Further modern treatments 
of A — 3 wave functions and 2-body exchange cur- 
rents have been applied to the process by Conglcton and 
Truhlik (1996| ) and others. The OMC results for g p 



et al, 1981a), and g p — 10.5 ± 1.8 for the world aver- 
age of all OMC experiments. These updates increased 
the g p extracted from OMC, as compared to the original 
analysis, so that now both OMC and RMC give results 
larger than predicted by PCAC and in fact agree better 
with each other than with PCAC. When uncertainties 
are taken into account however, only the TRIUMF re- 
sult is clearly inconsistent with PCAC, while the OMC 
results are only marginally inconsistent with PCAC. 

Our strong prejudice that PCAC is not violated at such 
a level, makes the situation for fi capture in hydrogen 
very puzzling. We therefore have examined some sugges- 
tions for solving the puzzle, such as modifications to the 
fi chemistry and the role of the A resonance. By 'tuning 
parameters' the discrepancy between the hydrogen data 
and the PCAC prediction can be reduced somewhat, but 
no clear-cut solution, which makes OMC and RMC si- 
multaneously agree with PCAC, has emerged. 



are completely consistent with PCAC, and the value of 
g p = 8.53±1.54 from Ackerbaucr et al. (1998) is probably 
the best individual determination of g p . 

The extraction of the coupling from the measurement 
of the fi+ 3 He — > 3 H+v statistical rate is unfortunately 
limited by theoretical uncertainties in calculating con- 
tributions from exchange currents. Further theoretical 
work is definitely worthwhile to quantify the contribu- 
tions arising from radiative corrections and nail-down the 
uncertainties arising from exchange currents, etc, but im- 
provements in extracting the coupling may be difficult. 
A precision measurement of the recoil correlation, which 
has enhanced sensitivity to g p and reduced sensitivity to 
exchange currents, would be extremely interesting. 

A significant body of new data on exclusive transitions 
in ordinary capture on complex nuclei has been collected 
recently. This includes the measurement of 7-ray corre- 
lations in /i 28 Si and hyperfine dependences in /i n B and 
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/z 23 Na. These data complement earlier investigations of 
polarizations in /i 12 C and rates in /i 16 0. Furthermore, 
modern well-tested models of nuclear structure in Op, Is— 
Oci nuclei offer improved multi-particle wave functions for 
interpreting these experiments. In most cases, with the 
exception of 28 Si, the values of g p that are extracted from 
the experiments with such wave functions are consistent 
with PCAC. Generally the major uncertainty in extract- 
ing g p originates in the interplay of the contributions from 
the pseudoscalar coupling, arising from the space part of 
the axial current, and the axial charge, arising from the 
time part of the axial current. Unfortunately it is difficult 
to precisely quantify such theoretical uncertainties. 

Additional experiments on exclusive OMC could be 
helpful to our understanding of the interplay of the in- 
duced pseudoscalar contribution and the axial charge 
contribution. Clearly if a wealth of data were available 
for nuclear OMC a better assessment of model uncertain- 
ties would be possible. However, the experiments gener- 
ally involve the use of methods and observables that can- 
not be applied to large numbers of exclusive transitions. 
Alternatively a measurement of observables in a transi- 
tion such as 6 Li(l + , 0) — ► 6 Hc(0 + , 0), where highly accu- 
rate wave functions and related nuclear data are nowa- 
days available, is definitely interesting. 

In inclusive RMC on complex nuclei the application of 
pair spectrometers has yielded data with good statistics 
and little background, and enabled the systematics of 
the RMC rate across the periodic table to be mapped- 
out. Unfortunately the situation in regards to the model 
calculation of the inclusive rate is less satisfactory. In 
general the models employed for inclusive RMC have too 
many assumptions and too many parameters in order to 
reliably extract the coupling g p . Furthermore questions 
remain regarding the effective Hamiltonian for radiative 
capture on complex nucleus. Therefore we believe that 
earlier claims for large renormalizations of g p in inclusive 
RMC on complex nuclei were premature and that, at the 
level of the uncertainties, the data are consistent with 
PCAC. 

In conclusion we hope we have convinced the reader 
of the importance of the coupling g p , both as a funda- 
mental parameter in nucleon weak interactions and as 
an important test of low energy quantum chromodynam- 
ics. In recent years an impressive body of experimental 
data has been accumulated and most data are consistent, 
within sometimes large experimental and theoretical un- 
certainties, with the PCAC prediction for the coupling 
g p . Unfortunately the results from /x capture in hydro- 
gen, which should be the simplest and cleanest process, 
are very puzzling. We believe that the resolution of this 
puzzle should be a high priority, and that, until the sit- 
uation is clarified, the accurate determination of g p , and 
implied testing of low energy QCD, remains an important 
task but an elusive goal. 
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